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method of fan sound mode structure determination 

FINAL report 


bv 

G. F . Pickett, T. G. Sofriti and R. A. Wells 
I.O SUMMARY 

I his report presents a method for the determination of fan sound mode structure in the in- 
let of turhofan engines using in-diiet acoustic pressure measurenier.ts. I he method is mised 
on the simultaneous solution of a set of equations whose unknowns are modal amphtude 
and phase. By deploying microphones etpial in number to the number ol modes to be calcu- 
lated tlie set of simultaneous eiiuations can be solved for the amplitude ami pliase oi each 
mode A computer program tor the solution of the equation set was developed. An addi- 
tional computer program was developed which calculates microphone locations the use oi 
which results in an equation .set that does not give rise to numerical instabilities, burther- 
moie. by calculating a set of inlluencc cocfticicnts for each case considered, the met.iod can 
be u.sed to determine the sensitivity of the modal calculation procedure to errors m the 
knowledge of the measuring location and in the mcasuretl phase angles and pressure amph- 
tudes. Also, an estimate can be obtained of the effect on the modal calculation procedure ol 
extraneous modes not included in the modal calculation. 

An experimental assessment is presented that shows that tlie method is reasonably accurate 
if either the number or the amplitudes of the extraneous modes are not too large. Based on 
this experimental assessment and an orderof-magnitude analysis ol the mnuence coelti- 
cients, a procedure is presented that can be used to estimate the number of microphones 
and the accuracy of experimental measurements that would be required to achieve a speed ic 
accuracy in calculating the amplitudes and phases of the propagating coherent sound modes. 

In addition to the development of a method for determination of coherent modal slructure, 
experimental and analytical approaches are developed for the determination ol the amph- 
tud.' freciuency spectrum of randomly generated sound modes for use m narrow annulus 
ducts. l\vo approaches are defined: one based on the use ol cross-spectral techniques and 
the other based on the use of an array of microphones. 


2,0 INTROnUCTION 


All modem liin'i-bypnss nitin, commerdul lr;mspoit cnpiiies ineojponile noise reduction 
fe;itmes in tlieir ilesipn ;is evidenced hy the hl:ide-;ind-v:nu‘ lUimher selections in tlte liin 
sections these engines. TlU'se selections iilonn with other noise reducint; fealnres such 
as larpe lilade-vatie spacinj', cycle selection, and improved acoustic treatment designs 
have resnltcil in sipnil'icant reductions in aircraft noise compared with that from firsl-nenera' 
lion aircraft. 


Despite liie.se reductions, the aircraft industry is faced with legislative rct|uiremcnts to 
meet yet lower noise limits in the future. Also, because of the equally urgent requirement 
to lower fuel consumption, the impact on efficiency of new noise reducing features must 
he minimi '.ed. 


I-’aii noise is predicted to he tiie dominant noise source for most advanced gas turbine en- 
gine designs. More specifically, fan tone noise at blade passing frequency (BPF) and its 
harmonics continues to be the most prominant spectral noise component from engines 
tested without inlet and aft duct acoustic suppression liners. 

Fan tone noise can be generated by four different sources: 1) rotor wake-stator, or 
stator wake-rotor interaction, 2) rotor-strut, pyioii, etc. potential field interactions, 3) 
rotor-turbulence and steady iiillow distortion interactions, and 4) pressure field of the 
rotor alone. F.xperimental evidence suggests that randomly phased tone noise generated 
during ground static testing is eau.scd by turbulent flows in the fan inlet that have very 
long axial length scales. These turbulent flows are believed to change in llight, resulting 
in greatly reduced turbulence intensities and length .scales. These reductions significantly 
reduce the importance of rotor-turbulence interaction as an in-llighl tone noise source, and 
is currently thouglit to explain the observations that flight fan tone noise levels from iiigii 
bypass ratio turbofan engines are significantly lower than those measured during static 
engine tests. The remaining contribution to llight fan-tonc-noise levels of randomly generated 
noise sources has yet to be determined. Broadband noise is, of course, generated by random 
noise sources and is becoming more of a contributor to overall flight fan noise levels as tone 
noise levels arc reduced. 


The dominant noi.se source in flight, however, is thought to be due to coherent acoustic 
fields generated by the pressure field of the rotor and the interaction of rotor wakes with 
stator vanes ami vice versa. The former will be classified as “mtor locked*’ and is a significant 
source of multiple piire-tone inlet noise from fans operating at supersonic tip speeds. The 
latter noise source will be classified as “interaction tone” noise and can he significant, parti- 
cularly for typical low-bypa.ss-ralio lurhofan engines, even at low subsonic fan tip speeds. 


I■an.s for engines currently designed for conventional takeoff and landing (CTOl.) aircraft 
(l(tr exam|ile, the JTdD) exhibit acoustic iluci moiles generated by all of the above source 
types, These fans typically operate at supersonic tip speeds at takeoff (generating inlet 
railiated multiple pure tone noise), and at about sonic tip speeds during landing approach 
wlierc, at least for llie secoiu’ harmonic of blade passing frequency and liiglier harmonics, 
the interaction tone mechanisms that lead to coherent modal structure are well cut on and 
ilomiiiale over the tone generated by the rotor/lurinilence interaction, which leads to inco 
herciit modal structure. 


PAGiO Tb.Ai'ii- r*- 


Ill nil iK'w lilt) ilcsifiDs, hliKU’-Viinc iiitiTiU'tion Ihc'Diy is cxploilotl to tlif I'Xtvnl possilili* lo 
ii)itiimi/o noiso ;it hliiilo piissinp Iretpioiicy. I'liis llu'oi y, hiiseil on >iii nniilysis tluit is con- 
stiiiilly lu*inp rcliiu’il, iik-ntllk's llic‘ nioik's thnt am pi'opnftiilt', Init it Ims not yet lu’en 
devolopi'il to tile extent lliiit tlie iiinplitiKles iiiiil plnises orciicli of the propiipitlinp inoik’s 
am he preilieteil leliiihly. Aii:ilytie;il models also exist for nindomly generated tone noise 
iiml broadhand noise due lo rotoi-iiinow Unl'idence inleraclion. Again, these motlels 
identify the modes tliat ean prop..gate, hut are not suflleiently developeil lo prediet reliably 
tlie expected anipliliides of the propagating modes, 

I'urtlier tone aiul l^roiulhaiul imise reiiuctions could lie achieved through iniproveil tieoustic 
treatni .‘111 design if the propagating modal structure could he c|uantified. Since reliable 
analytical predictions of modal structure arc not likely lo he available in llie near future, a 
need exists to he able to measure reliably fan modal structures. These measured modal 
structures could then he u.setl directly as input to analytical acoustic liner preiliction sys- 
tems, resulting in improved liner efficiency. The measured modal structure could also be 
used to check the developing theories for fan noise generation and to evaluate noise reduc- 
tion design concepts. Once a valid fan noise generation model on a modal ba.sis has been 
developed, it can be coupled with acoustic liner modal prediction systems .so tliat fan and 
inlet treatments can be jointly optimized to minimize engine weight and performance 
losses while achieving sped 11c fan noise goals. 

The specific objectives of the work describeii in this report, and performed under NASA 
contract (NAS3-20047) at Pratt & Whitney Aircraft are to: 


1 ) Develop and check out experimental and computational methods for the cal- 
culation of the amplitudes and phases of propagating coherent sound moiles in 
the inlet ol a turbolan engine using in-duct acoustic measurements as the cal- 
culation input. 

2) Develop experimental and analytical methods for the defermination of the ex- 
pectcii amplitudes ol propagating randomly generated sound nioiies in the inlet 
ol a turbolan engine using in-duct acoustic measurements as the calculation 
input. 

The report is organized as follows: In section 3.0 the overall program is described. After 
an overview of the program (section 3. 1), an analytical method is presented for calculating 
coherent Ian sound modal .structures (.section 3.2). In section 3.3 an analytical approach 
is presented for calculating random sound fiekis in ducts, and this is followed by an experi- 
mental assessment of the accuracy of the method for determining coliereni fan stnind 
modal structures (3.4). A discussion is tiien presented of the accuracy requirements for 
genei.il application ol the method lor coherent sound fields (.3..^) leading to an analysis of 
(he tnefliod utility (3..S.(,). The overall program results are reviewed in section 4,0 followed 
by a conclusions section (.S.O) anil a section recommending future work ((i.O) 
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„„.|lu,.ls l„r c,.lc„l„liim am|.lil,i,los ,„„l l'li„s« »U„lica„l “'''''77"''?;! , !,„ 

tuiboi'an t-'iiniiu's iisinn iinlud nu’asmvmcnls as llic lakulalinii in) u(. • < 

to lK^ ,|™1 „|h,I r„r dolc„,ii„l„y m. “riK-lna' „r in-i-c.l .U.cl 

ami coherent discrete-treqiiency sound fields. 

SiSS 

Uowever, for ,lc.«,nini„. lUc Hmiclua „f co icaiil 

irhasc is desired. Tliese eom)n.ter programs were deveio)K-d with sullicient storage to 
aecomtnodale a greuip of Hfty coherent inodes. 

A omiv was ,„.rlOr,„e.l lo .Is'llna llic om.r ii. Ilw rakolaad "»"k' ail.l'lil'"los ,m,l l.hassw 
doe to errors c iused by tolerances in tneasiiring microphone locations and maeeuraeies 
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lions wi‘iv iliM-iw’i! lor tlu- iiilliiom c coi'lTiut’iils or iiiirlial ilcrivalivi-s ol lliv inoilal co- 
dTiik’iils with ri-spcft to rirsl-onU'r errors in )>ivssmv mcasiiivnu'iU ami microphone 
location. These iiiriuence coeflicients were incorpoiateri in the secoiui computer program 
anil are evalualeil as a inirt of the computational proccilure I'or each ilata case that is uni, 

Upon multiplication hy a particular incasuremenl error, the corrcspoiulnm crronn the moik 
amplitude and phase can he determined. In addition, errors in the calculated modes were 
assessed due to the el fects of a nuillii>Iicity of e.xtraneous modes that propagate at the 
frequency ofinlerest, hut are not included in the calculation procedure. 

Tliesc two computer programs are written entirely in l'Oim<AN IV and are fully operative 
on the NASA computers. A user's manual was written for each program, 1 hese manuals 
define the algorithm and corresponding suhroulines that were developed to meet the com- 
puter program ohjeclives. Also, the input imraineters and output format are discussed tor 

each option. 

The method and computer programs were assessed in an experimental program approved 
hy the NASA Project Manager to determine the modal structure of the inlet duct souiu 
field in the P&VVA 2.‘5 cm (10 inch) diameter compressor rig. During the course of these 
experiments, measured acoustic pressure signals were obtained lor six dillcrcnt modal 
structures, four modal structures were generated hy the wakes from a 34 rod stator nUer- 
acting with a 32 hladed rotor, for two rotor speeds, different modal structures existed at 
both blade passing and twice blade passing frequency. The final two modal structures 
simulated steady iuHow distortion hy placing a single rod assembly upstream ot the rotor, 
for tills configuration, the rotor was oiieratetl at two speeds amt data were acquired at 
only BPI-. The measured acoustic pressures were filtered and then signal enhanced to pro- 
vide amplitude and phase information at each microphone location. Then the pressure 
signals were used as input to the calculation procedure, yielding amplitudes and phases of 
the specific inodes. 

The method utility and measurement accuracy were assessed hy using the calculated mode 
amplitudes and phases to predict the overall sound field measured by microphones at other 
duct locations not used to provide input to the calculation procedure. This study demon- 
strated that the method could he used within a desired accuracy when both the relative^ amp- 
litudes and the number of those modes not included in the modal calculation were sullT 
ciently small so as to yield errors in the mode amplitude and phase due to only measure- 
ment inaccuracies. The required accuracy of experimental measurements to ensure a de- 
sired accuracy in the modal coefficients is demonstrated to he a lunction ol the rel.itive 
amplitude, axial wave nuinhcr, and circumferential order of the modes. 

Monthly technical reports were completed am! submitted as stipulated hy the program 
schedule for the purpose of reporting work conducted under Contract NAS3-20047. In 
addition, this final report and accompanying user's manual for operation of the computer 
programs ilocuineiit the analytical analysis, experimental evaluation, and computatioii.il 
algorithms that lead lo the accomplisliineiil of the overall progiam olijectives. 
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3.2 APPUOACH FOR MliASURlNG COHERFNT MODAL STRUCTURES 


3.2.1 Uackgrouiul 

The bjisic work upon which fan spinning mode theory is based was initially pre.sented in 
I%1 by Tyler and Sofrin {ref. 1) following extensive analytical and experimental studies. 
Extending this work to include effects of axial How, Sofrin and McCann (ref. 2) showed 
that the general fonn of any coherent acoustic wave (such as those generated by rotor-stator 
interactions) in an infinitely long cylindrical duct can be written as the real part of 


p(r. e, X j w) = .2* ^ A 

' M ^ 







(3.2-1) 


u/Kere 


kr ~ 




|i:y 

I- Mt 


and where notation has been chosen so as to be consistent with reference 2. The notation 
description is also tabulated in Appendix A. In reference 2, it was shown that the number 
of modes that can propagate at a given frequency is bounded by the condition that k^^ must 
be real. If k^ becomes imaginary, it can readily be seen that as x increases, the term eikxX 
produces a term that decays exponentially. 

In fan inlets and fan ci . . t • ivpical of current turbofan engines, the assumptions underlying 
the acoustic theory of rti^icncc 1 are only partially met. In general, the duct geometry 
varies with axial distance, and the ducts are finite in length. For typical engine inlet geo- 
metries, however, these effects arc iirobably small, because the variation in geometry with 
axial length and the finite length effect is only significant for duct modes near cutoff. 

In spite of these difficulties, the analytical consideration of rotor-stator interaction tone 
noise, in terms of the general duct pressure, equation (3.2-1) has proven extremely success- 
ful. b.ach term in the simimation is considered as a separate acoustic mode with amplitude, 
'^mp’ phase, and it has been shown that rotor-stator interaction tone noise gen- 
erates a finite set ot these modes. The advantage of considering noise fields in terms of acou- 
stic duct modes is that the generation, propagation, and radiation to the far field can be 
considered separately lor each mode. The overall pressure at any desired location can then 
be determined by summing the individual mode contributions. The prime coherent noise 
fields that are (or can be) conveniently considered in terms of duct modes are rotor-.stator 
interaction, rotor-steady inflow distortion interaction, and multiple pure tone noise at that 
location in the duct where the shock waves arc sufficiently weak so as to be practically con- 
sidered as linear acoustic waves. 
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Till' iiuitliciiiiitical coiicopt ol iliict iiuhIcs (l(.’scrilHMt uhovc, has iiiotivatcd a luiiiihcr olT*x- 
pciimcntal stiulii's Jo measure modal structures and to verify theoretical results. Mca.sure- 
ment ol the amplitude and phase for low circumferential lohe numbers (m) were made by 
lyler and Sofrin using a rotating microphone on a probe mounted axially in the inlet of a 
10-inch <liamcter model fan rig. These measurements were made by essentially cross-correl- 
iiting the output acoustic pre.ssuro signal with a once-per-blade passing trigger signal obtained 
Irom a proximity pickup over the blade tips. Subsequently, other investigators (ref. 2 & .1) 
measured low m-lobe number modes by cross-correlating while others (ref. 4, .S, (>) u.sed 
cros.s-spectra tcchnitiucs on .signals from pairs of microphones in (he duct. 

A very thorough investigation by Plumblce et al, (ref. 7) on sound propagation in acoustic- 
ally lined ducts witli flow required the generation of low m-lobe number acoustic duct 
modes with an array of speakers and the accurate detection of the amplitudes and pha.scs 
of these modes with in-diict microphones. In a set of well controlled experiments, they 
circumferentially traversed a duct probe in discrete steps and cross correlated the acoustic 
signal output with a reference microphone in the duct. Specific duct modes were thus gen- 
erated and measured with a high degree of accuracy. This experimental method was then 
used in the verification of an analytical model that predicts the acoustic decay due to 
specified acoustic treatment designs on a modal basis. 

The previously cited works dealt with modal structures generated by laboratory equipment 
under carefully controlled conditions. More recently, Kraft and Posey (ref. 8) measured 
the relative amplitudes of sound energy in circumferential modes in a duct upstream of a 
model fan stage under conditions representative of an actual jet engine (i.e., large axial 
flows, large numbers of generated acoustic modes some with large lobe numbers (m). etc). 
They used cross-spectral techniques on the outputs from two microphones, one on a fixed 
probe used as a reference and the other on a circumferentially traversing probe. In all 
tests reported, they kept the pair of microphones at the same fixed radial location. Only 
the relative distributions of sound energy in the circumferential modes were obtained, and 
so the accuracy of their measurements could not be evaluated (e.g., by predicting sound 
field at another duct location). Also, they were not able to relate all the measured modes 
to the expected modal distribution that would be predicted by application of theTyler- 
Sofrin theory to the fan design used for the test. 

Tlie approach for determining coherent fan sound mode structures, adopted by Pratt & 
Whitney Aircraft, differs from the previous approaches described above in two principle 
ways. First, the resultant amplitudes and phases of the coherent sound field in the dud 
is obtained by signal averaging the outputs from a number of fixed microphones in the 
duct. From an assumption of the set of duct modes considered to contain the bulk of 
the acoustic energy propagating out of the duct, it is possible to form a set of simultaneous 
complex equations relating the measured resultant amplitudes and pha.scs to the amplitiules 
and phases of the duct modes. The modal amplitudes and phases then are calculated by 
matrix inversion, assuming that the number of microphones available is equal to the number 
of modes propagating significant acoustic energy. 


The second difference is that microphones are not located arbitrarily in the duct (often 
eciuidistant in one axial plane). Microphones are, in fact, placed at specific dud locations 
to ensure that the data | rovided by the microphones can be suitably processed to yield the 
component modal amplitudes and phases. 
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riK- mi-iluHl for ik'UMiiiininn coherent Ian soiiiul nuMlal slruetures is presented under three 
headinjrs. In (he lirsl, Method ol Moile Stmetnre Determination, the theoretical analysis 
inulerlying the method will l>e discussed including a hriel description ol the compntei 
pronram for calculating the mode structures. Under the second headini’, Selection ol 
Microphone l.ocalions, a procedure is pre.sentcil lor selecting microphone locations in 
the duct which ensures the slahility of the mode calculating procedure, hinally, under 
the heading Sensitivity t)C Method an analysis will he presented in the calculation of 
sensitivity coefTicients that can he used to ileterinine the errors in calculated modal ampli- 
tude and phase caused hy errors in acoustic measurements and microphone placement. 
.Section 3.2 is concluded hy a brief description of the User’s Manuals for the two computer 
programs discus,sed in this section. 

3.2.2 Method of Coherent Mode Structure Determination 
Ohjective 

To estahlish, program, and check out an analytical method and to establish test procedures 
for the determination of coherent fan sound mode structures. 

Approach 

Any attempt at measuring amplitudes and phases of coherent modes at tone-frequencies gen- 
erated hy a fan (or turbine) in an annular duet must addre.ss two major problem areas. First, 
during typical static fan tests, the inflow contains considerable turbulence and quasi-steady 
distortions so that any coherent acoustic modes arc generated in the presence of a multi- 
plicity of randomly generated acoustic modes. A method must thus be found to separate 
the purely coherent sound field from the total sound field. Second, the nunibcr of modes 
that a duct can support increases rapidly with frequency. For frequencies typical of most 
fan designs, there can be many hundreds of modes, both coherent and incoherent, propa- 
gating at BPb' at typical approach or takeoff operation and, of course, many more for har- 
monics of BPF. 

In considering the first problem, it will be recalled that previous mode tracking metliods 
(ref. 1 to ()) used cross-correlation or cross-spectral techniques on the output signals of pairs 
of microphones to suppre.ss the contribution of randomly generated modes. An alternate 
metliod, and the one adopted by Pratt & Whitney Aircraft under this contract, is to en- 
lumce by signal averaging the purely coherent portion of the output signal from each micro- 
phone used to obtain data. This signal enhancement at a desired frequency can be achieved 
by filtering the signal with a phase preserving fiitcr at the frequency ol interest (such as at 
BPF) and then averaging the signal over a period corresponding to that frequency. A 
oncc-per- blade reference signal obtained by a proximity sensor over the blade tips can he 
iisLul to separate the pressure signal into time blocks of one period, and a c(»nimercially 
available signal averager can sum and average the pressure signal over that period. The once- 
per-blade reference signal is also reipiired to establish a time base ln>m whicli relative phases 
at all microplume locations can be otitained. 


It) iiiMilioii lo i‘iili;ma’im'iil, llio coiilrilmlioii ol' nimlnnily lUineiiiU’il mculcs lo IIk’ 
ovL'nill si(!ii;il fan In- iviltifial liurinj! slalif U-slinji wilh a luihiilfncc siippivssit^; scrfcn ami/ 
or lu)tieyfoinl> |>lacfi) alunil Ihf iiilel. Uii|)iihlislicil results from full scale etu'ine tests con- 
iliifteil at I’ratt I't Whitney Aircraft have shown si(',i)ifieanl reiliictions in both inllitw tiir- 
tnilenee ami the siibse(|iiently {teneraleJ liistortion lone noise at HI'I'. 

As will be liisciisseil more fully in section 3,4, signal enhancing in conjunction with the inlet 
screen was foumi to be the most satisfactory inoile of ojieration liuring the experimental 
portion of the contract, providing greater repeatability and accuracy of coherent acoustic 
pressure amplitude ami phase measurements. Tor the remainder of this section, it will be 
assumed that the random modal structure is effectively .suppressed by the use of signal 
enhancing and the inlet screen. 

It remains to consider the second problem of large numbers of propagating modes. Using 
the general form for the coherent portion of the acoustic pressure in the duct as given by 
equation (3.2-1 ) I'yler and Sofrin showed that, at a given frequency, the duct can support 
only a finite number of propagating modes. If, in addition to this set, the slowly decaying 
modes that may have an impact on the fan sound field in tb.e duct arc added, the double 
sum of equation (3.2-1) can formally be replaced by a sum over that finite set of modes. 

This set of imnles. however, will, in most practical cases, be very large and a large number 
of microplione measurements would be recpiircd to provide input to a modal calculation 
procedure. 

If information were available lluit identifies the inodes likely lo contain the bulk oi the 
acoustic energy, the number of mierophones reciuired for obtaining data to provide input to 
a modal ealciilatioii ]>rncedure eouhl lie significantly retlueerl. One inelhod ol reducing the 
number of required miciopboncs is to review the fan design (blade and vane mtmbers) and 
determine from tlie kinematical I yIer-Sofrin theory (ref. 1 ) those acoustic modes likely to 
control tlie inlet and aft radiated sound fields, based on this information, locate mierophones 
in tlie duet to determine tiie structure of just those modes. As shown by Kraft and Posey 
(ref. 8). sucli a ileterminafion. however, does not always identify all of the dominant moiles. 

based on tlie above diseu.ssioii. il will he assumed that there are a finite and experimentally 
manageable number of colierenf propagating modes (or, in exceptional eases, .slowly decay- 
ing modes), iiieludiiig the possibility of reneeled waves that make up the pressure, whieli 
can be expressei! for a given frequency oj in tlie form 


P(r, e,xiw) = 


Set 

of 


(3.2-2) 


(- n/‘ 

and wliere (he plus aiul iiiimis sign aeeoiinls lespeelively lot iiu ident and rellei led w.ives. 
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ir the ampliliulc ami phase ol’ Ihe purely enherenl portion ol' llie pressure sijuial is measured 
at a mimher ol duct loeatioiis (either at llie duel wall or at any loealion in the duel) and 
inSLMivd into et|uation (3.2-2) a set o('e(|ua(ions will he generated with unknown mode 
amplitudes and phases It should he noted that the only rea.son the railial 

orders associated with a particular circumrerential mode can he determined using only wall- 
mounted pressure transducers is (hat the axial wave nmnher, in ei|uation (3,2-2) is a 
function of k,,,,^. The value of k,„^ is dilTercnt for each duct mode, I his difference gives 
rise to different spiral angles at the duet wall as each mode propagates in the duct. In some 
practical cases, for example: high radial orders of a particular circumferential order, these 
differences in spiral angles arc sufficient to allow for tlie discrimination of the radial modal 
structure using only wall-mounted microphones. 

The equation system can be written in matrix form where the matrix coefficients will be 
functions of the particular modes comprising the sound field and the microphone locations. 
If the number of microphone measurements is chosen c(|ual to the number of modes 
comprising the sound field, then in principle, the equation system can be solved for the 
mode amplitudes and phases by standard matrix inversion methods. 

Mathematics of the Metliod 

To illustrate the method of fan sound mode structure determination, the case will be con- 
sidered of a sound field consisting of t coherent modes at blade passing freciuency (BPF) in 
combination with contribution due to an incoherent source. If the incoherent portion of 
the sound field at BIM- is given by f(r. 0. x. t; cj), the pressure signal from a microphone at 
a duct location (r, 0, x) and filtered at BPF', can be written as 


xi (*>') c ^ /\ 

S m«ft«j 

-f- f ( r, » , X , ^ i w) 




wt -t 




£ 






In this equation k^j has been written as to denote that the axial wave number is a 

function of modal order, and it can be considered to be representative of incident and 
reflected waves as well as a decaying wave if k^ is complex. 


I'rom a knowledge of the fan design (bUule and vane numbers, axial flow, elc.), kj^. m. g. 
km^. and l-.,,,^ can be calculated. atui 0|„p. the modal ainpliUulc and phase, 
respectively, represent the unknowns. If (he period associated with BPF is t and the pres- 
sure signal is averageil over N revolutions. Ihe signal enhanced pressure defined in the 
time interval (0, r) becomes 


< P(r,e, . «*>) > 


JL_ 

BN 


„ akr’'* 

fliO 



-I- 


_L 

6N 


% f(r,e.x,t iw) 




1 1 


and B is the iiumbcr of roloi blades. 
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If f is corrcluteci over time sc;ilcs less Ilian some value t^. say, and if N is chosen large enough 
so tliat N > tJ^, then contributions from the second term in equation (3.2-3) can be 
neglected to a good approximation. The signal enhanced, measured output < p(r, 0, x, t;co)> 
is a sinusoidal function yielding an amplitude and a phase relative to a fixed once-per-btade 
signal that is used to “trigger” the signal averager. The filtered and signal averaged pressure 
from a microphone at a duct location (Xj^, r,,, 0^^) can then be written in the form: 


< 5w)> = j8„ e*" 


It can be seen tliat at each measurement location, two parameters (fi,^ and can be 
obtained. If measurements are made at N locations in the duct, the following set of 
etpiations can be obtained from coml.ining equations (3.2-3) and (3.2-4) for the unknown 
A|np ‘^in/L/’-'^' 




MoJtl 


iLkx Xu + >« + 

e 



(3. 


where n - 1 ,2.3, ... N. 

I'he e()uation system (3.2-5) forms the basis for the calculation of the coherent fan smnui 
modal structure. Since N was chosen to ecpial the ninnher of propagating modes, a .set of 
2N eriualions (for both the real and imaginary parts) is ohtaineil from which, in principle, 
the 2N unknowns (A,„^ and for set of N moiles) can he found. The sohition for this 
system (,'f C(|uations was programmed and is iu)w discussed. 

Computer Program 

The calculation iiroccdure to solve the system orv-qiiations (3.2-5) was programmed in 
horlran IV coile to I>e run on a tiigh speed digital computer. I he code ami logic which are 
ilescrihed in a Users Manual (rel. d). were checked out using ten analytical test cases com- 
prising a variety of fan designs and imnial structures. I he test cases used are shown in 
■faille 3.2- 1. It can lie seen that the modal calculation procedure was checked out tor cases 
(hat included up t(. 15 propagating modes. 
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I'ABU-: .1,2-1 


MODI- STRUC 'liJRi;S l‘OR TIiST CASHS 


Total 

Number 

Modes 

(u\u) Modes Used 

Input Modal Amplitudes 
(dynes/cm^l and Modal 
Phase (degrees) 

Engine/Rig 

Application 

Blade/Vane 

Number 

Frequency (Hz) 

1 

(-2.0) (-2.1) 

(-2.0): 1725 
(-2,1): 345 e‘« 

P&WA Rig 

32Z.34 

BPF 

2987 

1 

(2.0) (-.3.0) 

(2,0): 690 e'®'’ 

(-3.0): 10.0 X 10‘*^ 

I’&WARig 

32/5 

BPl- 

2 1 33 


(-4,0) (-4,1) (-4,2) 

Analytical Prediction 

P&WA Rig 

32/34 

2 BPF 

6000 


(-1.0) (-1,1) (4,0) 

Analytical Prediction 

P&WA liig 

32/5 

2 BPF 

3200 

3 

(-4,0) (-1.0) (2.0) 

Analytical Prcdiclioi 

P&WA Rig 

32/3 

BPF 

2667 

5 

(5.0) (.S.l) (5.2) (5.3) 

(20.0) 

NASA Suggestion 

NASA Rig 

15/25 

2BPF 

4470 

7 

(3.0) (2.0) (1,0) (0.0) 
(-1.0) (-2.0) (-3.0) 

All: 345 c'“° 

P&WA Rig 

-^2/1 

BPF 

2 1 33 

7 

(-4,0) (-4.1) (-4,2) (-4.3) 
(-4,4) (-4,5) (-4,6) 

Analytical Prediction 

P&WA Turbo- 
fan Simulator 

IH/40 

2BPI 

27.600 

7 

(11.0)01, l)(ll.2)(ll,3) 
(-12,0) (-12,1) (-12,2) 

Analytical Prediction 

JT8D Rcfan 

.34/23 

BPF 

2210 

IS 

(-16,0) thru (-16,14) 

Analytical Prediction 


46/108 

2 BPF 

3220 


The Modal Calculation Program (MCP) was checked out as follows: The amplitudes and 
pliases of tl)e modes in given case were assigneil eitlier arbitrary values or values from an 
analytical deck. The resultant sound field was then calculated by the MC'P from the given 
modal structure at specific locations in the duct. Hor this checkout study, the locations 
were specified by the Microphone Location Program (MLP) that will be discussed in section 
.T2.3. Ihe amplitudes and phases at these locations were then used as input to the MCP so 
that the constituent modal structure could be calculated. I'tie computer program was 
checked out by making sure that the input and calculated modal structures were identical 
within values (hat could be attributed to computer roundoff errors. In fact, all Ihe analytical 
test cases considered reproduced the input modal strueture to at least three signilicaiit figures 
in amplitude and phase. 


1 .^ 


1 
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3.2.3 Solcctioi) of Microplionc Lociitiim.s 
Objective 

To cs!alvlisli iiiul dicck out a methoii for ileterniining microplioiie locations that will pro- 
vide data input to the Modal Calculation Program .so that the numerical stability of the 
calculation procedure will be assured. 

Criteria to L’n.surc Method Stability 

It was shown in the last section (3,2.2) that the fan sound modal structure is determined 
by solving for the A,^^^’s and from the equation system given by eciuation (3.2-5). It 
can be seen readily that the solution requires a matiix inversion algorithm and, thus in order to 
get a numerically stable solution, the matrix of the equation system needs to be well con- 
ditioned in a mathematical sense. Since the matrix coefficients are functions of micro- 
phone locations, it is possible (for a .sound field consisting of a distinct set of modes) to 
select microphone locations so that the resulting matrix leads to a numerically stable 
solution. 

In order to establish a procedure for the selection of suitable microphone locations to 
determine a mode structure comprising C modes, it is necessary to consider the following 
criteria: 

1 . Microphone locations must be chosen so as to yield independent complex 
equations. 

2. I'or large number of modes, microphone locations must be chosen so that experi- 
mental and round-off errors do not accumulate significantly in tlic Modal ( alcula- 
tion procedure. 

3. Some mode combinations will require that acoustic measurements be made in a 
fan inlet with microphones on radial or axial probes in addition to those located 
on the duct wall. Viscous wakes from these probes will interact with the rotor 
to generate duct modes at DPI- and harmonics in addition to those predicted 
based on the fan design. 

4. locations clurscn can only be a function of the set of modes (m, p) considered to 
be propagating, and not on relative amplitudes atul phases (A,,,^ and of the 
modes, 

5. The method for locating microphones shoulil be sufficiently general so as to 
accommodate all possible types of nuulc structure. 


Procedure 

Addressing these criteria, a inoccduie for selection of microphone locations was devised 
and computcri/ed. I his method will now be described followed by a discussion on how 
this method addresses the above criteria. 
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(!ivcn that the soumi field in the duet is assumed to he composed of N piopapatinp duct 
modes, and given an axial region in the iluct wltcre microphones can he placeii away from 
the innuence of the direct rotor field and effects of the inlet termination, a reference micro- 
phone is located on the duct wall at an arbitrary (or convenient) a/imutha! location and 
within the specified axial region. The N modes are then ordered so that the highest radial 
order for the highest m-lohe is considered first and subsequently for that m-lobe the modes 
arc ordered in terms of decreasing radial order. This ordering is then continued for 
successively ticcreasing m-lobc numbers, 

For the first two modes in the ordered sequence, a second microphone is located on the 
duct wall and in the specified axial region using a two parameter (x, 0) stochastic search 
routine. This routine searches for a microphone location that maximizes the absolute 
magnitude of the determinant associated with the resulting 2x2 matrix equation system 
corresponding to the determination of the amplitudes and phases of the first two modes in 
the ordered sequence. The conditioning number - defined (ref. 10) as the ratio of the 
magnitudes of the largest to the smallest eigen value of the determinant of the matrix - is 
then calculated. As discussed in reference 10, this number provides a means for estimating 
the maximum number of significant figures lost in the solution of the simultaneous equations 
associated with the matrix. If the conditioning number is less than 100 (which, according 
to reference 10, translates to the loss of two or less significant figures in the resulting 
solution), or an alternate specified value if more or less accuracy is desired, the second 
microphone location is fixed and a third microphone, corresponding to the third mode in 
the ordered sequence, is located on the duct wall and in the designated axial region. Again, 
the two parameter (x, 0) search routine is used. In this case, the search routine locates 
the third microphone so a.' to maximize the absolute magnitude of the determinant 
associated with the resulting 3x3 matrix equation system. The conditioning number for 
this 3x3 matrix is then calculated and if it is less than the specified value, the third micro- 
phone location is fixed. 

Tliis process is repeated until cither all the modes are exhausted or until the conditioning 
number is greater than the specified value. If the latter event occurs, it signifies that a 
numerically stable solution is not possible for the set of modes under consideration using 
only wall mounted microphones. In this event, a microphone on a probe within the 
duct is necessary, Since radial probing (rather than axial) is generally more practical dur- 
ing both static and flight fan noise tests, the microphone is located by a two parameter 
(r, 0) stochastic search, assuming a radial probe at the furthest upstream location within 
the specified axial region. In this way, the wake from the probe interacting with the fan 
rotor and the subsequent tone noise generation will be minimized. The remaining micro- 
phones will then be placed on this probe (if possible), using a one parameter (r) stochastic 
search. At each subsequent microphone, the absolute magnitude of the detemiinant of the 
equation system is maximized, and the conditioning number is calculated. 

If a location cannot be found on the radial probe that provides a conditioning nuinber less 
than the specified value, then a two-parameter (r.O) search is made for a location in the 
duct on another probe (also (ilaceil at the furthest axial location). In the event that a 
position cannot be found that yields a satisfactory conditioning number, the last location 
of the microphone is fixed, and the program proceeds to locate the next microphone 
associated with the next mode in tlie tirdered se(|uence. 




IS 


Ihivinn iuvst'iiial (he |»rc)ceiliire, we will now tiiscuss to whiK extent the proeeduie com- 
plies with the eiileri;i est;ildisheil iitul listed etniier for loeiiting microphones. 

1 ) tty construining the conditioning inimher to he less th;in 100 tor :i siiitidily spicified 
iilteriKilive viihie) for the N x N e(|ii:itit»n system, the imignitndc of the determimint of the 
etpiation system is assured to he siiffieienlly well away from zero. Tluis, this constraint 
makes sure that the microphones are located so that each of the N complex equations 
derived is independent. 

2) The maximizing of the determinant for each additional microphone is compatible 
with generating low conditioning numherj and makes for a convenient criteria on which 
to ha.se the step-hy-step application of the stochastic search routine. Ideally, a search 
routine that maximizes the N x N determinant of the entire equation system by simultan- 
eously selecting all microphone locations would be desirable because the determinant 
always appears in the denominator of the final result. In this way, the sohition of the 
equation system woukl be independent of first-order errors in the matrix elements due to 
errors in locating microphones that will inevitably be met in practice. Such a procedure 
would be impractical however, since generally it would require a 3N parameter stochastic 
search which is not available for the range of M being considered under this contract 

(M ^ 1 5). The procedure presented relies on maintaining a fairly low conditioning number 
to ensure that small experimental and round-off errors do not accumulate significantly. 

As discussed in reference 10, a low conditioning number ensures that the number of 
significant figures lost in the calculation method is kept to a minimum. 

3) The procedure has been set up to avoid the need of radial probing where possible, and 
when it is inevitable, to minimize (he adverse effects as.soci;ited with the introduction of 
radial probing. 

4) The micropbo.'ie locations are selecteil based only on the modes considered to be 
contributing to the sound field and not on the relative amplitudes and phases of the modes. 

5) I'inally. the method is quite general, being able to locate microphones for any reason- 
ably finite set of modes. In fact, the procedure was checked out with a number of test 
cases including a 1 5 mode set (that had a conditioning number of 5,7) and a 2 mode set 
that was successfully compared with a closed form solution. 

Analytical Illustration of Method 

In order to illustrate the underlying approach for selecting microphone locations that 
assume the subseriuent numericai stability of the modal calculation method, the relatively 
simple case will be presenteil t)f determining the amplitudes of two radial orders of a given 
circumferential mode using two wall-tnounted transdueers. Incorporating the value of (he 
I'i-funclion at the wall into the A|,,j^,’s, the average pressure at a duct location of(X|,d|) 
bccotnes from otination (3.2-5). 

I A „ e t kj C ^Mt.7 
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riic pressure ill iinollier duct lociition (x-». 0 1 ), given by (x | + Xij, 0 1 + Oq) where (Xq. Oq) 
reiuesents the reliitive displiieement of the two microphones, can be written similarly as 
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(3.2-7) 


wKere 0 ■i't: 


Defining 


OfvM ■ Ami 6 
" Arti 6. 


t. [ k* y) X, + m 6i + 4 ’mi 3 

i [kx(x,^) (x,+ X«»)+ M (®>^ I 


Equations (3.2-6) and (3.2-7) become in matrix form 
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Tlius, from the measurements of the amplitude and phase of 6 / e ** *f/ and , the 

amplitudes and phases of Djj^j and can be obtained since tlie elements of the matrix 

are all known and Ap,|, Ajj^2* ^nil' ^ni2 t-'^ithtated. ('learly. the values of 

(xq. Oq) must be chosen so as not to make the exponential terms in equation (3.2-8) equal 
to one, since this would give rise to a zero determinant for the matrix and to two depend- 
ent equations. In that case, no solution for Dj,,] and D^2 ‘■'i*” obtained. Assuming that 
the determination of the 2 by 2 matrix is non-zero, the matrix equation can be inverted in 
the usual way yielding 
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Notii'i’ llint llu’ vmIiio oi ilu* (Irii'iiiiiniinl in llu' ili'iiomiiiiilor is 


t. [, Kx ( 3 ■* ( k(^3 J Xo 


:iiul llio only ii';ison lluit IIk' rii(ii;il iiunli's Ih* lU'ts'nuiiK'd iisiii}’ only duel w;ill iiiicro- 
plioiu'K is tluil tiK' ;iNi:d wave nun her is a lunelion of the eireumrerential and radial 
inodes llial are propajialinp. No'e also ihal even for this relatively simple ease, poor 
resolution due to the ''enoiniiialor o( ei|ualion (,V2-d) approaching /.ero is unavoiilable if 
the axial wave numbers (k^.) lor the two inoiles are almost eipial aiui Xg is boiiiuled by 
piactical constaiiits. Such eases olten are met when several raiiial orilers are propagatinji 
for a particular ni-lobe lircuinlerential moile at high fre(|iiencies, In particular, if m is 
fairly high, the axial wave lUiinbers of the first few radial orders for typical practical 
situations are indistinguishable. In such cases, the detailed mode structure can best be 
obtained by probing within the duct. Poor resolution is also experienced if the distance 
between microphones is chosen such that 

e ^ ^ 

iToni inspections of the form of the determinant, etpiation (3.2-10), it can be seen that a 
maximum value ol the moihilus and the real part ol the ileterniinant is obtained when 
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(3.2-11) 


Two criteria are satisfied by this choice of Xg: 1 ) both the modulus and the real part of 
the determinant are well away (rom zero, and 2) both the modulus and the real part of the 
determinant are at local maxima and will be insensitive to llrst order errors in Xq (that is, 
the microphone locations). Having established the microphone location accoriling to the 
form of the denominator of equation (3.2-')), the actual values of l)„g and D,,,-, depend 
r)ii the form ol the numerator and the measured amplitudes aiul jihases of the pressure at 
the two microphone locations. However, as long as (he matrix-vector product in the 
numerator is bounded, singularities and gross sensitivity to llrst-order errors in measure- 
ments and microphone locations can be avoided. 


l or this simple case, (he stochastic search that maximi/es (lie ileterminant for each addi- 
tional microphone location in the Ml l» would determine one of the values of Xg satisfying 
etpiation (3.2-1 I ), aiul within the specilieil axial range for locating microphones. Iftme ol 
the values of xg is iiol in the specified axial range, ihe stochastic search will find another 
value ot Xg Ihal will inaxiini/e the delerininaiit in that range. Since the eigen values of the 
rieterminani in the ileiiomiiiator of cipiatioii (3..’ ')) are (1 + i) aiul ( I - i) when a value ol 
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is v'liosi'ii l)i;il sMiislii's i*(|uati(>ii I 1 ), ilii- i iitiililionini' luimiH'i is UMilily m\'ii In )n< 
imily. I'lu' sysu-in ivstilis tioni iln' dioiiv u| inii ioiilioiu' liu aiinns liy \[) 

in ^'(jtialion (.^.1 1 1 ) is lluis si-oii lo Ih' wvII i omiiliniu'il. 

('ompiiU'i' IVoi'i'iim 

riu‘ mill nplioiH' liualioii svK'ilion iiiiiii-iluiv was ptut'iainincil in I oKian IV auk' lo lu* 
run on a hipli-spi'i'il ilij'.ital loinpuiov. Tlu' I'oik* atui loi'io, wlncli aiv ili'si iilu’il in a Usi't's 
Manual (ivT. 1 1 ) wi’iv fhaki'il oul ky loi-aliiii', Iwo luinoplioni's loi a aisi- wlu'iv tlu' 
souiul lidil 1 ’onsisU‘il ol' tsvo propapatin}’. inoiU's. riiv annpulm oulpnl was i inupaivil with 
llu' exact analyliial solution olnainai hy ai>|»lii‘alion oldiualion (.1.2-1 1 ). Apivi'iiu'iil lo 
williin lour sipnillcanl lipuivs was oklainal. As nolat in sci lion .1.2.2, llie Mierophoiu' 

I ocalion I’ropram (MI.IM also was run lo obtain niicroplioiu' ku'alions lor Itie U'n 
analytical (csl cases that were uscil to elieek oui llie MCI’. I'inally, once Hie Mi l’ was 
elieekeil out. il was useil lo .seleel mieroplione loealions lor (lie exporiinenlal portion oC 
this program ileseriheil in section .1.4. 

In this seeiion a inelhoil has been eslablisluul aiul cheeketl oul for tlelerminiii!’. microphone 
locations lhat will provide input lo the Modal Calculation I’ropram (de.scribed in section 
3.2.21 tlial will ensure the numerical stability of the calculation procainre. I he method 
for delcrniininp coherent fan sound mode siructures is thus complete. In ilic next 
section. Ihe sensitivity will be considered ol (be modal calculation method lo errors in tlic 
knowledge of the nlea^alrinf^ loealions ami in the measured phase anples and pressure 
amplitudes. 

3.2.4 Sensitivity of Method 

Objective 

l o deiive a set of influence coefTicients from which the sensitivity can be determined of 
the calculated mode structure lo errors in the knowledge of tlie measurini; location and in 
the measured phase antxles and pressure amplitudes. 

Derivation r»f Inniieiiee Coeffieieiit.s 

The accuracy v>f the modal coelTicients, as computed from the measured pressures, is 
affected by errors in these measurements, and also by Ihe precision of the microphone 
loealions. Spi-cillcially Ihe amplitude and phase errors associated with a specifieii mode 
are rapiired as functions of amplitude and phase errors in the measured pressures, and 
also as functions of errors (or loleraiicesl in Ihe radial, axial, and circumferential 
coordinates of Ihe pressure measurinj’ locations. In this subsection, exi'iessions will be 
obtained for Ihe inlluence coefficients or partial derivatives of the modal amplitude and 
and I'll ISC with respect lo pressure and location comdinalcs. I hese can be used to 
estimate (he sensiliciiy of the calculated moiles (o firsl-order measurement and position 
errors, by multiplyint; a I’iven inlluence coefficient by an appropriate measurement 
error. Ihe correspondini', error m Ihe amplitude or phase of the mode due to that 
error is nblained. In addition. Ihe standard deviations of calculated mode aini'liluile 
and I'hasc can be obtained in leinis of assumed or exi'ciimentally deiived 'laiulatd 
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ik'viiilidiis of tiu' mi’iisinvmctit uml pusilion errors. 'I'lie deriviitioii ofllic iiiniunce 
eoerikiiMits iiiul their use in estimating the stanilard deviations in calculated modal 
amplitude and phase are discussed in the following paragraphs, 'I’his discussion is 
then followed by a physical interpretation of the intluence functions. 

The analysis is simplitied and will he easier to follow if all appropriate complex variables 
are denoted initially by single symbols rather than exprcssetl in rectangular or polar form, 
until such time as is necessary to specirically obtain separate amplitude and phase expres- 
sions. b'urthcrmorc, matrix notation is employed as the most elTicicnt and revealing way 
of manipulating the large number of variables involved. 

In matrix form, the relation hetv/ecn N measured pressures and mode coefficients (equation 
3.2-5) can be expressed compactly as 


P=Hz (3.2-12) 

where is the column matrix of N measured complex pressures, z is the column ol N 
desired complex mode coefficients, and H is an N x N matrix, the elements of which arc 
functions of microphone location and modal eigenvalues and eigenfunctions (but not 
mode amiditude or phase). Tlie expressions for the H matrix elements (denoted by M,^^) 
can be written from equation (3.2-2) as 
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ns 
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(3.2-13) 


where (x.,, r,,, 0,^) denotes the nominal location of the nth microphone and s denotes 
the sth mode in the ordered seqi ence of modes discussed in the iast sub-section and m^, 
and kj^j. are respectively the circumferential lobe, radial order and axial wave-number 
corresponding to the sth mode. Solving equation (3.2-12) for the modal coefficients 


where 


or 

Or Q 'S -bke mVftfic. of i U H 


(3.2-14) 


It should be understood that although the elements of II have closed form analytical 
expressions given by (3.2-13) the elements of the inverse, 11'^ cannot generally be expres.sed 
analytically because of the complicated form of the determinants involved in matrix 
inversion, However, in Ihe modal calculation procedure described in sections 3.2.2. and 3. 
numerical values of the II matrix are obtained, and the computation procedures included 
in the MCI' give corresponding numerical values of Ihe inverse, Q “ II * . I'urthermore. p 
is known from N duct microphone measurements and Ihe resulting compuled quanlilies 
z (the amplitude anil phase of each mode) are known from the output of the modal 
calculation procedure. 
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riu* ivsul) on Ihc calcii!:ili‘ii modal lOi’ITidcnts of errors in the measurement ofp and 
errors in mieroplione loealion, whieli are manilest as errors in II, tmiy he seen by replaeinp 
p, 11, and /. in ecpiation (.1.2-1 2) with p+ Ap, II + All, and / h Az. .Snbslilnlint>. and sub- 
lractin{; p - Ilz from eaeli side gives 

Ap= H4Z + AH-E + 

ir the errors are small, the second-order term All Az ean be ignored, and to denote this 
approximation, dp, dH, and dz will he employed in the usual sense of differentials, giving 


cif? = (.1.2-K 

Using standard matrix algebra the error in the mode coefficient can be expressed in terms 
of pressure measurement and microphone placement errors as 

= H'' Jp ~ 

or 'Al Q - Q aH-i 

Since Q can be numerically calculated from the known matrix II, and the z are 
calculated using the modal calculation procedure, it can be seen that the errors in the 
amplitudes and phases of all the modes dz can be expressed in terms of the known 
elements of Q, the calculated mode structure z, and the errors in pressure measurement, 
dp. and the microphone placement errors, that give ri.se to errors, dll in the elements 
of the matrix 11. 


In particular, from equations (3.2-1 7) and the form of tlic expression for the elements ot 
H given by equation (3.2-13). 

where, since II is a function only of microphone location, dll has been replaced by 





l-rom this equation, tlie itidependeiit effect of each error type in pressure and loealion 
coordinate can be evaluated, b'or example, tlie error in the modal coefficient /|^;|. due to 
an error only in x, say. is given by: 

“ Q XT. ^ 

Multiplying out the last Iwolerms as is shown in equation Iil3 of Appetulix H. 

.U 
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wlKMv U , . is tlu- clum'nt in Ihc n“' row :«m! s‘'' volumn of llu' in;ilrix H, :m.l VS Hh' 
st!i modii^ i-odTii'ii’nt ns orilnvil in tlu* modal t-nkulalion mothod. 

TIh’ (.'xpi L-ssion 
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is ddinod as an inriucncc aK’lTicicnt and is simplified still furtlicr. alonp with :dl tlie 
other inniience etieillcients. in Appendix B. It'the measured aniplitiidc and phase ol the 
coherent sound Held at the t/'* microphone at location (x„. r,,. 0„) «s given by H„ and i^„, 
titen the iniluence coelTicients can be i>resented in terms of amplitude and phase and aie 
eouivalent to partial deviations of tlie modal amplitude and phase with respect to 
acoustic pressure and microphone location coordinates. As shown m Appendix 13. these 
partial derivatives are as follows 
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The remainder of the notation is fully described in Appendix A. 

It sliould be noticed from these partial derivatives that, in general, tlie iniluence 
coefficients are functions of llv: microphone location (x„. r„. 0„) and the acoustic pressure 
mcosurciiKiils (»„ alul *„), ami also ac|.oml on Iho calculalml iiioJal “'"I’ 

(Am a,ul Thus II.C Inllm-m-o om-fridonls arc- .llllcscnl In, aaoh .nodal slnKluK II...I 

is determined. 

I xpressions for the standard deviations of A^ and 0^ can now be obtained in terms 
of the standard rieviations in pressure measurement and microphone location errors. 

Since the resulting amplitude and phase errors of the Mth modal coellicients dA]^^ and 
are linear fimclions of the input errors, they can be written as 


j A - 1 iii: ^ "k ^ * a 

it A = ^ ( H: ■*“* * " if; ^ if: ^ 


I lu‘ii, if il is asMtnu'il ihul llu- pivssiin.' and position riioiN t(l|j, il^^. d^, ilj., and d^y I aiv nnU- 
poiulnil variahli's havini' ( laussian di-.lril)nl ion linu lion > dial aio indi']H'iidrn1 of n and if 
these staiuiaiil deviations are eharae(eri/ed liy zero mean ami siaiutard deviations (i|j. 

and 0 (f. the resniiintt standard deviations of ilie eoinpnled mode amplilmles and pliases 
are simply 



Both influence coeffieienls (eiiuations .1.2*lX)aiul standard deviations ami 
are caleidated as standard procedure wlien the MCP calculates the amplitiules aiul (diases 
of the acoustic duct modes. I luis, a signillcant feature of the modal calculation methoil is 
that the sensitivity of the method to errors in measurement and microphone location can 
be assessed for every moiial structure that is calculated, fhe basis for this assessment is 
made by inputting to the MCI* estimates of the standard deviations of the errors in 
pressure measurement ami microphone location, lisliinates of the errors in the amplilmles 
and pliases of each of the calculated modes arc obtained as an output of the M( I*, l or the 
remainder of this section, the form of the influence coefficients will be discussed and a 
physical interpretation will be pre.senled for those cases where it is reasonably clear. 

Phy.sical Interpretation of Influence ('ocfficienls 

It is not seen readily from the expressions given in equation (,1.2-IX) for the influence 
coefllcients, how errors in calculated modal amplitude ami phase are inHueuced 
by acoustic measurement and microphone placement enans. Il is apparent, however, that 
a significant term in all the inlluence coelTicienls is q The (j in turn are com(ili- 
cated functions Ilf the particular mode set that is being determined, along wilh the 
microphone locations used to obtain data foi' input to the MCI*. 

In detail, the inlluence coefficients e.xhibil ilifferenl characlerislics depending on whether 
they arc associateil with microphone location errors ot acoustic pressure measurement 
errors. 'Ilie inlluence coefficient charaeterisiics will thus be eonsidereil separately under 
these heailings. 

Miero|)lu)ne t.ocation lirrm-s 

If the mode amplitude and phase sensitivities to angnlai or detrots iidj, and ‘•0,,,. dd,, 

aie cottsidered first, it can be seeit from the ei|u,ilion ( t,.' IS) that lhe\ tlepend diieeth 
ott the values m^ which are the set (d ititcgen. giviiij’ the iiumhei id eireuinferenlial e\clcs 
of the sih spinning mode patlern. If all oi most of the in,^ are laij’c. the seiisitivitv is high 
coinparetl wilh cases where all the m,^ are relatively small nuitibcis. A liij’h value oi m in 


Jo ^Sh;:";::;::^ 

location, dx„. ^m c"o« in axial 

1 lie physical interpretation of sensitivities to errors in radial position 3 A /3r and 3rfi / 

k T- '(k 1 tT' ““'"‘"‘"S ‘lie "xpresion to k„ J, *'"= 

«.n Thtrctorc, pressure measurements at or near the duct walls are verv in ’ 

(ni. 1 ) and (m,2) moTelt in "r* " ' (m.O), 

(for liigh hub-tip ratio) would be near tlip'^^ r” *? ^ ^ neighborhood of the mean radius 

no.lo Of ,l.c (m.r) .uol yL.,lL"V„ige X ^ 'oc«„„ would he vc^ near ,o ll,u 

depend on the relative amplitudes of the (m. I ) and"('m,2) iod™ * *"^»Wty would 

The eigenvalue, k,, also enters the sensitivity as a factor -ind i • 

relalively large pre.ss„r^C” X:;e!l' t^ S r;.d;a?r'"“' 

more gradual, as is the ease to low order radial modes (i.o low value* “ 

Prc'sure Measurement Errors 

errors ill nieusured'ir'res^^^^^^^ ampllnide aVl''' 1'"''““^^^^ eocITicienIs giving llie effects of 
suiting modal 00^;! hie amplitude and phase of the re- 

solely of the fact that the ,|uaulilies involved Ire auupkx vaSto ^llv,^d”’^ 

r;id .-sun pS:::; 
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Wlioiv A iiul li aio u\ tlu* m"’ imnlc ami moasuivil pivssiiivs al liie n‘ 

";M.ccliv;'ly, ;.n,l llu- iIk- mwnilu.lcs ,.r tlic m»mx 0,„„ oK'.m-nls, ,o , 

ihc error in a typical resultant can he ohiamal Irom the ililtoienti.il as. 
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The inHuenee cootTieients. are seen to he simply M,„„. or just the elements of the 

Q matrix. 

In the ■letinl case etiiuitions (3.2-14) are seen to contain q,„,^ as in the above hypothetical, 
real variable ease but they also inelmie other factors, incliulinp harmonic Uinetions ol an 
mi tat 1 coml natiot of the phases of the measured and calculated pressures and the 
of 0 ' Tl. Vlu,uld be uoHcea. however, .ba, the l„n,.enec e«IT,eiems u. this seehr.h 
are not functions of the calculated modal amplitudes). 

To get a general appreciation for the nature of this eoinplieution, eonsnier the ease of two 
eninplex nioJes amt two complex pressures, for illnstration equation (3-- Ml becomes 


(3.2-21) 


(3.2-22) 


2;, i (?« f. + ^1*- P*- 

2^ ^ Qn p. + P*- 

Tlie (complex) error reiations again can be found from the dillerendal; 

c[2. = 0.1 J Pi + 
c| = <5t( Pi Qii ^ ?>■ 

Now these expressions are just as simple structurally as the liypothetical real variable ease 
given by equation (3.2-20). It is only when we obtain the separate eltects on moda am i- 
tude aiul phase of both measured pressure amplitude and measurctl pressuie phase that c 
eira factom (apart from q. the amplitude of Q) in the inlluence coetlicent expressions 

appear. 

To suggest the nature of this complication, let us examine /., and d/, . in the above equations 
(3.2-21) and (3.2-22). I'or simplieily 

2, - 0,1 Pi + Oii ft 

£ , = tJi + Ut 

Ihen the complex error in /., due to errors in the complex eointmnents. U , - U, ,1’, . am! 

(.It - (J| ilS can he ex|nessed as the dillerential: 




cii, flu, + clUt 


I■inlllv J-l shows/.] ;is iIk’ ivsultiinl of U I ;itul Iliich t.-ompli-.x I'lior liU j aiul iIUt is 
shown R'solval into two components: 

1. A component colincni witli the appropriate U vector. I'his component is the 
ilifferential of the amplitude of II or the error in tlie U amplitude. 

2. A component perpendicular to U which is etjual to the error in the angle of U 
multiplied hy the magnitude of IJ. 

lugiire 3.2-1 shows the projections of these two different components of the errors dll | and 
dU -1 upon the resultant factor / 1 . fhe sum of these errors may be seen to give the error 
dAj in the amplitude of the resultant / |. The trigonometry involved in describing this 
projection process mathematically can be shown to account for the factors for 9A^/3B,^ 
and 3An,/9i/^„ in equations (3.2-18). 

A similar process (not illustrated) may. be used to tlnd the error in the phase of the resultant 
/] . This time the same four U-component errors are projected on an axis perpendicular to 
Z| . The sum of these perpendicular components is equivalent to the error in the phase of 
Zj multiplied by A j . the magnitude of Zj . On dividing by A| , the phase error of the 
resultant is obtained. This procedure may be shown to give the mathematical expressions 
for the pair of inlluence coefficients and 30 ill/s'/'ll equations (3.2-18). 
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The main conclusions from the above description of the physical characteristics of the 
influence coefficieiits for both microphone location errors and acoustic pressure measure- 
ntenl errors are as follows: 

I'or microphone location errors: 

i) Influence coetTicicnts are functions of the calculated modal amplitude and 
phase as well as of microphone location. 

ii) Calculated modes with large circumferential, axial or radial wave numbers are 
sensitive to errors in locating microphones. 

For pressure measurement errors: 

i) The influence coefficients are functions of measurement position and the cal- 
culated modal phase angle only. 

3.2.5 User’s Manual for Computer Programs 

A separate Usc’-’s Manual has been written for the MCP and the MLP that were developed 
and checked out in section 3.2.2 and 3,2.3, respectively. They are both pre.scntcd in a self 
contained report in references 9 and 1 1. 


3.3 APPROACH FOR MEASURING RANDOM SOUND FIELDS 
3.3.1 Introduction 

In addition to coherent discrete frec|ucncy noise resulting from rotor-stator interaction, fan 
noise contains significant random or broadband content produced by rotor-turbulence 
interaction, for example. The objectives of the analytical investigation conducted in this 
part of the program were to derive a means for specifying the essential properties of a 
random noise field in a cylindrical duct, and then to evolve instrumentation and data re- 
duction technicpics that permit experimental determination of these properties. 

To simplify the problem, which is more complicated than the analysis of coherent, discrete 
frc(iucncy fields, attention is focused on thin annular ducts, where radial variations in 
the .sound pressure are negligible. The content is divided into these parts: 

Section 3.3.2, "Field Structure of Random Noise in Thin Annular Ducts”, formulates the 
concept of the imula! pot er spectral density function, its expression in terms of cross 
property functions of local pres.sures, and interpreialions of the function including spinning 
tnodes. The elements of data ac(|uisition and processing systems for experimental deter- 
mination of the field structure are discussed. 




SL'i'liim .'..^.3 lU'si’i ilH’s :i special syslcm I'oi incasuriiip i.iiulom I’icMs in tliin aniuilar ilncls 
an ana\ nl' ciicnmrcivnlially spnceil liclcciors, iiuoiinnalinp progressive linear Jelay limes. 
The rnnelion ol lhe array is lo provide a ilirecl measure of ihe set of nuHle-lieiineney pairs 
havinp eommon eiilolT ralios. Array pei lbrmatiee eliaraeleristies are ileriveil tor hollt ran- 
dom and eolierenl. diserele-rieiineney sound Helds. 

3.3.2 l•'ieltl Slriicturi of Kaiulom Noise in Thin Anmilar Ducts 

liuckgromul 

In order to appreciate the ohjeelives and approach of this analylieal stiuiy, it is lielpful 
lo mention certain feaUiivs of the current theory for dueled coherent sound Helds. 

The theory of coherent, discrete wave propagation in cylindrical duels lias been eslahiished 
for some lime and has been validated by a number of experimental investigations. A 
centra! item in both formulation of this theory and also in determining the lype of measure- 
ments reipiired to acquire meaningful experimental lesitlts is the appropriate specification 
of the pressure distribution in a iransverse plane perpendicular to tlie duel axis. The 
radial structure of this distribution is expressible in terms of, sums of Bessel functions of 
order corresponding to the cycles of circumferential pressure variation. Thus, the trans- 
verse pressure siniclure depends crucially on a very clear understanding of tlie character- 
istics of the circumferential pressure variation at a fixed radial location, typically along the 
circumference of the iluct wall, flic features of this circumferential behavior are embodied 
ill the case of a very thin annular duct having, at any instant, a constant pressure across its 
very narrow amudus. 

By focusing on this simpler case it is possible to work with relations that are unencumbered 
by Bessel Hiticlions and other coniidications that lend lo obscure certain basic features. 

In this way it has been pos.sible lo recogni/c clearly that so-called spinning nuKies area 
useful aiul convenieni way of ilcscribing the coherent pressure Held of a rotor aiul the Helds 
generated by rotor-stator interactions. In particular, simple interpretations of propagation 
direction, cutoff, and other properties follow from this thin annulus approach aiul are used 
to guide and interpret the results of studies of Ihe more general cases where Ihe duct 
annulus is of finite size. In the current study of random noise, the approach of formulating 
the |>robleiii for the case of a thin annular duct has been retained. 

An analysis of ratuiom Helds in cylindrical ducts (ref. 4 aiul M has been done previously 
by postulating Ihe existence of a spinning mode ampliliule spectral density Hmclion. 
However, to exieiul the eolierenl spinning mode formulations to describe random pressure 
Helds requires that the concepts of “amplitude" and "phase" of such a mode be modilled 
lo accoiiimoilale random time variations. Providing, an ;i priori physical basis for siieli a 
modificalioii is difficult, and it was therefore decided in the present study to adopt a 
formally different model lor the space-time ilescrii'liim of a random pressure Held one 
that did not emluuly s|>inning modes as its starting point, fins model s]u'cifies tlie Held 
by pioihiels ol i*airs ol liiiu lions fiiiu lions ol space only that arc strictly periodic, ami 
fuiu liotis ot lime alone that salisly onl\ Uie usual requncmeiils lot random lime series. 
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tunction. 

n,vUly the modal psd function is resolved into two component 

ilSn'o?thL^ep^^^^^^^ 3.3.3) that these random field spinning modes '‘exist m the 

operationally definable and measurable sense. 

Mathematical Model of Random Pressure Field 

A coluTCiil discrete pA-ssurc nd J in » thin anmdar duct may he specined by ll.e super- 
position of spinning modes in tlu. usual way by: 


pCe t) ~ ^ X ^os ' n uxi} + 

r > /)-< 
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(3.3-1) 


ruis form cannot conveniently be modified to describe a random field iKcanse f',,,,, and 
0,,^,^ must incorpor 
provided. 


JllVtllXIlUY i-'v inv/vi*t.-v .4 1,, 

ate time dependence, and a continuous fraiucncy variation must be 


Instead, the following 
in a thin annular duct: 


fonn will be taken as the mathematical model of a random field 


pCPjt) = ^ Co5(ma) + bn,Ci^J 


(3.3-2) 


where a,„{t) and b ,,(t) 


are random lime functions characterized only by their cross ami 

, Ml V • 

direct correlation and power spectral functions. 

Kquation 3.3-2 is seen to satisfy the following rcMuiremcnls for a model of the pressure 
field: 

1 Time dependence at nil positions, 0, that is arbitratily random in liic eonvenli.nial 
sense of the term as used in the context of time senes (rel. U). 

2. Spatial variation at any instant of time that is strictly periodic in 0 . with period 


Zir. 


2d 


i 




I’owi'r Spi'cti'iil Density of I’rcssiiit' I’idtl 


rile ohjeelive is to lierive ii sniliitile rmielion. slarliiip with i> (fl,l ( j'iveii hy etiiialion 1 ), 
that will speeiCy itie Jisliihiilion ofaeoiisiie pressure as a I'unelioii of mode miinher, 

(M . ami rreiiueiicy. to. Miaiis for olMaining this I'unelioii in terms of e.sperinieiilal data 
are also ret|uired. I liis rimetion is called the modal power spectral density, ,S(m, to). 

Two hasic steps are used to ohtaiii an expression for S(ni , to), l•'irsl. the cross-correlation 
function corresponding to etiuation (.V.^-2) for p (t).l) is rornuilaled. Secondly, l\y virtue 
of the relation that power spectral ilensity and correlalion functions are I'ouricr transform 
l>airs, an immediate statement of the funclion S(-w. to) is tibtainod. 


file correlalion function of a function tliat is periodic with re.spect to 0 and random with 
respect to lime can he written as: 






27, 

fKe,t)‘ p(et<r^ ■b-ft') 





It will he useful to define the partial correlation function. R^tOiff-.T) hy means of; 

7-->oo r ^ P(e,t) - P(e + <r jtf 'if) Jt 


~ p‘pc<r,r) 


so that equation (.V.V.l) can he c.\pres.sed: 


^2.77' 

Rc(r,r) - j R^(3y(rr)do 

o 

I he integral K^^U);0'.T) is simply the lime-correlation of the jiair of pressure signals at 
0 and t) + (T; 

the resulting modal powei spectral density, St»i. to) is the telouhle) I tHirier transform ol 
kl<r . -JT), which hecause ol IheO^- periodicity (ris the tiitference in 0 coortlinatcs) Iiecomes; 






/; 

i? '«• 






On climinatinp R((T,a) botwecn c«|uations 
modal power spectral is obtained: 


s (3. .1-4) and<.1..1-6). a {-.eiKMal expression for the 


Scw.io-) J j J J dt^0M<r 

This BxpKSsioh for S(« . «) will be cyaluated ih Icnhsof Ihc sirccific model for (O.D 
given by equation (i.3-2) by executing the following sequence of steps. 

The explicit form for p (d.t) given by equation (,t.3-2) is s^rstithted in equation 
(3.3-4) to obtain the partial correlation function Rg(9;cr,'l). 


1 . 

2 . 


^ ^ 4- is applied to 

The Founer time transform 

Rfl(0 ; CT','2) to obtain a local cross spectral function expression m terms of moda. 
spectral functions. 

3. A space-averaging operation. ^ -* is applied to the result of step 2. The 

result will be a spectral function in the form of a Fourier series with argument O'. 

4 By applying the Fourier transform for <T-^ ^ in equations (3.3-6) or (3.3-7). 
■ ihe coefficients of this series are obtained. These coefficients are thejesired modal 
power spectral density functions, and the set of these functions is S(i7t . oi). 

Thus, the four integrations involved in obtaining S(wi , oj) by equation (3.3-7) art executed 
in the order: ©> (T2 

step I Substituting equation (3.3-2) for p (O.t) in equation (3.3-4) for the partial 

correlation coefneient, Rp, and using indices p and q for m in the double 
summation gives: 


fte Ce; <r, 2-) f f^t :[ | + b^Ct)^ pe]} • 

m 

co4jn bnJ3^) f 


pn 


n 
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Nvl.c-iv )!k- I p,, IKIV.' IIU' trnon.1 sinulmv 






S10|. J Al'I'lyiW 

V 

{R& C^; cr;r)} = ^ ^ 

/M ^ ^ 

ni ^ 

|.l„s sum or l ourior mmsromis of tl.c I',,,, k-rms in o.|uallo„ («1. 

. ,- .• I I n .« ilu< iiot-ition for tho 1-ouricr tnmsform imirs in L-qualions (S) 

TIu' followinp listing ilelincs tin. notation loi 

ami 


CoRRELflriOf^ 

PUNCTlO^i 


an^a^ct) 

UTbZc^) 

dm 

b^*«T2'» 


SPECTRfii- 

fuNcno^J 

■Pi 


In cMiuition (,r..V‘i) only llm lorms in Ilk- 

::;t;u.rt "Xui'i^ ,X',«uUinr fro.'.: .i,; 1™,.*^,™. «. ,i,y ^ 1-.,^ U, y.,ui,Uo„ 

(.^.3-S) involvi.' fl-iU'potulLMK'c. 

S|C |1 ,r A|i|'ly ll'O O-iivyrapini! o|H'r;ilion iviiniiva l>y olimtion (.l..r-«) lo Mlll.ilion 1. .. 
•Iiu. ,un hiuui si.ir oi r.i-i'i- 

::;;:xr:;::= 

connoctions. IWcaiisi- of its importamr H .s ivinvsciilml b> Hu ssmi U t. , J 


Wilv’iv 






'.n 

) J cl' i Qn (0/ . '^ )} 


App!ii;ilini. o( (hi- .'iH'nilion X lo llio liphl siilo nro(|M:ili«Hi ( <..<-'>} 

m.ikos iilt li’iius luvinj’ luimumii.’ W-ili‘piMuk'iia' vunislj. I^.•:lvill}', only llu’ followiiii' siinplf 
Slim i'lom Iho lirsl lim* (»r I'qmiliim 




t rill’ symbols '1^ cimjiqMlo (oj) have replaeeii 

Since I he I •tmrier Iransrorms 

sliown to be conju^uile. 



atui Vr« ‘'1'” reatlily be 


hitroiliicitiu (be ileti'iitioii: 

i 6 ~ C C^) ] 

l-,i]uatioii (3-M n becomes: 


Crerl ) 


(3.3-12) 


^ r 

= 2 1 i + ^ r/ (3..M3) 

/w-o 

Step 4 l oiirier transrormint; equation (3.3-13) with respect to 0" completes the 
process for obtaining the modal power spectral density. Sfiri , oj). 


r' -*/w0^ 

S tv) = ^ 

e» -» q f 

~ j ^ f e J(r 

nrt> ~yr 
nro^ 


(3.3-14) 


The dummy index t) has replaced m as the index of summation in equation (3. .’-1 3) in 
keeping with conventional practice. The f - functions required in the f-outier transforma- 
tion of harmonic functions are defined to unity operalois when their argument is zero and 
otherwise arc annhilating operators. Collecting tcmis with and \(jvi+h) laetois. 

and replacing the pure imaginary S(iM . to) can be expressed in terms 

of real functions: 


I'OVii 




t.; 


(3..M5) 


Ivllinf- 




(3.3- 1 6a) 


/ (3..3-16b) 

U ('<*>) = irj M " 


Uk' moiial power spoctral ik-nsily bocomos: 


5im,i.iy) UC<^) SC'^-'^') + -4 Vn(^) S (W^n)J (3.3-17) 

f}=a 


I bis eqiKilioii indicales (liat S(vn. w) is tk-rined over both positive and negative values of 
^ aiul oj; I’be U and V functions, being combinations of functions icsulling from I'ourici 
transforms arc perforce delineil over the entire u> ilomain. Also the (unction ^ (yn ■hP') 
directs to fill the negalive^-domaiii of S(-jrt . co)jusi as tills its positive 

domain with {Jn (i^)- 

Now by combining these characteristics with the symmetry properties ol the constituent 
and functions, S(yji , to) can be simplilietl. t or this purpose it shoukt be 

noted that since results from I'ourier translorming the auto correlation functions of 

the a«Ce) andt*(-t/ series it is an even function of to.^ Similarly the combination ol 
cross spectral functions of these series that comprises (sec eci. 3.3-1 2) resullsjji 

being an odd function of to. ( This can be shown by recogni/ing that Cftwb« 

) is an odd function of2C I'lien converting the factor exp-itor to rec- 
tangular form when writing the I'ourier transibrm to obtain it will become evident 

that is pure iniaginarx' and ati oild function of to.) 

I'he lollowing itnpoi tatil relations are consei|uence oi the above ptoperties: 


or 
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Un ' Vn 0^^) 

\Jn (-U ') =■ On 0t>) 


(3.3-lSi 


aiul 


or 


5(r>H,~Ur) = 


( 3 . 3 - 19 ) 


way, as shown i,,Fig„a^3.'lufxl^^^^ 

Lth qLadra u n anriv “l''™' »P«incatlon of m 

sen, a, Ion on, yin qnalL, I V til: of’™ 

values of >71 in Vi component functions for negative 

possibility orcaluajing ,ho fnneions for na^tl.e J u”lf » V° 

p;ior^£;:?sS,^^tiJ iiret^^ r? 

^ H 5^m H • The + and - superscripts are used to indicate opposite and similar 
pairings, respectively, of the signs of m and oj in the corres]ionding modal psd, S(m.co). 

-"Vl "r'' -vcy 

p-'ts of forward -,n interpreted as the random field counter- 

colrerent case. PC^,4 ‘wSeleT' 

- ToXw? ^ *P'" t*Ps mode is therefore Jl = /10/Ui 

. . Similarly, a reverse spinning mode, in which Tn and to taki» 

ha same sign, has argon, an, ^n-lindica.ing that pressnra a, AS distant from e L 
received earlier or l^t s ^ ^ uisiant troni 0 was 


.15 





t igtirt’ MotlaU'owcr Spci iral l)vn\iiy. S (m. ui I Point Wilui-s 
(HI (IS a subscript is itlways jiositivc) 


m. to HAVE SAME 
SIGNS IN S (m. w) 


m, to HAVE OP 
POSITE SIGNS IN 
S (m, to) 



l iyinr yt ih 


Point I '(lilies oi ,v^^' / 


(jJ I I iind .S 


in 


till in 


i to 1/ f oinpoih nis ol V tin. to 


(/ siib\, iipi IS nhvitvs positive ! 


ORIGINATE PAGK 1:' 
OP P(M)H \ 


,l,n.K- 1 P ovionsly ill IlK' Jovdopmcm ..I S w,«)l .... IIk' ruiulon. IlK ' ' 

K. I, su„lc pair ..r S- Ipnclio.- c»d. 

- IT h is or cLun^c the welUnown expression lor the “power" or mean-square 

pressure of this iorwanl moiie. 

IXspita llK similaritk-s .Icscrihcl above in Hu-so ran.lon. ami 

be t'ouiul a physical or operational basis for use ol tliese terms. In i.ondudmg iis ^ 

,.r or kTiv I rtlio modal power spectral .lensity. RatI.er, it emerBcl as a component of 

nt: moL. ps.1 lunction tbat satistled attain riTs ,own 

With coherent, discrete frequency tleld spinning modes. But lurthcrmorc. 
in section 3 3 3 random field spinning modes turn out also to liave a very dose pliys u 
4e,Sce (cienm-cl in ter,na of measurable atlributcs) to ti.e.r col.crent pressure- neld 

counterparts. 

Tbc followinf e„ualionssummari/.e tire expre-ssions for tbc spi.minp n,ode^comp..nenls. 
s tfl.ol) ami . of llie moilal power spectral ilenslly lunctron. 6( . 


I'orward spin component. 

luive opposite 
sign in i 


SmO^^O ~ SCtT),'*^) 5C-W, w.?) 

:ir 2 . 

= (|U>/) (^tUx 


S(r(\ , lO) 

Reverse spin component. 
1 M ^ u? liave same 
sign in SOi^j * 


5^0 ( A ?/) = SCth^lo) <-0) 

s 2 SC/W/, /u")/) 

^fVm 0*^0 C- i (-iuii)) 


(nole;*ff? 


as a siil'scrii'l is always a positive intej'ei I 


.V 


Mcasm- 1 -Miu.’iU of Modal Power S|)eelral Density 
in I'erms ofCross Spet linl Iniiulions 

I'uKeiliiral Allematives 

There are a number of ways that can be followed to determine experimentally the 
modal psd of the pressure field in a thin annular duet, (‘hoice of method depends on 
several factors, including the use to be made of the linal residl and the availability ol 
data actiuisition aiul processing eMuipment. Two procedures are described in this report. 
This section gives a method that is essentially pallerneil alter the analytical procedure 
that was useil to derive an expression for the modal psd in terms ol cross spectral limctions. 

Section presents a different method, using a circumferential array of detectors incor- 
porating progressive time delays. In a thin annular duct, the array system maps the set ol 
modes having a common cutoff ratio and eliminates the reiiuirement for obtaining cross 
spectral properties. Additionally . the cross spectral method will provide a basis for com- 
parison with the delayed array measurement system described in Section .T.T.v 

Before proceding to this major topic, a brief description of the use ol the cross 
])ioperties method in measuring random pressure lields is lirst presented below. 


Modal i^sd I 'rom Cross Sped rai Measurements 

Since cross-properties are involved in this procedure, simultaneous information from 
pairs of microphones must be processed. II Ibe Held is stationary, these signals may lu 
obtained from a reference detector, and a second detector, separated A0- (T Horn the 
reference, with the imdcrstanding that the 0-location of both microphones must assume 
a succession of values. Atlernalively, signals from ^/ detectors spaced A0 = ai>art may 

be selecled for processing. I nis procedure is reiiuired if tiic Held characteristics are not 

stationary. 

I qualion (.C.Vl.'i) gives explicit instructions for computing in terms ol 

and aiul e<iualions (.C.^-20) specify clearly how the modal components 

and :"c to be computal in terms of the same aiullV^’J funclions. I herelore. 

it is sufficient to describe a proceilure for obtaining ami 5fTX. in terms ol measured 

data: 

l irst, (he 0 -averageil cross spectral deitsily is obtained, bor ei|Uisi'aced 

detector locations obtain 


- v|. 




B; -t 

} J 






wlu'iv SC^' is llu- iiu'iiMiivtl I’loss nil di'itsily lu'lwi'on iltii'flois 

III ^ iliul O ; +^07 . is t'v;ilii;iU‘it Im 0,-0, d7f/M , fk-. lo 

licl'iiu' ;is;isc( nl N liiiu fions of . V;iliii's o( :i'"l •1'"“' 

((T 7 , I'kiy I'o loiiiid ltii a si’( (»r rii,'(iiii‘iiL'ii.‘.s, ti) . in llii.' rani^v o! iiiknvsl. 

Typically for such a fivinicncy. sayLi4 . Hie known value of 5^ ilcleniiined liy 

Ihe above proceiliiie has been expressed by means cd ecpialion (d.3-13) in lenns o( the 
desired modal functions J^‘^i)aiul C” ^ 


co^yn(n + 

m-o ^ 


(3.3-22) 


The right hand side of ct|nation (3.3-22) isjnsl a I- on tier .series representation of (he 
known lefthand side (iinclion having (numerical) coefficients 

These coefficients are then found liy staiulard algorithms ft)r N-point 1-oimer analysis; 


M-l 

i (^) - -^'E Sff- C<Td ; ui) 0x171107 

/V tjo 

and 

v-< 

i- Y (.o^) ^ -it % S(T CtJj ; ^4.) Atit W<T2 


Cw^ o) 


In practice since StrC*^' is a complex mnnber ^fn is real, and 5T»*' is pure 

imaginary, it follows lhal^«can be obtained by proce.ssing cnily the real part ofS in the 
cosine algorithm and ( (he magintude of ) can be obtained by using only the 

imagimny part (d S in the second Ibrmula above. It Ihe complex (|uantilies. ,S, are useit in 
the algorithm, vanishing of the aiternale parts (d the (wo results will serve as a check on 
the cotnputations.) 

Summary' of Steps 

1. Uefore taking data tlecide on highest frequency. at which modal psii 
information is needed. 

2. ( aleutate (he ciivum feiential mode nmnbei at or slighlly above cutoff tor 


3. Select N, the number ol microphones or niicro|)lione locations, such llial U 
to avoid aliasing. 

4. Obtain data in the form of N/ cross spectral density functions l»elv.eeti 

miciophones pairs separ.ded by Q; (Oo^\ a/-() . tlnil is obtain : 


' . I> . ! . 

J ’ / / 
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NoU-: OnK N milosiK’.tni ;>iul nm'-li:ill o)' llu- ( spoc-tni iK'oci lu- nuMsiiivd- 
lilt- ivinniiiiii)' i rnss s|K‘i Ii:i convspoiuliiu', lo valiios CT >7f aiv nhlamahK- as 
i-omplcs cDtijiijiak's of Ww \\n \wr siiuv 5 (&j ^0^ ^ j u?) - 

5. Avi'ia|!i,' lliosi’ I'lmctioiis nvi’i' © by i'oniiula (3.3'“2l) • 

b, Sok'ct l iviiunu ios O) - ul ,10, u?„y , ami proa-ss tlio ivMills cd’ 

step 5 at oach ol‘ these w> by eipialions (3..V23) lo obtain (he limelions 

7. ('ompule the Sji (i'^O ;>'hI 5^ of the nuKlal psd. 

Iroin the sum ami dilTeieiiee of 'a ami '6 obtained in step b by the simple 
fornuila (3.3-20). 

3,3.3 Delayed Cireular Array Systems for Measurement of Discrete Freiiuency and 
Random I’rcssiire Fields 


Introduction 


A prime objective of the background material prc.scnted in section 3.3.2. “l ield Structure 
of Random Noise in fliin Annular Ducts Uackground" has been lo define the modal power 
spectral density function, its properties, and a method lor its ineasurement. 1 his inlorma- 
tion will be needed now lo establish the operating characteristics ol the delayed circular 
array for pressure Held measurements in a thin annular duct the subject ot this mam sub- 
section rite background matter has formulated the concept ol modal power spectral den- 
sity S (m. to), and itscomponentsS + (|co| amt S,;,(lto|) I'his function, shown in 1-igure 
3.3-2 for an illustrative random field, gives the distribution over circumterential mode num- 
ber. m. and frciiueiicy. to. of the mean-square pressure. Only integer values ol m apply 
because of the requirement for circumferential pre.ssure periodicity. Associated with each 
value of m is a continuous distribution over fietiuency. Parenllietically. a cohercnl. discrete 
frequency field may be characteri/ed by a corresponding modal psd in this case, the 
frequency distribution for any mode is m>t a contiimoiis function but rather consists ol one 
or more “spikes” or (effective) 6-funclioiis. 


The methods used to evolve the analytical formulation of the modal psd. S (in. to), in terms 
„f the mathematical model for the pressure iK'Id. P WA) involved cross correlation and 
|•■ourier transform operations. The resulting expressions lor S tin. w) and its componeiils. 
S+ (loji) and S,;, tlcjl). e.iualions (3.3-1 7). were given explicitly m terms ol cross spectral 
fuliclions. Ill conclmling this backgiouml. section 3.3.2. details ol one procedure lor cal- 
culatiiig S (III. lO) from ineasurctl values of cross spectral limelions til ilata olUained expeii- 
mentally were specificil. In essence, S tin, to) was mapped as a function of the cooidmate 
pairs m and to. using cross spectral fiinctituis of pairs of microphone sienals. 

Now in the (111 to) plane over which the modal psd is defined, lines of constant slope, 
to/m correspond to a common value of cutoff ratio, oi equivalently, circimilerenlial 
spin velocitv il. bu wave spiral angular tlirection, as illustrated in I igures 3.3- 1(b) ami 
M- ' It should be possible to evolve an alternative method for mapping S tin,, to), using a 
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A drcumtVrontiat array of micropliones, incorporating a sdcclabic, progressive time delay 
between elements provides just siidi a measureineni system. Hy adjusting llie time delay 
rate, the array is tim.ed effectively to follow all waves having a common cirenmferential 
velocity or spin rate, 'riuis the array acts as a filter in mode-freiiuency space to enhance 
tlic set of waves having the common properties of the .selected spin rate, spiral direction, 
or cutoff ratio. 


An outstanding feature of this array system is the complete elimination of the need to 
execute any cross-properties signal processing operations for each array rate setting a 
single, onlinaiy frequency spectrum is used to obtain values of the modal psd spinning 
mode components at the selected spin rate and the overall power of the set of these modes. 
Another desirable feature of the method is that it applies equally well to the mapping of 
discrete frequency pressure fields. 


As a guideline to the array analysis that follows, the material is organized as follows: 


Section 


Features 


(1) The Array Response Function Using a single spinning mode discrete frequency 

field for convenience, the main features of the 
array response arc iletermined. This case is 
preliminary to analysis of array response in 
ratrdom fields. 


(2) Array Application to Discrete 
Frequency Field Measurements 

(3) Array Response to Random Fiekls 


(4) Application to Distortion Tone 
Measurement, and 

( S) Application to Mroadhand Noise 
Measurements 


An example of array use in mapping discrete 
frequency fields is given before proceeding 
to random field.s. 

Hie power spectral density of the (single) array 
signal is expressed as the Fourier transform of 
the signal autocorrelation function. This result 
is a combination of the pressure field character- 
istics anil the array response function. Using 
only the simpliest property of the response 
function, a method is given for determining 
the modal psd of the delayed array system. 

A basis is found for the existence of random 
spinning modes as a physically meaningful 
concept. I iirtiier properties of array response 
and its use in random fields are deferred to 
the subsection ’'Details of Array Function”, 

These two subsections, “Application to 
Distortion Tone Measurement" and "Applica- 
tion to Broadband Noise Measurements" pre- 
sent essentials of how the array method may 
be used in ineasineinent of two random field 
cases. 


(6) Dctiiils of Array Response i-unction A more coinpi’eliensive examination of the 

response function is given in order to allow eval- 
uation of the sources of error in array measure- 
ments. Possible means for sharpening the 
response are mentioned. 

(7) Hxtentions to Cieneral Cylindrical Mention is made of the po.ssiblc uses of arrays 

Ducts in general cylindrical rather than narrow annular 

ducts 


Ihe Array Response Function 

In this section are presented the basic features of the performance characteristics of the 
delayed circular array when deployed in a narrow annular duct and subjected to the pres- 
sure field of a single spinning mode. This simple discrete frequency single mode input is used 
to determine what amounts to the filter characteristics in (m, to) space of the array, and 
corresponds to the practice of analyzing ordinary (to) filters under harmonic excitation. 

This procedure is a helpful preliminary to considering array detection of a plurality of 
coherent modes and the detection of random fields - the main objective of the investigation. 
After finding some of the key properties of the array response, a more complete presentation 
of its features is deferred until section “Details of Array Response Function”. 

Two purposes are served by this break: 1 , use of the array in measuring a discrete frequency, 
coherent field containing a plurality of modes and frequencies is illustrated in the section 
"Array Application to Describe Frequency Field Measurement”. Then the next section 
gives the derivation of the array response to random fields, followed by illustrated applica- 
tions in the next two sections. I'he material in these sections, “Application to Distortion 
lone” and “Ajqilication to Broadband Noise Measurements'’, requires minimal information 
about the response function and is. tlierefore, presented immediately after such information 
is obtained. This procedure should help to convey the essential concepts better than if 
more details of the response function were given here. 

After the basic principles and usage of the array method have thus been presented, a more 
detailed examination of the array response function is resumed in section "Details of Array 
Response Function”. 

The basic features of the system may be perceived by reference to Figure 3.3-3. which 
shows an array of N detectors equi-spaced AO = &?r/N around the circumference of a 
narrow annular duct. An incident wavefront traveling in a general direction, a, is shown, 
together with its phase velocity \q = 12 in the 0-dircction. Note that, since the concept 
of wavelength is not even presented in Figure 3.3-3, all wave trains, mgardle.ss of frequency, 
that spiral in the direction a have the same circumferential phase velocity or spin rate, 12. 
Further, at any or all frequencies it can be seen that there is a one-to-one correspondence 
between direction and spin rate. In addition, cutoff ratio of a wave train is also uniquely 
related to or or 12 since it equals the tangeniial Mach number, 12r/('. 
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ARRAY DELAY RATE = At/Afl = 1/n = ?? 


Figttre Delayed Ciniilar Array Developed I7fu' 

Imicpendcntlyof tlic pressure field existing in the duct, the array system is designed to 
incorporate a progressive time delay for detection of a signal by successive microphones. 

(This time delay is usually best provided in the data reduction phase of the array system 
operation). Tlie array delay characteristic is adjustable and is conveniently specified by a 
parameter, , called the array delay rate. The selectable delay rate is given by 
seconds per radian. Tims, for a selected value of , signals received simultaneously by 
all detectors are delayed in processing by time increments: O , ^^0 , 2^^ " " ’ , / s 

... CA/“|)t7 ^0 corresponding to the microphones at <9-C , 

d. progressively delayed signals are adtied to form the array signal. 

Returning now to l-igurc 3..T? and in particular to the illustrated wavefront trading in the 
Q -direction at JfV radians/scc, it is easily seen that il the array delay rate, ^ scL/radian, 
is selected according to the expression , the delayed array detectors cflectivcly 

sense the wave simultaneously instead of progressively as would be the case lor /.ero delay, 
Thus the array is atljusted by means o\'y[ U) “follow” the wavefront. The array signal, 
being the sum of such progressively delayed detector signals will, therefore, enhance all 
waves traveling in the direction (K (and hence having a common spin rate, JL ) compared with 
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otiKM- wnvi’s tiiivdiii}* in (lilTcicnt liircctions or with spin rntiis (lilTcrcnt from M. . Notice 
that tliis enhancement applies to the set of all wave trains (imiepciuiently of frequency t 
having the spin rate JLr , in a complicated pressure field slnictiirc. 

It should therelore he anticipated that in a complex pressure field (regardless of pcrioilic or 
random time dependence and independently of the degree of modal coherence) that the 
ideal array signal overall mean-stiuare value is propcrlional to the total acoustic power of 
all waves traveling in the direction corresponding to Jl = . hurther, it may be expected 

that an ordinary narrowband analysis of the array signal will give the spectral distribution of 
the acoustic power propagating in this direction. 

These expectations are indecil realized, as will be seen in the following material. At this 
point, it should also be made clear that use of the delayed circular array completely 
eliminates the need for employing cross-correlations or cro.s.s-spcctral signal processing 
operations. Thus the need for executing a large number of relatively complicated cross 
properties operations required by the method previously described in section 3.3-2, “Field 
Structure of Random Noise in Thin Annular Ducts”, is completely eliminated in the 
determination of the modal psd function of complicated pressure fields. 

The remainder of the report on this subject establishes the validity of the above-mentioned 
array properties. The key questions to be resolved are: How well does the array enhance 
the desired set of modes or wave trains having n common spin rate, Jl~y^ ? Also, how 
well does the array suppress contamination of the signal by waves traveling in directions 
differing from the “target” direction, oC ? 

Answers to these questions (which are really two aspects of the same question) are provided 
by deriving the army response to a single, coherent, discrete-frequency spinning mode. Such 
a mode plays the counterpart role of the pure liarmonic signal used to define the transfer 
or window function of an ordinary (time-series) filter. For, as will be seen, the time-delayed 
circular array is in essence a filter in the (G,t) or domains. 

Consider a single TH -lobe pattern spinning at A radians per second in the 0 -direction, 
lire pressure field will have harmonic variation at frequency u)^ w-A. and can be represented 
in conventional notation by 




(3.3-24) 


where the complex number Cm 6'ves the amplitude of the mode and its phase wilii respect 
to an arbitrary reference harmonic signal. 

Let the array consist of N detectors spaced AG =2"/^ . With a delay rate W 

the pressure at the T)^ detector located atGn=^^^ will he sensed at a time '2^ = 
later than by the detector at 0«> s O .If the pre.ssurc at On is 0^= the delayetl 

signal will be + . giving 




(3..V25) 
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(3,3-29) 


Now let the army delay rate be set to the value 7f * w/<*^ . This makes 

and makes each of the (unit) complex terms exp(in^ 2 ??/w ] equation in (3.3-29) a unit 
real number. Their sum will be A/ and 2 is thus exactly un'ty, Consequently, from 
equation (3.3-26), the array signal reduces to 




Tluis, as was explained previously in a qualitative way, the array signal is identical in 
amplitude and phase with the complex modal pressure. Cm the array is tracking perfectly. 

But suppose we leave the array tracking the mode Cm undisturbed and add to the pressure 
field a new mode C|< at the same frequency, u) . Wiiat happens to the array 

signal as a result of superposing the new mode? The answer turns out to be - absolutely 
nothing. The array continues to track Cm perfectly as before, and its signal is completely 
undisturbed when the mode is superposed. 

This complete lack of signal contamination by a mode different from the target mode (and 
thus spinning at a rate different from the target rate ) is exactly what had been 

hoped for. To see why there is zero contamination, examine equation (3.3-29) for ^ . 


The terms are all unit complex numbers with successive angular orientations in the complex 
plane differing by radians or p/n revolutions. Now ^ was set to make 

~ m of the target mode. So ^ for the new mode will be . If M is 

one more than tm or one less than , it is clear that the terms in equation (3.3-29) 
consist of N complex numbers in asymmetrical “wheel spoke” arrangement. 

llic sum of this arrangement is obviously null and the average, which is the array factor, 
is therefore zero. Hence the array signal response to the added M -lobe pattern is null, 
and its response to the original target w -mode remains unaltered and uncontaminated. I'or 
other values ofM , it is shown in section “Details of Array Response Function” that the 
terms compri.sing Z also make symmetrical, null-sum spoke patterns, provided that N ^ 
the number of microphones in the array, is “sufficiently” large. 


This elementary analysis explains the previous statement that the array enhances the response 
to the target mode. If the acoustic pressum field consisted only of the target Tfl -mode, 
the array signal (before averaging) would be N times the output of a single microphone. 
Consequently, local contaminating pressure Huctuations. due for example to hyilrodynamic 
pressure nuctuations of airflow over the detectors (“noise”), are very poorly correlateil. It 
is easily shown that the signal-to-noise ratio of the array is also erhaneed in this ease by a 
factor of 'j/Ji/. 
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It IS W<>nl) mi>iuioninB that, for the case where several modes exist, all at the same frequency 
the ah. ity of the circular array to reject totally all modes except the target mode depends 
trititally on the fact that the wave number '»i in the quantity d cJ-^) is restricted 
by physical considerations to integer values. In conventional linear array theory there is no 
corresponding requirement that the number of cycles of signal variation along the array at 
any instant be an exact integer. Consequently, the equivalent wave number can vary con- 
tinuously and waves arriving from other than the target direction cannot be completely 
prevented from contaminating (lie array signal. 


Thus the success of the delayed circumferential array system depends both upon the system 
design per and also upon its use in a field having inherent pressure periodicity in the array 
dircc ion. While this result has been obtained for the coherent modes at a given frequency 
It will be shown in section “Array Response to Random Fields” that the same integer 
restriction on implied by ^ -periodicity is an important element in the use of the array 
in measuring random pressure fields having such periodicity. 

The next feature to examine is the ability of the array to track not only single or multiple 
modes at a given frequency, but^so an entire set of modes having a common spin rate Vl 
T o do this, the array delay rate 7/ is made equal to When the array signal is 

spectrally analyzed the spectrum levels at frequencies^yt. . ^Jl, SJI, . . . correspond exactly 
to the amplitudes of the successive modes , 2 , 3 - .This case is examined in 
more detail in section “Array Application to Descrete Frequency Field Measurement.” 


The last matter to be examined in this section is the problem of mistuning. That is, 
suppose the pressure field consists of the pair of modes characterized by and at a 
common frequency A7 , as before. The delay rate ’’^remains set to track the An -mode 
and the array signal is exactly . Now suppose that the frequency is changed by 

an amount ZtA . Two effects result: 


1 . The array response to the target /n mode is slightly reduced by the drop in the 

value of the array factor for the target mode, Z (^j.Z y) 

2. Secondly, a slight amount of contaminating signal is introduced into the array 
signal from the off-target yt-/ mode. This contamination results because the 
array factor for the /f- <^7 is not exactly zero unless 
The amount of contamination depends not only on the magnitude of 

but also on the proximity of to the target mode number 7n. 

These effects are disclosed by the graphs in Figures 3.3-4 and 3.3-5 which present the 
magnitude / o\ the (complex) array factor. /«?/ is called here the array response function 
Figum 3.3-4givcs/.?/ versus . When (on target l./r'/ is I , and Figure 

3.3-4 shows that for-^^c/-^ /^/ differs trivially from unity. 

nie/V-mode contamination is harder to understand. Figures 3.3-5(a) and 3.3-5(b) show 
/^/ as a function of the separate parameters /m and uJ for a fixed value of ^ . In 

particular, Hgure 3.3-5(b) shows what happens to/V /wlien the target mode fretiucncy 
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dumbos hy ;i stnall ammint, A , Caivru! I'xaminalion ol' tliis liniirc shows tlie coiitam- 
inatiii;', contribiilions Irom modos in ihc nd^Oduirliood ol tho tarj^ol iiioilc, ni. I'Or ui ■■ (Oju 
t Aco, llu‘ laincsi foiUamiiialioii ivsiilis iVoin mock* m+l , wilh siitvi'ssivdy lessor amoimls 

Irom inodes m i .1, m+d and m I , ni 2, . . , when co - Au), tlie ureafcst eonlamina- 

tion is prodiieed hy llie m I mode, wilh proc.ressively smaller eontrilnilions I'roin more 
distaiil mo«les. 


The properties of the array response I'unetion, /t? / , revealed in tliese ilniires are derived in 
section “Details of Array Kesponse !'nnelion“, Uel'ore takinn np tliis detailed snbjeet, the 
array response to a random Held will be given in section "Application to broadband". 
However, in section “Application to Distortion Tone”, an application of tlic array system for 
mapping coherent, discrete frequency field is given first. 

Array Application to Discrete rreciuency Held Measurement 


The delayeil array system can be used to map tlic complex modal coefficients, 
of the pressure field of the fundamental and liurmonics of fan blade frequency noise, 
o.g. 


P = Ms. ^ ^ ^ 




where C^o is fundamental blade-passage circular frequency, t7 - 3 is the highest harmonic 
index of interest ami is the greatest propagating wavenumber at this harmonic. 

Values of /V ^ for and n--3 will be decaying modes, but such modes can be 

detected if their amplitudes arc sufficient. Also note that negative/’’ subscripts are now 
used and identify reverse modes sinee cu~ /?4t^is taken with a negative sign. (This usage 
is consistent with conventional spinning mode formulations.) 


A typical modal amplitude (not power) spectrum is shown in figure 3.3-6. Also shown is the 
first target delay rate setting, This ratctargels the modes: m-/ fortes; u/^ 

A?7- 2 fox foxuj- 3 uy^ \ 


As has been mentioned several times, the number of microphones, /V, in the array is 
.selected such tliat A/^/rTwX u/=3u^ to avoid aliasing, riicn the following steps are followed: 


1 . 


Set 



■ 321'^ u. 
— CU iyy'o 


a. niter the array signal at oJo 

b. filter the array signal wXzuJo 

c. filter^fi^J the array signal at 3iPo 


giving ZC/// 
giving /o^z! 
giving / Cj3l 


3. 


Set i‘"Tel 2 


for 


lor <? /’7- ^ /aR. S OJ ^ 


4. 


a. filter array signal at 

b. filter array signal at zu/q 

c. filter array signal at 3uJo 


giving ! C^il 
giving /C4Z( 
giving / C «3 / 


,S. Repeat tlic proce.ss for other ilelay rale settings, including negalive rales for 
reverse iiunles to completely map IlieZ'^j a/jllekl. figure .V.Tti can be used to 
see what values of 'y '^ou required. 
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DECAYING MOOES 


REVERSE SPINNING MODES 

t I 

I f 



FORWARD SPINNING MODES 


n = — FOR CUTOFF = 


FIRST TARGET SET 
1 m 1 2 3 

n <jj Wq’ 2wq 3(a)q 


Figure 3. 3 (i Sample Mode AmpUmde Spectrum of Fan and First Target Set of Array 


Notes; 

a. If the random broadband noise spectral density is “low”, the effective post- 
array filter bandwidth need be only narrow enough to resolve components one 
octave apart. 

b. If the random broadband noise levels are high, the post fitter bandwidth must, of 
course, be reduced. 


c. 


If the array signal is low-passed so that frequency information beyond J<^^> 
is suppressed, the signal total mean square value is the “power” of the set of 
modes traveling at jt.- For example, for the target set illustrated in Figure 
3.3-6 the array signal mean square reading equals ^ f C„ • ^ C3I>Q3iy 


d. Tlie relative “phases” of the modes at any one harmonic, U/Of oa.3U/o 
can be found if reference signals e-/ - and <2 j - cos 

arc provided and if the Held is perfectly coherent (frequency, amplitude, and 
phase of each mode are substantially constant over the analysis time). This is 
accomplished by the usual methods for phase determination, operating suc- 
cessively on the filtered signals at a given array setting corresponding to the 
modes/rj.^^. eWo -y) . 77 • The fact that the complex array factor re) 

is exactly real and unity at these cu) combinations makes possible this 
phase detection, or equivalently, determination of the complex modal co- 
efficients, t/mm as well as simply amplitudes / C/r>/?/ 
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i‘. ir ilii-iv is soiiu’ sil'iiiil JilliT/ 1 ^ ' 

nimlomly v;iryii4i i-ompomMil , IIk' ;in:iy iiiiiy lu' iisisl in al Icasl iwn ways. I Ins 
sij’iial iillt-M may 1 h’ liiu' )o Ian spi'i-d lliu lnalions, or il m;iy ivsnll lioni iiilliiif; 
of npslivam waki's llial havi* soiiu* ilopivo ol unsU'atliiu’ss. 

1-1 . I'irsC irtlh' ami)lilinK-sX-m/7/il‘>"‘' wanti-d, llu'iv is no pioldmn 

tlu.' post nik'ivd array will lU-loil coinponoiUs havinn rn'iiiK'iu-y liamlwidtli 
small compaivd wiili ludli tin.' array (!iV(|Ui‘niy ) I'aiulwidlh al a parlii.’nlar 
'9 and to tin.' ptisl Ullm clT(.’(.'tivc Inmdwidll!. I'iiiinv d..^-5(l)) shows how 
small rivciuoncy sliifts ivsnit in only slight loss of array ivsponso and small 
amlamination from nciglihoring nnuK's. SoLlion “Dotails ol Array Response 
I'nnction” explains tins ’‘on* target” response in terms of detailed properties 
of the array response fimelion. / <?/ . 

e2. Second; If the fan sj>eed is constant, but signal jitter occui-s. phase as such 
has no meaning and only the amplitiules/CAr?/? / can be loutul. Howevei. 
if frequency changes result from deviations in fan speed, the portion ol the 
field that is coherent with respect to the rotor may be determined by 
sy nchroiujus iletection techniques. Here the releience sign.ils G /~) - 
cos nC (^) are generated by a rotor shall transducer 

that follows shaft speed or fundamental blade frequency. /t 

By sampling the array signal ft") at time inienals controlletl by 
the instantaneous value of <co, ^CJ«) the coherent portion of the array 
signal is enhanced. Processing of this synchronously detected signal com- 
lH>nenl by the methods described above will give mean values of both 
colierenl modal coefficient amplitude and phase, again assuming that the 
range ofA<^^> is not loo large. (A more sophisticated procedure wouhl 
involve using a variable delay rale ^^’t.) 


Array Ro.spoiisc to Raiulom PieUls 

This topie is the major item of concern in the investigation reported here. It will be shown 
that the delayed circular array inoviiles a relatively simple and effective system lor mappikc 
the modal power spectral density function. which was established in section d.d--. 

"I'ield Structure of Random Noise in I hin Annular Ducts”, as an appropriate lunction to 
specify the characlerislics of random fiehls. 

Although the linal result of this study a method for use of the array system is simple, 
the derivation of the arra\ response when deployed in a random field is somewhat compli- 
cated. I hese comi^lexities are just a consequence of having to employ cross proiierties 
analytical techniques I not measurements) and of having to use some ol the l.tirly intricate 
re.sulls of the background malerial in section 


Despite the comi>lexities ol executing the requiied algelua. the analytical method can be 
specilled in one senleiice. the analytical model of the pressure field given by etiualion 
t,V.k2). / ‘‘n m is used to obtain 

the array signal, alfer which the mmlal jiowei spectral density ol the signal is obtained by 
1-ouriei Iraiisfonniiig the signal auto corielation timctioii. I he following material results 
from compleliiiii these instniclioiis. 


s: 




As hi Ills- iiivviiihsly iiniilyml ilm-rolo In-.iiu'iK-y iw. Hu- array silinal is rJvi'ii by 

M’l ^ 


where r,7 -• isilu array ilelay rate. 

The eorrespoiuling auto-correlation funetitni is 

■7^ 

f-i 

Tliis can he reduceil to 


(3.3-31) 


/€ /T) $1 ^j< (Z^£ ^-Jj 7 ) 

J‘0 ar=0 


(3.3-32) 


where T) is t!.c ordinary cross correlation hetween the undelayed pressures at ^ 
and^^ ■ 

The power spectral density of ^ A ^7 ^ * is ‘•'^''sfonn of 

equation (3..T33). It is known tluit if the transform ol a lunction ot t 
d) , then the transform of the function with argument »s ei ^ > ■ 

us with C. - £ z< Jj fj 2 Tf^ /^ ^ 


(3.3-33) 


Thus 


^7 7^; - 




/ 





yr 


£! ^ ^{Ku)u/^2.y^ 


(3.3-34) 


✓ I 


onr.TKAi. ® 
HI.- pnnll ' 


s.t 


wliori' i1k‘ lioss spi'ctral iK’iisity hi'twi'on Hu* iiiuldayt'd sij^nals at ;iiut 

6^ 1 is the IraMsIorm ol ^ jti/'T). l'.(|iialioii (d.3-34) sliows that the singk' output 

psd of the array automatically “p^’rlorms'’ the operations of mcasiiriup /V 

separate cross spectral densities. adding them after appropriate phase shifting 

hy the factors 

Let tile /V^ terms in etpiation (3.3-34) be collected into groups for which tlie combination 
K j is constant and is denoted by the index, -b . Also, for-^,<'/<x<^ the cross spectral 
density between detectors at and put ). suppressing the ^^-depcndence 

for simplicity. Tlien (S^^can be written as and 

be: 






Collecting terms for which S is constant, and distinguisliing between those with positive 
and negative exponentials, the A' terms can be arranged as follows: 




A' 


A ( 
_5.^3 


_/ 


U-^ 

AJ 




(3.3-35) 


'Iliis ratlier complicated expression is a mixture of spectral functions of the pressure field 
and components of the array factor cf . One reason for the complexity is that it 
involves a double summation wl-icli in turn is a consequence of representing a power spectral 
density rather tlian an amplitude spectral density. 

The form of equation (3.3-35) contains fragments of (finite) L'oiirier transforms on the 
twoj -sums of the field cross spectral functions. If the two exponential weighting factors 
happened to be numerically identical, the expression would be greatly simplified. This 
simplification can in fact be achieved if attention is restricted to integer values oi the 
quantity . Invoking this restriction will be shown not to impair use of the array 
method in measuring the modal psd. 

Operationally this restriction means that for a given array delay rate setting, the 

array signal/>s</<S'->)fu;j will be evaluated or used only at the successive values ol ^ • 

It is important to recogni/c that this procedure ensures that 
the array factor t? is exactly unity or zero since is zero lor and 

(3-- A tion-zcro integer for -^ / Jr 'fhe use of3y/^<.^;in evaluating -5 //t 7, oz; at 
intermediate frequencies will not be discussed at this lime (see section ‘ Details ol Array 
Response l imction”). l or the present consideration will be liinitctl to obtaining -5/rn,u/) 
from ) for this specific value (d . Accordingly let an integer. 
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'I’Ik'ii till’ SL'cniul cxpniK’n'.iill ill t*(|u;ition ( 3 . 3 - 35 ) hoioiiu's 




<r 


- TT/y^ 


<r . 




wliich is just the factor of the first set of terms in eiiualion 3.3-45. Accoiiiiiittly 


JS-yf^r^) (3.3-36) 

Tile quantity in braces is simply tlie (finite-element) space average of the cross spectral 
liensity of detector pairs spaced (T* -S^;^apart. This (juantitv was denoted in section 3 3 
equation (3.3-10) as Hence 


^ i i*') 

_S-cP 


(3.3-37) 


Replacing the above _S -sum average by the exact equivalent integral representation (ref 1 2) 
gives; 


^ (3.3-38) 

-77 

In tliis form, -5'7/'=^ycan be recognized as the I'ourier transform of S^/eT , ^ ) with 

the transform variable - replacing/n . Since the transform of 

was found to be the modal power spectral density of tlie random pressure field, the'link 
between and-S/^vzi./'^has been made. 

-S^/o/Jis obtained m terms of by replacing with -M in the previously 
obtained expression for -S^ z^/^equation (3.3-17) g-’diig 




(3.3-39) 


l•.valuatmg ei|uatioii (3.3-39) for the positive and negative values of , namely AY‘/ay/yl 
' contributions will result from (y^ and y only for '■ 


resulting in 


no 


i 


4 ‘^'1 H"^/) I'/y/jj 
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ORIGINAL PAGE IS 
OP POOR QUALITY 


i iictlly, siiK’i' hy LH|iiatioiis (.V3-20) iho fni wiml -lobe mode power speetra! density 
eomponenl j is piveii hy^ i/^C-/uj/) or < IIhtc follows: 






(3.3-40) 


Now Sj^ ( Itol ) was seen in section 3.2.2 to be the component of the modal psd, S(n, w), 
associated with values of n and cj having opposite signs. It was described there that the 
choice of notation, involving + and - superscripts was based on a formal correspondence 
of forward and reverse spinning modes between the random and coliercnt field mathemat- 
ical results, and that it also anticipated information that would be (and has just been) 
obtained in this section. S„ (lwl) from equation (3.3-40) is a term in the expression 
for the array signal psd. The array delay rate ^ has been positive tlirougliout this entire 
derivation and it was seen in the previous section on arrays in discrete frequency fields 
that this selection will cause the array to track coherent forward s, winning modes. We 
have just found out tiiat with positive delay rate the array in a random field tracks only 
the portions of S(n,co) with opposite n and to signs, which were called the forward spin- 
ning modal components, Sjt (Jcol). Because the same array in a coherent field with 
positive tracks only forward discrete frequency spinning waves, it ^merges that the array 
in a random field recognizes values of sj (|oj|) in exactly the same way as it identifies 
colierent forward spinning modes. 

S+ f|w|) is the power spectral density of the forward .spinning mode with circumferential 
wave number n, and Sj, (|cj|) is the random field counterpart of the modal power of the 
reverse spin n-lobe pattern in the coherent field case. Tlierefore, it can now be stated 
lliat spinning modes ilo “exist” for the random field case in the sense that the array 
provides an operationally definable way of obtaining their power speetral densities which 
comprise the field modal psd, S (n, to). 

Returning to equation (3.3-40), the array signal p.sd, S_(to), is the set of the numerical 
values of successive forward-spinning modal components to ), each evaluated at a 
specific frequency |tol =fn/rj| . This is the equivalent of the class of forward modes 
having a common spin rate Jl. = % = . Values of the array psd at other frequen- 

cies are undefined, since the derivation of equation (3.3-10) specifically restricted the 
results to integer values of toTj. (By similar means it can be shown that if the delay rate 
is made negative, the array will track the corresponding set of reverse spinning modes, 

I'iguro 3.3-7 shows the sample random field spinning modal distribution illustrated in “I 'icld 
.Structure of Random Noise in Thin Annular Ducts.” section 3.3.2. In tiic (m. oJ) plane, the 
line with slope I] m/w represents (he target sefting of the array. At the inteisections of 
this line with the lines m - I, 2. 3 . . , ordinates of the component modal psils. s| (|u)| ^ |2 /t?|). 
etc. aie shown. I’he set of these ordinates is (he meaningful part of the array signal psd. as 
given by etpiation (3.3-40). 
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^ CirriiiliJ't-iviiiiiil Arnty rtir!:t-t Set I'Sli in l<iimh>in l ichl 

In pracliLV. the valuos of llu* array psj Sj^icJ) at llu' spccil'ic IrctiiicnciL’s lo - • 

etc. arc ohtaini?d simply by compuliiiji or nu'asmini; iIk’ siiinal psd at Ihcsc fri.'iitii.‘iKii.'s" 
only, ignoring inti'niicilialc lirqucncy inrorinalion. As will ho soon in section “Details of 
Array Response l•llnction’^ tlie inlernieiliate rieqiiency inrorinalion is eontaminaled by 
array signal contribniions I'roin neigliboring. olT-largel nuules as was seen lo be Ihe case when 
measuring coherent Helils. Ibis eonfaminalion is a conseqneiiee of tlie incomplete vanisli- 
ing of the array I'actor /. lor bordering modes at I'reqneneies in the neighlnn hoods of the 
target rreqnencies. ) 

In the subsections immediately below, two exami'lcs are given ol' the use ol' the delayeil 
circular array system in mapping the^;^ spinning mode components ol' the 

modal psd lor random llelils. 

Application to Distortion Tone Measurement 

The term distortion lone denotes narrowband eneigy in the neigliborhood ol' rotor blade 
passage treqnency Tor esamide. due lo Ian cutting ol' unsteady inllow disturbances, l ei 
|CU,| be a circular rrequency in Ihe neighborhood ol blade I'assag.e treqnency, lull not 
esaclly eoineideni with HI’l . The array ssstem can be used It' measure the imulal com- 
pt'silion of the lieU al)cO»/ ui the rollowing. \\ay: 


ValuL's or tlio power spectnil ilensilies of all propapaliiifr I'orwarJ ami ivveisi* spinniii!’. 
iihkIos L'ontrilnilinp lo the sipnal at I’retiueney |a)«| are re(|Uireil. Ixt the array sipnal 
he IHtered (digitally or eleelronieally ) with a narrowhaiul litter centered at It will 

he convenient to suppose the niter elTeetive bandwidth to he unity; if some other hand- 
width is used, division ol' the llltered signal power readiii|! hy the actual haiulwidth will give 
the power spectral density at |cOol ■ 

(In this and the Ibllowing suhsection, the same remarks concerning selection of a sui'llcient 
numher of microphones in the array to satisfy the Nyquist requirements apply as well to the 
random Held cases as were made in connection with the discrete frequency fields discussed 
in section, ‘‘Array Application to Discrete t-'roquency Kield Measurement”,) 

Since equation (3.3-40) for the array signal psd applies only when cor/ is an integer, the 
delay rate, 'p , must he selected to correspond to successive mode numher n 
That is, to determine the power spectral density at |£0»l of the forward mode, Sj 0^.1 ) ^ 
lw»lip s n requires ips"/l«J*l. Reverse modes are detected by setting corresponding 
negative delay rates. Tlie following tabulation indicates the delay rate settings and the 
resulting output of the filtered array signal: 

Freuuency ltJ«| (Array signal filtered at )( 0 ol ) 
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Delay rate 

V 
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Modal psil 
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Reverse mode 
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1 

3 

0 

Delay rate 
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V|o/u| 

Array output 

■ 


SlC'Wot) 

S] (|Wo|) 


Modal psil 

SCrt1,(j) ■■ 

S ( Ij Wo) 


S(3, U>o) 

S(o,w.) 


Uy repeating this process for a set of/w*|in the neighborhood of blade passage frequency, 
a map of the modal composition of di:stc'rtion tone noise is obtained. In performing this 
data analysis it is assumetl that for the nearest to Itl’l’ the filter effective ha.iidwidth 
is sufficiently narrow ami the rotor speed steady enough so liiat coherent HI’I' signal is 
excluded from the filtered array signal. If such is not the case, it may he necessary lo 
obtain a reconstructed signal in which (he periodically sampled coherent blade passage 
component has been eliminated from the total signal. (If a discrete frequency component 
is conlaitied iti the filter battdwidlh, the array will of course respond to this signal as well 
as to the haml-litniled lamlotn power.) 




Appiicalion lo Hi'oiulhaiul Uaiulom Noise Measiircincnt 


TIk'IV arc applications where knowledge of the }>eneral shape of the modal power spectral 
density function, S (r'ljio ) for a broadband random noise field is desired. An alternative 
to the methods based on cross spectral measurements ilescribed in .section .^,3.2, ••|■ield 
Structure of Kaiuiom Noise in I hin Annular Duels”, is provided by the array system. 

This procedure once aj^ain is based on the property tliat the array si(!nal provides a corr- 
ect measure of the power spectral density of each mode only at those frcciuendes for 
which 0\ is an intejier. 

For efncieiit mapping of the modal power spectral density the delay rate rj is set to a 
(small) value, say ip, . The re.sulting array signal, S|p^(ftdO is narrowband filtered or FFT 
transformed at the following seipicnce of center frequencies, giving the indicated modal 
power spectral densities: 

Delay rate ~ 7 , Array signal narrowband filtered over 0 ■< lvo| < 


St 

S(-n, V^,) 

^ Yi, ’ vs *“t' V.I ^ 

cor»‘esponding to the value of XL for which cutotf occui's. Negative values of rj produce 
the reverse mode spectral densities. 

Figure 3.3-7 illustrates the process for sample value 1 ?. It is dear that this procedure is 
an alternative mapping of the Sj and Sm components of S ) where values of S 
are determinetl at the intersections ol lines of constant with grid lines at integer 

m in the (w ,10 ) plane, as shown in 1-igure 3.3-7. 

It may be noticed that the mode So ( IWI ) has not been specified in the.se tabulations. 

If the system described above were extended to the case n = o , the corresponding freq- 
uency |t*>l = would be lw| ~ o . In practice, information about the spectral content 
ol the 0=0 mode is obtained by setting the array rate ^ at zero and simply perform- 
ing a narrowband fieciuency analysis of the basic array signal. The case of 5=0 corres- 
ponds lo enhancement of all waves arriving simultaneously in the )dane of the array 
(plane waves), iiulependently of fret|uency. 
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Details of Array Response Inmctions 

„o,e,U »a'un,cy of .l,c array sod, » ■ 

r':;;’arc rindX:: „ nra, ,.c notary .o yo 
IZS Xk nlydiraic detail bdorr. au inter, nv.alU.n of the reaulta can hr made. 

|.„r eenvettienee the ,rrevi„.,dy obtained eameaaion for the array fartor. /., is re,.eated here: 

I tnfl *«*/m 

-z -- N 2^ ® 

n< 0 

uik&re & - (w “ 


The properties of Z can be found by considering two eases: 
Case A: /fl = X .an 


,n „„a ease. Z is generally the -b" ^ , 

I rZ. NOW for 7 : S; h-h. ait such — 

(Subcase 1) — 

The nrst hrteger value 1 = O is espeeiahy signifieant. I;;_>dis ease, o 

results when ’’“J ^ for all modes having the ratio «/..>* /-fi- the 

track the spin rate of . ; composed of a plurality of modes and frequencies, 

rt5ir Wii. by ...e army^ Alt that is re.,nired is 

It for any mode the ratio -/m be er,ual to the set delay rate , . 


(Subcase 2) 

To disenss the other cases where a«-I i. » ‘"li'e'' 

m par aid, . U^herr- be It, m!t. odes indicated by m • ai- P-mt in I m press,, re 
Mat this freduency, to . Then when ft has value Irld IK do ) 
follows; 


^sk'N m-M 


of) 


Ilfiicf m * fitl^Nwil! iilso in:iki‘ llu' iirniy Im lor / 1 , lliis iihmiis thiil il tiu'ii’ I'sisis >i 

liuiilr li;iviii)t , ils >.ln'ii)'|li will mniiilmh- liiily Id IIu' nnay siy.nal, even llu)in-li 

ils spin rate is laiMeally dilTeienl Irmn (he value ol tlie iarnel nnule "/(,(. In ilisenssions ot 
signal proeessint’, lliis type c»l eDniaininaliDu is ealh'il aliasing. 

To avoiil aliasiii)', a pmeednre analoj'Dns Id llial employed in ordinary sijtnal proeessinp 
is indieaU’d: At the liipliesl lrei|ni'ney ol inleresl, . lei m* he the laiy.esl mode 

ntnnber that can propap.afi'. II I he innnber ol deleelois is made larper than Zw . then 
lor any target M less than m* . |m|s|Mt.kMl will be larger than w*^and will corres- 
pond to a mode that cannot exist in the liucl. The reriniremenl N >Zr* may be recog- 
ni/.cd as the e<|uivalent ol the Nyipiisl crileiion in signal processing. In this way lalse 
signals due to aliasing are eliminated Irom the array output. 


(Subcase 3) 

In the following il is now assumed that the Nycpiist criterion. , has been satislied. 

When - I, an integer that is jio] a nudtiide of N, the sum of the N unit complex num- 
bers spaced apart by T2ir/n is easily seen to vanish. Thus the array signal vanishes for 
all inoiles having integer values (d'ytf‘«*-w^ ilillerenl Irom zero. Ihiis, il at a given frequency 
is set to track the stieeillc nuule>n‘f^ , t m , and/*<"*t^ . All wilues ol m 
other than make 4 a non-zero integer. Therefore, with N > Am . if a plurality 

of modes, m , exists at a Irciiuency co and tlie array is set to (rack a specific mode. 
M, the array signal will completely reject Ibe coniributions Irom all the other moiles. 

Case B : a non-integer 

When the array late ^ is set to target a stiecific mode-frequency combination, it has 
been seen that the array signal gives tlie amplitude of the target mode, rejecting com- 
pletely all other modes that may be present at tins freiiuency. It remains to detennine 
what happens to the array signal, if the delay rate is not precisely set equal to . 

and also how the array rcspoiuis to the presence of other modes with this mistuned delay 
rate. This information is obtained by evalualing the array lactor, /,. for noninteger values 
of . In particular, it may be convenient to sui)posc tliat is targeted to some 
Vco ratio so that /& is initially zero, and tt) then suppose that/^ increases as a conse- 
quence of a continuous change in the fre<|uency, W . of the target mode, rvt . This 
freriuency change, for example, could be proiluced by changing tan speed, aiui its ellect 
on the array output is obviously imprrrtant to determine. 

/.is evaluated for noninteger 4 by obseiwing in equation (.T.3-2‘f) llial NZ is the sum ot 
a geometric series of N com (ilex ternis, starling with I forn^o ami with subserpient 
terms having a ratio of exp i&z'lf/H. Using (he staiuhml expression for the sum of a 
geometric series, after some reariangenienl. Z can be shown (ref. 13) lo become 
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itic m:innitude oJ /, |Z| (sometimes called tlie array (lattern) is significant in determining 
the array response. In tiiis report |/)is called either the array response or the array res- 
pon.se limclion. It is not likely to Ik* confused witii the (generally) complex quantity 
/, which is called the array factor, hecaiisc of tlic notation used. Since the exponential 
laetor in erpialion (.^..^-40) is just a unitcotnplex number, the array response function is 
simply: 


2 - 


ii'ft /3it- 
N sift (SttAj 


wlicn the array is on target = Q and Z = I. Also it is dear that for /3 - 1,2,3... 

Z - (). When N is reasonably large (as will be the case when the Nyquist criterion for 
aliasing suppression is satisfied, the array response becomes c.ssentially independent of N 
and equation (3.3-43) simplifies to: 


\Z[^)\ 


I 

I /S'n- 




By means of Tigurc 3.3-4, tlie array response, (Z|, may be seen in a simple manner. If 
there is but a single mode at some frequency and if ^ is selected by ip » "/o to track this 
target, theii/^* o and Izj =1. As the frequency, to, is changed, the response falls grad- 
ually to zero when.tf = I and then beliaves as indicated. Fuirthcrmore, the effect of otlier 
modes in addition to the target may be seen. At the target frequency other modes given 
by<»< generate values of^= 1. 2, etc. and zero response. When the frequency of other 
modes differs Ironi the target the array response will be given by noting the behavior of 
\Z(j3 )l. 


To illustrate the array response as a function of the discrete variable m . and the con- 
tinuous frequency co. Figure 3.3-5 may be helpful. The array is illustrated when opera- 
ting with a specific ilclay rate^ so tluit the entire class of moiies for which are 

on target and the array responds fully. Mode-frequency combinations departing from the 
target ^ contribute responses as illustrated, by the shaded curves. 

With array .systems incorporating just a selectably constant delay rate, Ip = 9 = co/irt 

and incorporating no indepemlent frequency filtering, it is not possible to avoid signal 
contamination by mode-frequency combinations that ilepart from the target setting 

■ I his contamination is a result of the limited ability of such arrays to discrim- 
inate among waves having directions slightly off the target value. In linear array techno- 
logy (ref. 13) methods for shariJeiiing tlie array response have been developed. These 
include unequal detector spacings and unequal weighting of the individual detector sen- 
sitivities. It i-i reasonable to assume that corresponding methods may possibly be de- 
veloped to sharpen the res]ionse of the circular array. I he iileal response would be the 
S -luncfion; |Z ( ^ ))= S t/5 ) ^ (m - Wp ). For an array approaching 

such a response runction. the output for a discrete fretiueiicy field wouki be a set of 
^ -lunclions ol strength located at (reipiencies co - , and mode-1 lequeiiey 

comlunations departing from the target set would contribute trivially lo the array signal. 
.Similarly the random field signal psd would automatically Ik- the set ol inodal psd's 
Sjli or S« (iwlj given liy ei|uati(m (.F.f lO) and shown in I ij'ure 3. .1-7. 
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As far as using arrays to map llu- modal |>sil in ( ) spaa*. i1 is obvious that 

more information is rcipiiivd tinm can lie iibtaincd troin a circular array o| miciophoncs 
at a single radius in one transverse plane, It may be possible to extend the technic|ues that 
were explored in this section to provide the recinired information in general duct geomet- 
ries. Such extensions would involve the deployment ol array detector elements in either 
or both radial and longitudinal directions. To arrive at ellicient data act|uisition-processing 
systems to handle this problem would reciuire further investigation. 

3 . 3.4 Summary of Section 3.3 Results 

1 A simple mathematical model was obtained for the random acoustic pressure field in a 
thin annular duct, p«J,t), in the form of a sum of products of random functions of time 
alone with harmonic functions of space, fl, only. 

2. Starting with this model, an expression was produced for the modal power spectral 
density S (m, co) - the distribution over wavenumber - frequency spac’ of the riiean 
square modal pressure. This function is expressed in terms of the power spectral densi- 
ties of the random time series comprising tlie mathematical mode! of the pressure fiel . 

3 From the symmetry properties of tlie modal psd, S (m,w) was found to be resolvable 
into two unique sets of components, (|wl) and (Iwl). having mathemabcal 
properties similar to those of forward and reverse spinning modes in coherent, discrete- 
frequency fields. Because of the formal similarity, these components are called random 

field spinning modes. 

4 Means for experimentally determining the modal psd in terms of measurements m a 
transverse plane of a thin annular duct were devised. F.ithcr two separately movable 
(in 0 ) microphones or a set of fixed microphones may be used to acquire cross power 
spectral density functions which are combined and harmonically analy/.ed jspatia ly) 
to obtain S (m,co) and/or its spinning mode components, S (1^1) and b ^Ucui). ut 
same tecimiques may he used in coherent discrete frequency field analysis or for pres- 
sure fields containing a mixture of random and periodic components. 

5 A data acquisition and processing system was devised, using a circumlercntial array of 
microphones and incorporating a selectable, progressive time delay. The time delayed 
outputs from the microphones can be filtered at a specified frequency and added. Tht 
resulting signal from this system can then be used to determine the modal psd (random 
or coherent) in two-dimensional annular ducts. 

6, In general annular duets, having huh/tip ratios too small to be representable by the nar- 
row annulus simplification, the delayed circular array system will enhance all (circum 
fereiitial-radial) modes luiving a given frequency and a common spin rate, as it docs for 
the thill duct case. Although there is no longer a one-to-one corres|K)ndence hctvvecn 
mode spin rate and cutoff ratio, the common spin rate property ol tlie enlumccul mode 

set may have value in certain applications. 
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?, Tlic ordinary power spectral density of the delayed array signal is used to map tlie 
aeiuistic modal power spectral density function, vS and/or its components 
and S~(|c.>|). 

H. When the sound field is concentrated in integer harmonics of BPF, selective information 
about the total modal power associated with all modes spinning at a common speed is 
available from the delayed circular array signal without the need to map S (m.w) ,xten- 
sively. For a thin annular duct this selection on the basis of common spin rate is vquiva- 
lent to obtaining the power associated with all modes having a common cutoff r: tio. 

9. An especially attractive feature of the delayed array system is that it eliminates the need 
to execute any cross spectral data reduction procedures. Thus, the large number of rela- 
tively complicated cross-properties operations required by the basic (non-array) method 
for obtaining S (m,ui) are completely bypassed through use of the delayed circular array 
method. 

3.4 EXPERIMENTAL ASSESSMENT OF METHOD ACCURACY 

3.4.1 Objective 

To conduct an experimental assessment of the method to determine coherent fan sound mode 
structures. 

3.4.2 Approach 

An experimental program was conducted on a small-scale, low-speed fan rig at Pratt & 

Whitney Aircraft. The rig was capable of generating specified modes under controlled condi- 
tions. Six fan sound mode structures were generated in the fan inlet under conditions carefully 
controlled to maintain constant fan tip speed and ambient temperature. For each mode struc- 
ture, the MLP was run to dctemiinc suitable locations for the duct microphones. Because of 
the simple mode structures generated, all the microphone locations chosen by the MLP were 
on the outer-diameter duct wall. The microphone measurements were used to provide input 
to the MCP which in turn was used to calculate the amplitude and pha.se of each of the 
propagating duct modes for each of the six generated mode structures. The MCP also was 
used to calculate the resultant sound fields at other locations in the duct where microphone 
measurements were made that were not used to provide input to the MCP. Assessment of the 
method accuracy was based on a comparison of the predicted pressure amplitude and phase 
at these locations and the actual measured pressure amplitude and phase. 


3.4.3 Facility Descriptiun 


I he Pratt & Whitney Aircraft 25.4 cm ( 1 0 in. ) diameter re.searcli fan rig was used throughout 
the experimental portion of the program. This rig consists of a low axial flow, 2 bladed fan 
that has an operating speed range between 3000 and 6000 rpm. I he fan tip Mach number 
operating range is approximately O.1 2 to 0.24, and axial flows of up to 25 meters-per-secomi 
(SO ft/secl result. I he fan is cantilevered from a flaired flange .38 cm (15 in.) downstream of 
the fan and. in its baseline configuration has no inlet or exit guide vanes. 
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In nnk’r to comply with :i niiiNiimim iiumhci of the ;issiim|Mioi)s usal in rormiilalinp llic 
TylcrSolrin analysis, a constant annular inlet ilucl was constniclcil liy positioning a 
cylimirical ccntcrboily concentrically in (lie duel. The diameter ol’lhe centerhoiiy (or ID 
duel wall) was I 1.2 cm (4.4 in ), and it extcmlcd 2.H meters aliead of the rotor, with a sup- 
port strut 1 ,S meters ahead of the rotor. At this tlislance, the wake from (he support strut 
was considered nenlittihle. 

To reduce both steaily inllow distortion and unsteady inllow turbulence, a wire mesh screen 
was placeil about the inlet at a distance sufficient for there to be negligible pressure loss 
across the screen. 'I'he screen was supported by O.b cm (0.25 in.) diameter rods. The loss 
factor (k = Al’Al) the screen was calctdated from a formula due to Schubauer (ref. 1 5) 
to be 2.84 for the Reynolds number range: 50 200. At lower Reynolds number 

(typical of those through the screen when placed ahead of the fan rig), the pressure drop 
cannot be measured accurately, but indications from available information are that the loss 
factor increases at the lower Reynolds numbers. Kxperience from as yet unpublished work 
show.s that screens with loss factors on the order of two to three are effective in reducing 
inflow distortions; thus the screen material used for the tests conducted under this i>rograni 
was considered suitable. A schematic of the rig is shown in l-igure 3.4-1. and a photograph 
taken from ahead of the inlet is shown in Figure 3.4-2. 
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I lie reason tor iisiiii; an inlet screen to leiluce sleaily inllow ilistorlions was to eoin)>ly with 
tile assumption of uniform asial flow in tlie ryler-Sofrin analysis. Reilueins! steady inllow 
ilistorlions in turn reduces (lie amplitude of modes peneraled as a direct cause ol tlie steady 
inllow distortion (see ,\ppendi\ (' for such a veneration mechanism I. These modes (called 
extraneous modes in this text) are not predicted from a direct application of the Tylei-Sotrin 
analysis and, as wilt he seen later in section ,V4.t). cause inacciiiaLies in the inetliod presental 
for (he determination of coherent sound mode structures pretlicted (o exist hy the I yler- 
So frill theoiy, 

I lie reason lor usin'.' Ilu’ scieeu to reduce imsiead\ inflow luihuleiiee was to in(, rease test 
repeatahilitv and therein to ivdtua' the anionni ol lime needed loi sienal arerayiiii;. I Ids 
objective was ,u liie\ed el lei livels ,is <. ,ni In- seen in I leiiie t .T I Ills I i"ine siiow s the 
der iaiion (in ileciheh) ol the enhaiuckl sienal .iinpliluiK- ek-m-iakai In tlu- wakes liotii T1 


RESULTANT AMPLITUDE DEVIATION ~ decibels 


nuls inlcnulinj! wilh a hlailoil H)U)i al HIM- (oblaiiK'd liDni I'ivi' sainpK’ niilpiits al a 
j’.iwn mii’mplioiu') as a luiu lioii o) iho avi'iapiiij’ limo (or luimivr of samples). A comiiari- 
soii ol dll' lii'vialioiis is pirsonii'il lor ilk' (wo lonrij^iiralions with ami wilhoiil llu' ink'l 
scivi’ii. I( fan Ik' si'cii dial tlk- fonvi'rju'ik'o lor dk* fonfipmalion wilh tlu' inli'i sitivii is 
ohlaiik'il widi shorli'r avi’raniiip link’s ihaii lor dk’ fotdii'.tii'alion willionl dk' sl’I’l'i’ii. A 
similar ivsiili was also ohlaiik'il al HIM''. This improvi’iik’iK in ri'pi'alahilily aiul in lln' 
afi'iiracy ol moasitii'i! cohoront signals is athii’vi'il willionl sigiiififanl fliango in (In.' tom.' 
or hroailhaikl soiiiul pri'ssuiv lovcls. riio two I’igmvs. .1.4-4 ami .1.4-.S, show spcflra ohlaiiu’il 
rioiii a wall inoimteii niifrophonc wik'ii die intel coiilrol sereen was respectively renioveil ami 
iiislalleil. A eonipari.son ol' the spectra imlicales that die tone ami hroailhaikl levels are 
almost identical and, thus, the predominant sound field in the duct could not he generated 
hy rotor turhulence interaction. 
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3.4.4 Data Ai’«i»iisitii>u aiul I’rmc.vsing System 

An ■•(.n-lu.c- clal» pma'ssin,, sysl.,,, was .IcvcIoiksI I.. |..o.kK' c<.Ikt.s ,1 prossmv siBiuils fni 
thsciiioill IISL- hi cllh calLiihiliim ptiKc-iliirc to iIcK'niiine aiuplllii.lc anil pliasp ol 111.- u.ii 
lnp.,1 to llK- POK-C-Ssinp syslom is the- oyarall siphal l-om cad, micropho 
a one %vlll. a rolor lockc.l lipu-rclorciicc sieiial. Illc iwcrall acoustic pressure sienal "'t ' 
rc,,^^ nhutioi, irou, lwoco,„|.o„c„ts: a pcrioclic colicrchl sigual all.l a ram out inu lie 
L a The priitdpal output of the .lata analysis system is the eolictenl amphtmle am pl so 
m ea!-!, tuk-ropl,.'i.o locatiou, lit a.lclitiou, the ...itpiit iitelmles Itllercl tone ttotse a. .lesi ad 
rrcciuencics (BPl- and 2 X HIM’ for the purpose of the present test) versus fan rpm or tin c 
amt spectral analysis. Aiialop reeor.li,i6 capahilily was also provitlcl, so that ilata e..uW Itc 
Stored for any subseeptent analysis that may be required. 

F-ieh of these system components is discussed in detail along with an accompanying de- 
sidption or how they were employed .lurittP the experimental evaluation 
conmonents of the data analysis system include: 1 ) signal enhancement, 2) one tracking. 
3 ?rctml amilysis. and 4) analog recording. Pirst, however, a cleseriplion w, I lx presented 
of the instnimentatioi, for acquiring the overall pressure signal that is input to tile data 
analysis system. 

Instrumentation for Acquiring Input Signal 

For each conngiiration tested, ten microphone systems were useil to measure the sound 
Held at ten locations on the duct wall. These ten locations incorporated 

phones and locations as specilied by the MLP for each contiguration. and he ruiuinder were 
iLd to provide data for checking the method accuracy. Locations lor each conliguralion 
were detennined prior to the test and provision for the different locations required were made 
bv drilling 14 boles in the test section of the duct, rhese locations are shown for re crcucc 
h '^ure 3.4-6. The microphones used were B&K, model 
condenser type They were calibrated for output sensitivity bclore and al ter each test r in 
with a B&K model 4220. pistonphono calibrator. The microphone and pistouphone calibrator 
were certified acceptable by the Pratt & Whitney Aircraft standards laboratory with the 
pistouphone output being certified at 124.0 * 0.2 dB. All microphone outputs were normal- 

i/ed at 1 .0 volts equal to 1 24 dB. 

Data Analysis System 

■flu- iuput lu 111.- .1.1111 ..Ulllysis xyslV" »» “'”Vr'’"'rU,?u-Thr“' 

liipc U. lx- pr,K«w.l lit .1 liilCf llmv Hy Hu' •'.ui lilK-- iltslr,u„c,il;,ll,.u. I lu- . ii lu, ,1.U. 

.millysls »ysl.-,„ uiw.l lluouglmul 111.- l«l proprum iv »l'>'w" 

tint in addition to providing signal enhanced amplitude and phase values, the on-lm 
"giliil all L s,k1i,„ iiimly/cl. 11II.T.-.I ill .l.-vi.c.l flcpi.-iuio. iill.t pl.itt.-.l wrsus tin, 

or time. 

In l ieure 3.4-b is a block diagram of the magnetic taiH- recouliug system. I his later system 
was developed to document the acoustic piessiire signal at each microphone location as 
stipulated by the experiiuental assessment luogiaiii. 
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Sij'Mill I iiliaiiivini'iil 


I o pmviik* i-oliiTi’iil acoiislii- prossutv sipiial . ai mii rnpliom’ loi-alioiiN iik'iililu'il hy ilu‘ 
MU’, a inloi-lockrii (imo ivI'iti'ikv was olXaiiii'il I'rnm a spci'd lri)T.i*i iiip sysiom (hat iisi-.l a 
uMUl'i fiuiiil to pi'iK-ralt' a piilsi' I’vt’iy tinu' a inpjr kkuk' passi'tl a pnim of ivlViriui', TIu’ 
liipfiiTiiip sysU-m opaali'tl on llu' Pit ivvoliilion spi'i't! sipnal amt m’lUTaloil a piilsr 

at Ilk' A’ld irossitij' of llu* sifjiial. This /.mo crossing occmroii at ftpial linu' iikiviiu'tils as 
mu-h hhuk' passoil. Tlk.* otiipiu was a transistor-transistor lopic Cri l ,) pulse that was time 
loekeil to the rotating! Ian, I he I IT pulse was used in sinnal enhaiieenient to deterniitu' the 
zero referenee time and repetition rate of the enhancement. 

nuriii}; the test, all microphone signals were nitered at cither Hl’l- or 2 X Ul’l' and stored hy 
the memoiy ol the signal enhancer. I his storing was achieved hy digiii/ing the analog 
input over a time period sulficient to aec|uire at least one lull cycle of the input Iretiuency. 
Suhsequcnt cyeles, triggered hy the ITL pulse, were digitized and added to the memory. 

A total or.12.7(>8 summations were achieved lor signals at Bl’l', and si.immations, 

lor signals at 2 X HPI-. The summations were chosen as optimum based on lepeatahility 
studies. 

The enhancer output was dte average of the smmnations. and thus the coherent portion of 
the microphone signals was obtained, in order to measure signals representative of more than 
one blade pressure pattern, the blade related ITL tiigger signal was fed to a preset counter 
that generated a puke on every third blade. The reXor assembly contained 32 blades, so 
triggering on every third blade caused the enhancer to sample all of the rotor-locked pres- 
•sure patterns every three revolutions of tlie fan. The enhancer output, representiim the 
coherent portion ol the microphone signal, was plotted on an .v-y plotter from which tiie 
averaged signal amplitude and phase could be obtained. 

Tone rracking 

1 he overall tone sound prc.ssure level is also an output of the data-anaivsis system i his out- 
put IS used to demonstrate the cut-on of a specillc mode by increasing the rotor speed and 
recording the overall microphone signal for the blade-passing tone or harmonic of interest. 

le tone level at a specitic trequeiicy coiiki also be compared with the coherent amplitude 
to determine whether most of the sound energy is either coherent or ineoherenl. 

si.uial loi signal enliancemeiit. a tracking ratio tuner and a tracked filter system was rc<iuiie*l 

I Ins capability was provitled as follows: 

The microphone signals were filtered with a 50 hcil/ phase-locked bandwidth niter centered 
the irequency ol interesttm’l- 2X m>lT. The tracking filter was continually tuned by 
the .t pulse per revolution speed signal which was passed thnnigh a phase coherent liackin.. 
i.ilio iiinei that had the capability ol lollowing any nxoi speed variation and to tune the track 
mg tiller to any multiple of the blade passing freiiuency. I he ,?2 pulse per revolution siun il* 
was input to the tracking ratio tuner and was inulliplied by one or two to prov ide the proper 
^luviiig licpicncv to the trackiue lilter f<n- HIT or 2 \ HIT', respeclivelv . lor plots of lone 
iio.se vers.K rpiii or time, the output from the tracking filler w,,s dirccled lo the V ,,xis of 
the \ y pl‘’iki I he "Dt inoportioual to lic<|uency" output Irom the tracking ratio l.mer 
piovulcd the .\ a\is (rpni) signal lot plolliiig diiiint' accel deeds 
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Spivtr:il Analysis 

I o assi’ss ilu' spi'i lral iiniU'iil of tiu' ovoralt miiroiilioiu' signal, a spi-ciial analy/i'r was in- 
foipora It'll in Ilu* ilata analysis systfin. Spectral plots coiiU! he ivailily assessed to deleniiine 
whetlier the sound pressure level of ilie hlade-passing lone or its harmonies ilominate the 
hroailhand noise component. Speciral ilala obtained during the lest program indicate that 
this was always (he case. 

A tnnallel output lYom the mieroplione was connected to the high speed speemun analyzer 
winch allowed simultaneous spectrum analysis during signal enhancement, riie analysis was 
made over a rrequency range I roin 25 hertz to I OK hertz with a 25 hertz hamlwiiiili lilter. 

I he analyzer provuled tiigital storage and averaging of the analyzed signal. In order to obtain 
data sampling consistent with the enhanced data, the analyzer was set to accnmnlate .S 1 2 
.spectrum averages. Stored averaged spectra were then plotted on the x-y plotter. 

Atialog Recording 

An analog tape recorder was employed to record the microphone output to enable analysis 
similar to that pertormed on-line. I he advantage of llte recorded analysis was tliat each 
microphone signal could be analyzed over the same time period; however, ilisadvantages in- 
cluded the careful calibrations required to presers-e phase relationships. 

Data trom the ten micropiiones were recorded on magnetic tape simultaneously for eacli 
conllguration tested. I'ive micropi.ones were recorded on tlie odd record head; and five, on 
tlic even head. The .^2 ptdse per revolutiot) speed signal was recorded simultaneously on the 
odd and even heads. Ml data were recorded on the wide band (ironp 1 recorder in the I'M 
mode at a tape speed of .fO inches per second, l igure .^.4-8 is a sketch of the magnetic tai)e 
recording system. 

( .ilibration signals were recordeil on all the recorder channels simultaneously, incluiling the 
speed channels. Signal enhancement of these calibration signals provided a measure of in- 
dividual data channel amplitude sensitivity. The enhancement also provided a channel-lo- 
channel and head-to-head relative phase correction. Data recorded on the mid head were 
enhanced using the speeil signal recorded also on the odd head and similarly for the even head 
data. 

3.4.5 Test Program 

I he matrix inversion melhoil tor determining coherent fan-.sound-mode structures was assessed 
for six mode structures. The assessment was made with measiireinent input obtained from the 
test program described in this section. Ihese modes structures were generated In' two sets of 
distortion generating rods on the Pratt Whitney Aircraft 25.4 cm 1 10 in.) diaiiielei research 
Ian rig. I our o( the six mode structures were generated by one set of hardware which 
simulated a stator wake-rotor interaction by positioning .14 equally spaced rods of O.Tt cin 
lO.dO in.) iti diameter. ,T7 I cin ( l.4(> in.) trom leading edge to trailing eih^e upstream of the 
rotor, for this conliguralion. the bladed rotor was operated at both .t lOO ipm ,uui sN 1 5 




iptn. Mii rophoiu* l(K■;llin^^ wi'iv ik'ti'i inltu'il for imc!i oI IIh'm' imuhil simi iiiivs hy iiso nf ilu' 
MKropljoiu' l.(K:iti()ii I'loj-nnn ( MU*), as disnissi'il in sri'ioii d.4.4. Al i-adi mini' s|uvd. 
liala \u-iv nhlain al kolli Ma«k‘ passing riv»|itoiuy (HI’I ) ami iwii i- hladc passing Iiv(1ik-iu y 
(2 X Ul’l ). In lliis way. lour modal slHK liiri.-s wciv calt iilali-d hy llu- Modal (’aU ulalion l*ro- 
y.iaiii (M('l’). 

1 Ik' <MIu'I sol ol haidwaro siimdatoil steady inllow dislorlion hy plaoinp a siiiftlo rod assoinhly 
in IrcJiil ol rlio lador, I wo modal siniotiircs wore evaluaiod hy oporaliny iho rotor al two 
speeds (.MOO rpin ami .M40 rpm) and aeiiuiriiui data at HI'I'. l-or this emiUjimatioii. the MM* 
was used to loeate mierophoiios lor the modal slriieture generaled at the .^400 rpm rolor 
speed, l-or the higher speed, which represents a ten perceni increase in rotor speeri. the same 
set ol microphone positions were u,scil as tor the lower speed. The objective erf this test was 
to determine d a set of microphone iivaiions as determined by the MLP for a given rotor 

speed and lret|ueiiey. can be used successrully to obtain data at a slightly dilTercnl Iretiuency 
and rotor speed. 

A summary of the lest program is provided in fable 3.4-1. This table includes: the geometric 
eonliguration, rotor speed, frequency, the total number of propagating incident nunies 
supported by the duct, (lie mode structures predicted by the Tyler-Sofrin analysis to be de- 
termined for each ca.se. and the conditioning number associated with the microphone locations 
as determined by the M LP for each case. 

For each of the six mode structures, the amplitude and phase of each of the modes was calcu- 
lated using the M( P. Input to the M( P was obtained “ondine'' from a set of microphones at 
duct locations rletemiined by running the MLP. as described in section 3. 2.3. .Additional 
acoustic pressure tneasuremcnls were obtaineil al other localions in Ihe iluct. and were used 
to eheck the method accuracy by comparing the measured sound field at various locations 
with the predicted sound field at those localions baser! on a knowledge of the calculalerl mode 
stiucture. Shown in lable ,v4-2. lor each lest configuration, are the microiihone locations 
used to provide data lor the modal calculation and the localir>ns for Ihe microphones itseil to 
check Ihe accuracy ol the moilal calctdalion nieilunl. 


In addition to the.se measurements, an evaluation of Ihe slaiidard deWalion for the errors in 
measurement of acriustic pressure aiul micro|ilu>ne location was eoiuiucleil for each of the 
six mode struc,lures. lable 3.4-3 presents, for each lest configuration, the values used as in- 
put to (he M( P so that the standard ilevialnnis for each of the ealculaled modes conUt he 
determined. 

I he standard deviations rif the errors in microphone location (as obtained from Ihe Pi'iWA 
.Standards laboratory) are relaled lo (he tolerance of the measnremeni lor delermming ihe 
axial and azimulhal c<iordinate ol (he microphone. I runs in the radial localion of flush 
mounted (ransdneers retied Ihe variation in the ronndness of die outer duel wall and ihiis 
is seen lo be small. 
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TABLE 3.4-1 

TEST PROGRAM SUMMARY 
32-BLADE ROTOR 
(25.4 CM FAN RIG) 


MODE GROUP 


Test 

Designation 

Stator 

Speed (rpm) 

Frequonev 

Duct Supported*®^ 

Mode Structure to^*^^ 
be Calculated 

Conditioning 

Number 

Stator-rotor 

interaction 

34 rod 






1 


3400 

BPF 

7 

(-3,0)(-2,0)(b)(.],o)(0.0) 
(3.0) (2,0) (1.0) 

2.1 

2 


3400 

2XBPF 

20 

{■4.0) 

- 

3 


5813 

BPF 

16 

(-2.0) (-2,1) 

1.3 

4 


5813 

2XBPF 

53 

(4.0) (-4,1) (-4,2) 

2.0 

Steady inflow 
distortion 

1 rod 






5 


3400 

BPS 

7 

(-3.0) (-2,0) (-1,0) (0.0) 
(3.0) (2,0) (1,0) 

2.1 

6 


3740 

BPF 

7 

(-3.0) (-2.0) (-1,0) (0.0) 
(3.0) (2.0) (1,0) 

3.1 


Notes: (a) Not including reflected modes. 

(b) (-2,0) is only one in the group 

predicted by Tyler-Sofrin Theory 
for a stator-rotor interaction 


(c) Predicted by Tyler-Sofrin analysis 
for a stator-rotor interaction 
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he standard deviation of tlie inaccuracies in acoustic pressure measurement are not only 
due to measurement system variations, hut also to the different sound fields at the different 
requencics considered during tlic test program. The accuracy of the microphones to measure 
acoustic pressure amplitudes was assessed by comiiariiig a precalihration and a iwstcalihra- 
tion ol tlie microphones for each test. During this lest program it was experimentally man- 
ageable to adhere to carelul procedures .so that the standard deviation of the errors associa- 
ted with transducer calibration was relatively small (about 0.2 to 0..1 dm. In addition to 
calibration errors, the standard error for pressure measurement includes an evaluation of 
measurement repeatability for each of the six mode structures. The variation in repeatability 
was as.sessed m a similar manner as the deviation shown in l•■igure 3.4-3. By this procedure 
an enhanced signal was obtained from five sample outputs at a given microphone at the num- 
ber of suminalions chosen as optimum based on such repeatability studies. The standard de- 
viations ol the variation in repeatability are typically larger than those attributed to trans- 
ducer response. I Inis, a combined standard ileviation (in decibels) of both calibration errors 
and variations in repeatability, as presented in fable 3.4-3. yields values that are controlled 
by the relatively larger observed variations in the repeatability of the enhanced signal. 

In light of this observation, it is expected that the magnitude of the variation in the enhanced 
.signal will depend on whether the random contribulion in the overall signal is comparable to 
Ihe coherent sound pre.ssiire level, fhe ratulom component of the overall sound pressure level 
was compared to (he coherent level for (he two mode structures that were generaleil at Bl'f 
(test ciniliginalion I ) and 2 Bl’f (test configuration 2) by the wakes from 34 rods interact- 
ing with a 3_ bladed rotor operating at a speed of 3400 rpin. 

*o'V I‘>I live .sample outputs at a given microphone 

was 0.3 dB. However, it was found that the cohereni sound pressure level at five microphone 
locations was tio.n T | p, | l.s dB higher than Ihe random component, flius. in this case the 
randomly phased noise structure was on an order of 10 dB below the coherent sign il so ihc 
variation in amplitude and phase of the enhanced signal is assigned Ihe smallest pressure mea- 
surement error of those given in fable 3-4-3, 


A( : :i Minilar ;m:itysis :il lliicv iiiiL ioplioiu’ l.u iiiioiis w;is lomim tod iiy vompni iiir. 

sifinal I'nhaiuvil ami ovvrall lorn- lovol data, l ids analysis iiulicali's tlial llu’ cnhoivtil noise 
iK'ld is iK’lwocn -I.: lo 5.-I ilU liiplicr Ilian (lie lamUmdy pliasml sirnclniv. 1 lie pressiiiv nu-a 
smvineiu error 010.5 dH in 1 able .^-4-} lor test eonfinmalion : relleels the larger variation 
in pressure anipliUide measnivnieiu due to the relatively liielier randomly phased eomponeiit. 


A llnal eoinparisoii ol the overall lone level and siunal enlianeed i>iessnie amplituile Uiis eon- 
iluelcd for tlie tlfth mode slnietmv eonfijuiralioii. In lids test eonfiimration a sinjiie roii was 
placed upstream of a M hladed rotor o|ieratini'. at 3400 rinii. At several mieropiione loea- 
tions the measured resultant pressure amplitude was small due to the eaneellation ol the eon- 
stilueiU modes, l lius. a eoinparisoii of the relative ami>lilude of the coherent and randomly 
phased sound fields varied .substantially from one microphone location to another. In parti- 
cular. at microphone location 5 the coherent souiul field was lO.b ilB liiiihcr than the am- 
plitude of the randomly-iihased component. However, at microphone .3 the measurcil en- 
hanced pressure amplitude was 1 b.3 dB lower than the sound pressure level for the random 
structure. .As may be especled. accurate measurement of the enhanced signal at some micro- 
phone locations was difficult so that in fable .3,4- .3 a relatively larger stamlard deviation was 
assigned to the pressure amplitude error. 


In the next seetion, .3.4.0, the amplitudes and phases and the calculated standard deviations 
are presented for each of the six modal structures investigated. In addition, the accuracy ol 
the calculated mode structure is demonstrated by comparing (be measured sound lieUl at 
various duct locations with Hie |uedicled sound field at those locations based on the calculated 
iiunle strudiire. Discrciiancies between prolicled and measiited ilala are discussed and related 
lo the calculaleii standard deviations associated with each ol the modes. A sutnmaiy atui .m 
assessment of the lest results arc presented in section .3.4. 1 . Keailers wishing lo avoid the 
detailed test description may proceetl direclK to section .3.-t.7. 

.3.4.6 Test Re.sults 

Isesults from each of the six inode structures will he presented under a separate heading. An 
assessment of the overall test results will be leseived until the next section. .3.4.7. 

Mode Slniclurc Configunilion No. I 

I'he iirsl modal structure was generated .it Bl’l ( l-Sl.t II/) b\ the wakes from the .34 rods 
interacting with the .3.’ blailed rotor operating at a speed ot .3400 rpiii. fhe overall acoustic 
spectrum generated al this speeil at microphone No. 10 is shown in l igure .3.1-0, and it c;iii be 
seen that BIM- is well alnwe the background broailbaiul noise level as rletermmcd Iroin a .V*' 11/ 
bamlwiilth filter. Al Bl’f. the fan inlet can support seven mciileiil imules I ,3. 0). ( .\ ()), 

(-1 , 0), to. 0). 1 1 , 0). ( 4. 0). and t .3. 0). Based on an .ipplicalion of the I yier Solnn analysis, 
it was expected that the sound pressure le\el of the ( ’. Ol woiiKI I'c much greater than the 
sound pressure levels generated by the lemailiing six incident tm.des. However. I'ccaiise ol 
the low number of imules supported In the iliicl. it was decided to locale seven microphones 
in the rind to provide data for c.dculaling the aniphtudes .nul I'b.ises ol all seven incident 
modes. It w.is assumed in this lest that modes letlec'.ed Iroiii the inlet do not conliibule 
significantly lo the sound field I he aclu.d Uv.ilions .is delei niined b\ the Ml T ot the seven 


iiiiciopIioiK's lire f'ivon in I ;ililc ;iiut ii‘sulk'il in n loiidjlioiiiii)! mnnhi'r (see v.‘i iion 
4. ) nl 1 . Also, ihidt wi-n- ohi:iiiu-(l :i( llnvi- ollu-i niii rojiiioiU' lornfions so IhnI ;i i lifrk 

i>ii Ilk' ;it\'ur:icy of llu’ i.';ili‘iilaloil imuk’ slniidni' would ho ohlainod. 



I 'i^ure .i.4-V On-rall Si>crlniiii ( n-nauicil ui J-ioo ,ii Mia\ i/i/ioiir l.in iition \'t>. 1 0 ,‘iUiilc 
Sinicdin' ('oiitiyiiriiiioit .Vc. / 

'I'lu.’ ooliL’R’nt acoustic pressme amplikulc and phase ,il all ten micio|)hoiic locations arc listed 
in I'ahle 3.4-4. The Hist seven aiiipi.uides and phases were used as input to the M(T. and file 
atnpliUules and phases of the seven p .ipa.ealine incident inodes were calculateil. In addition, 
the slandanl deviations lor the anipliiude atui ]>has.' ol each mode were calciilal d iisini! the 
standard rievialions due to anipli lude ami phase iiieasiireinent errors, I hese slamiard devia- 
tions in measurenienl errors weie deliueil in l ahlc 3.4-.t. The calculateti amplitudes and 
phases ol each moile aloiiii with the corresi>ondiiiu siamlanl deviations are listeii in I'ahle 
3.4-5(a). rite n odal ami>liludes and error hands niven hy plus and minus one slandanl devia- 
tion are illustrated in I'ifiure ,h4 It). 

The imuial structure was used to predict liie ivsiiti.ml sound Held at three iniei'ophone loca- 
tions: No. X. No. d, and No. 10 oT l ahie .H I .'!. 1 he ilisctepancy helweeii predicted and 
measured acoustic inessuie ainplilude (m dicihel.) atui plane im decrees) is listeii iii lahlc 
,V4-f). It can he sei'ii that thi‘ ph:i.-.e is predick’d sliehtly hiyh (3 d pi ii. ' "i at .ill Ihrei' loca- 
tions and that the attiplilude is u i ihin O..S dlt to I ' dH Sin.e three check lo. at ions are a 
statistically small population, it is dillicii!l lo assess .icciu.ilch the method loi this lest case, 
l o net a hetler estiiiiaie. advanl.iyc will he taken ot the l.iel that the ( ()( inoile is iinieh 
p.reatei in amplitude than liie olhet piopajMliii!' modes If die . onti ihiilion liotn modes 
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lo delerniine the amt'lilude and ph.ise o| the dominant mode I he ten mu tophone niea aiie 
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pliasi' ol the ( 3. I) i modi.- e.ni he i lei ei mined .ilonr vs it h i he te-.pect i\ e sl.mdai il de\ la lions. 
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ccluation system so as to yield a conditioning number of 1.1, Data from the three .selected 
microphones were used as input to the MCP and tlie resulting mode amplitudes and phases 
along with respective standard deviations were calculated. These values arc shown in Table 
3.4-S(b) for the (-3, 0), (-2, 0). and (-1 , 0) modes. Comparison with the mode structure ob- 
tained using seven microphones shows a discrepancy in amplitude of 0.9 dB for the dominant 
(-2, 0) mode and on the order of 1.5 dB for the other two modes. Also the calculated phase 
of the dominant (-2, 0) mode changes from 243® to 246® - a value further removed from the 
estimate of 241 .8® based on the average of ten microphones. 

Thus it can be seen how the omission of only the four smallest modes out of the seven incident 
propagating modes affects the predicted amplitude of the dominant mode. In this case the 
differences are small, showing that a good estimate of the dominant mode amplitude and 
phase can be obtained without considering all seven modes. Despite the small difference, 
however, it is interesting to note that the omission of the four smallest modes (more than 
20 dB down from the (-2, 0) mode) results in a 0.9 dB change in the calculated amplitude of 
the dominant (-2, 0) mode. 

Mode Structure Configuration No. 2 

The second modal structure considered was that generated at 2 X BPF (3627 Hz) by the wakes 
from the 34 rods interacting with the 32 bladed rotor operating at a speed of 3400 rpm. At 
this frequency, the fan inlet can support twenty incident modes although only the (-4, 0) 
mode is predicted by the Tyler-Sofrin analysis to contain significant acoustic energy. For 
this case, it was beyond the scope of the program to make a sufficient number of measure- 
ments to determine the amplitude and phase of all twenty incident modes. It was as.sumed, 
therefore, that the incident (-4, 0) mode contained the bulk of the acoustic energy in the 
duct sound field. With this assumption, only one microphone measurement was necessary 
to determine the amplitude and phase of the (-4, 0) mode and the respective standard devia- 
tions based on a multiplication of the influence coefficients and estimated mca.surement 
errors. Although the choice is arbitrary, microphone location No. 2, listed in Table 3.4-2 was 
selected to provide input to the MCP -principally, of course, to determine the standard 
deviations of the mode amplitude and pha.se. In order to have a large .sample of check micro- 
phones. measurements were made at ten microphone locations. The measured amplitudes 
and phases at each microphone location arc listed in Table 3.4-7. 

The amplitude and phase of the (-4, 0) mode and the corresponding standard deviation as 
calculated from the influence coefficients and measurements from microphone No. 2 arc 
shown in Table 3.4-8, Also shown is the amplitude and phase of the (-4. 0) mode along with 
the standard deviation in amplitude and phase obtained by averaging the ten microphone 
measurements. The main observation to be seen from this example is that, as was seen in Hie 
first test configuration, the predicted errors based on (he inDuence functions, underestimates 
the actual measured variations in amplitude and phase. A comparison of the average mode 
amplitude and phase with that obtained from microphone No, 2 has no real significance in 
this example, because of the arbitrariness in the initial selection. 
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TABLK 3.4-7 

MKASURBD DUCT ACOUSTIC PRliSSURB 
(RPM = 3400 2 X BPF = 3627 Hz) 

CONFIGURATION NO. 2 


Microphone 

No. 

(dB) 

Coherent Resultant Amplitude ^ 

(dyncs/Cni^) (lhf/in.“) 

Phase 

(deg.) 

1 

108.4 

52.89 

7.672 X lO*"* 

233.1 

2 

106.3 

41.47 

6.015 X lO''^ 

330.9 

3 

106.9 

44.35 

6.4.33 X lO'^^ 

320.6 

4 

107.8 

49.24 

7.142 X lO*'^ 

285.1 

5 

109.8 

61.59 

8.933 X-LO’^ 

38.5 

6 

106.6 

42.88 

6.219 X 10"^ 

325.8 

7 

105.5 

37.85 

5.489 X 10"^ 

133.2 

8 

108.2 

51.20 

7.425 X 10‘‘* 

145.5 

9 

107.4 

46.89 

6.801 X lO’'^ 

195.6 

10 

106.9 

44.09 

6.394 X 10*'^ 

64.3 

Avg. Resultant: 
Amplitude B^ = 

107.5 

47.24 

6.394 X 10'“^ 



TABLF 3.4-8 

CALCULATED {-4.0) MODE AMPLITUDE AND PHASE 
AND CORRESPONDING ERRORS 
(ONE MODE INPUT TO MCP) 

CONFIGURATION NO. 2 



Amplil.idc 

Error 

Phase 

Error 

Mode 

(dB) 

(dB) 

(deg) 

(deg) 

(-4.0) 

1 14.2 

0.6 

317.3 

2.1) 

Avg. (-4.0) 

I'rom 10 iiticrophoiies 

115.3 

1.2 

315.5 

10.4 
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MODE AMPLITUDE DECIBELS 



An ailililioiul (.•slimaU* ol' llu' nuMhoil ai i iirai j was ohlainoJ 1oi tins niiu' nuHle laso, 
ivpoalin.u till' U'sl wilti tinvo aiMiiiotial miiroplioiu' lo. ions ik'sifinali’il No. I 1. No. I’, aiul 
No. 13 ill I'ijuiro 3.4-P. l lu'so Ihivo mkioplioiu's loiu'llK’r with inii rophoiu’s No. 7 ami 10 
(iiOi usal to proiiOo input to du' l ak iilaiion) pioviilml data, so that tlu' iiK’asmvil sound 
fu'ki at tiu'sc s]K'v irii- liiK t loialioiis I'ould Iv (.•onipaivd with tlio pivditMod sound luMd hasi-il 
on ilk- (.ak'ulatod inodo stiik'tmv of niik’ monos. TIu' ditfoivna’ Ivtwoon pivdiolod and 
inoasiiivil ampliuuk’ and phaso aiv shown in l ahio .1.4 10. It can ho soon that tho piodiotions 
aro jionor.illy juuul. I ho aiuplitudo is proiliotod somowhat low by about 0..S dH to I .d dll 
with a standard do\ iation of 1.4 dll about tlio moan and tho phaso is within vO" am! -S" with 
a standard doviaiion of S.(f about tho moan. It is diffionlt to oslimato tho aoouraoy with 
whioh tho modos othor than tho (-4. 0) aro prodiotod. Howovor. it is soon from I'igiiro .1.4-1 1 
that Iho staiulard iloviation of modal ami4iludo oalonlatoii from tho inlluonoo fnnotion. is 
lariior for modos with low sound prossiiro lovols, 

Modo .Struoturo C’onfigiiration No. .1 

riio third modal struoturo oonsidorod was that sionoratod at Hl’i- i.lltH) II/) by tho wakos 
from tho .14 rods intoraoting with tho rotor oporaling at a .spood of 581.1 ii'tn. .At this 
Iroiiuoncy. Iho fan inlot duot oaii support siMoon luopagaiing inoidoni modos. I rom an appli- 
oation ol tho 1 ylor-.Solrin analysis, tho (-2. 0), ami (-J. I ) modos woro prodiotoil to oontain 
tho bulk of tho aotmsiio sound onorgv in the duot. I'oi this oaso. tho MI T was usod to solool 
miorophonos No. 7 and No. S Irom Ihoso lislod in lablo .1.4-7. and tho miorophono oombiiia- 
tion luui a oonditiouing luimboi of 1..1. .Again, inoasuromonls woro mado at a total of ton 
niiorophono looativms tmioroidiono No.'s 1 to 10 in I ablo .1.4-.?) to proxido oight looalions 
to ohook Iho aoouraoy ol tho oatoulalod modo struoturo I ho moasurod amplitudos and phasos 
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Ml ASUKi n DUCT Acousnc I’Kl-SSUUT 
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CONI ICURATION NO. 


Mierni'lume 


Colu'ieiit Kosiiliani 

Ainplitiuie 

Pliase 

No, 

fdlD 

(ilyiu's enr ) 

(Ibf'iii.') 

tde«.) 

1 

l.tl.2 

725.4 

10.521 \ i0‘-^ 

5.5.2 


12.S..? 

521.0 

7.5(,5 X I0 -’ 

.528.0 

> 

,y 

12S.*) 

555.') 

8.0(>.5 X lO'-^ 

(>1.0 

4 


441.2 

O..508 X tO -^ 

205. .5 

5 

1 ,U).8 

o‘)2.() 

10.0.5O X 10'-' 

288.0 

h 

125.5 

.577.2 

5.471 X 10'-^ 

144,4 

1 

12(>,l 

405 .8 

5.8S(> X 10-^ 

201 .8 

S 

i.b)..? 

(>54.1 

0.487 X lO'-^ 

107.8 

n 

125,7 

,58(>.2 

5.001 X UV' 

.5.5(>.0 

10 

12(vS 

42>8.0 

0..55.5 X 10 -' 

22(..l 


TIk’ •mipliliiitt.'s :iiui plinse.s of llio (-2, 0) ami (-2. 1 ) iiunlos ami iho I'ospm'tiw 

slaiulani iloviations laU-iilalod from llic inlluotK-o c-X’lUcicnls ami osliniaU'il moasiiri'iiu’iil 
(.•rrors tor this t.asL’ aio shown in Table .T4-I I'Ik’ preiliclal sonml lu’Ui based on the 
ampliliides and phases of the t-J, 0) and (-2. I I modes was eompaied with nieasnred data at 
the eijilu cheek mie.ophone locations and is piesented in Table T.4-I.V It can be seen Ironi 
this sample that the pi-edicted ampliliide is within - 1 .4 dH to 4.(1 dll witli a standard deviation 
of 1.5 dlT and the predicted phase between -24. U' and I 1 .4' . with a standard deviation ot 
14.5T 

In an elTort to determine the amplitudes and phases ol' modes other than the ( 2. 01 and 
(-2, 1 ) modes and. thus, improve the predicted accuracy at a set of check microphones, a 
lariter set ol' microphones and modes was selected to test the methoil in a manner similar 
(o that Used in the last test case, .\eain tiitie modes were selei ted so as to be eronped about 
the predicted dominant t-2. 0) and t-2. I ) tnodes. The nuules selected were: (-4, Ol. t-.T 0). 

(-2. 01, (1. 0). lO. 0). 1-2. 1). (-2. It rellecled. (-1 . 1 ) and (0, I). The renecled t-2. I I mode 
was incimleil in the selection because the inciiient (-2. I ) mode had a culotT ratio neat one 
ami a sienineanl portion of the acoustic enerey was espeeted to be rellecled. 

I he nine modal amplitudes ;md phases and the respeetiu’ standard devialit'iis are shown m 
Table .T4-I4. The amirtiludes ami ami'litude standard ile\iatioiis are shown nrat'hic.ilK in 
T inure .1.4- 1 2 It can be seen that Ihe l-.T 0 ) mode and Ihe I - 1 . I ) modes are quite comparable 
111 am|dilude to Ihe t-2. 0) and ( 2. I ) nu'des illusli.ilinn tor this v .ise the weakness ol the 
.issumplion that the sound field Is eonlrolled onK In the modes piedu led In the l\lei 
Soil in aiialv sis_ 
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A)',aiii. this lest was rirpcatiHl with check microphones at four locations different from the 
oriidnal ten. A coni|iarison was matle of the measured and predictcti souiul field based on the 
calculated amplitudes and pliases of the nine modes at these four microphone locations and 
a fifth from the orininal test that was not used to provide data input to the M(T. This com- 
parison is sliown in Table 3.4-15 along with the standard deviations in amplitude and phase 
of the resultant sound field at the five microphone locations. It can be seen that the accuracy 
obtained using the calculated structure of nine modes is similar to that obtained using the 
calculated modal structure from two modc.s. It is possible that no improvement was obtained 
because one or more modes were omitted from the modal calculation that actually contribute 
significantly to the duct sound field. 

Mode Structure Configuration No. 4 

The fourth modal structure considered was that generated at 2 X BPF (6200 Hz) by the wakes 
from the 34 rods interacting with the rotor operating at a speed of 5813 rpm. At this fre- 
quency the fan inlet duct can support 53 propagating incident modes. From an application 
of the Tylcr-Sofrin analysis, the (-4, 0). (-4, 1 ). and (-4, 2) modes are predicted to contain 
the bulk of the acoustic energy in the duct. For this case the MLP was used to select micro- 
phones No. 3, No. 9, and No. 10 of those listed in Table 3.4-2. Measurements were made at a 
total of ten microphones locations (as in previous tests), thus providing seven locations to 
ctieck the accuracy of the calculated mode structure. The measured amplitudes and phases 
from these ten microphones are listed in Table 3.4-16. 


TABLli .1.4-15 

DIFFHRtNt'E BETWEEN PREDICTED AND 
MEASURED DU(T ACOUSTIC PRESSURE 
(PREDICTED MINUS MEASURED) 

CONFIGURATION NO. 3 


M icmplutiic 
No. 

Ampliuidc Della 
(dB) 

Phase Della 
(deg) 

.3 

2.4 

* 

17..1 

1 1 

• 1.7 

1.7 

12 

1.4 

1 1.8 

1.1 

0.2 

-2.9 

14 

0.1 

-.10.2 
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■l lu' caln.latc.! ampliuuio aiul pliaso of iIk- tluvo imnlos pralioU-a to Ih> dommanl 
respective standard deviations are shown in Table T4- 1 7 and the ami>htnde eo.nponents arr 
ilh sinted in iMtjnre 3.4-1 T I Ik- predicted sound field, based on the calculated moda strue- 
aulle sev^t check nucrophone locations was co.npared witl, .neasu.d a.nphuu^s and 
phases and are presented in fable ,V4-1«, I he standard deviations are also shown o Ik d.l- 
ference between predicted and measured ainiditiides and phases. It .an Ik si.n tha t _ 
errors are larpe. 17.: dH to -4.: dlf in ainidiinde resiiltinp in a standaid deviation o .. 
dH about the mean, and between and -.<1 in ph^se. resiiltnip in a standard deviation 

of alnnit the mean. 
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TAMl.r v-» ts 

nii i i-Ki Nci- mawul-N i’Ki:nicTi:i) and miashki o 
inirr acoustk' 1 'ki;ssuki- 
(i*KiM)in i:n minus m!;asuki:i)) 

CONI ICUKATION NO,>l 


MicropliDiK' 

No. 

AinplitiKle ndla 
(dll) 

Phase Della 
(ileg) 

1 

2.7 

.1.1.4 

2 

-.1.4 

.1,1.11 

4 

1 2.2 

-.12.1 

5 

0..1 

17.6 

(> 

-4.2 

-0.8 

7 

..1..S 

16.0 

« 

5.6 

-8.5 


=.1.5dB 

*11 



Ad».liHoiKil predictions with more than three modes were not conducted for this case. How- 
ever. It is reasonahle to assume tliat the inaccuracies at the clieck microplioncs is due to the 
significant acoustic energy levels in some of the 50 modes not considered in tlie modal 
calculation. 


Mode Structure Configuration No. 5 


The fifth modal structure considered was that generated at Bl’h' 1 181.^ H/) by the wake ol 
a single rod interacting with the rotor operating at a .speed of .1400 rpm. An acoustic spec- 
trum at microphone location No. 10. for this configuration is shown in Figure .1.4-14. This 
case was included to evaluate the method for predicting all the incident propagating modes 
in the duel where the mode amplitudes are expected to be comparable. I he propagating 


nuxles at thi‘ trcMuency are the saute as those for configuration No. I . and these seven 
microphone !ocatit>ns for this case are preilictetl by llie MLl* to be at the same Utcations as 
in configuration No. 1. Also, the same three check microphones were used for this case. 
Tabic .1.4- 1 presents the measured amplitudes and phases at the ten microphone locations 
and the calculated modal amplitudes, phases and standard deviations are presented in 
Table .1.4--0, The mode amplitudes and ainpliUule standard deviations are illustrateil in 
Figure .1.4-15. The accuracy of the method at the three check microphones is illustrated in 
l able .1.4-21 in the same manner as in previous ca.ses. 
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l Aiii r .'.'I M 

DIM I Kl NC I- m rwi I N I’HI Din 1 1 ) AND Ml ASDKI D 
DUn ACODSl IC I'KI SSI'KI' 

(I’KI DK ITD MINUS MI-ASUUI DI 

l ONI'IUUKAriON NO. 


MuMophniK' 

(No.l 


Amplinulo Di'llit 
(JID 


IMiasi' Di’llii 
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S 
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10 


5..' 

O.l 

0.: 


-.'-'.7 

0.1 


TIk' I'irsi point to olisorvo is lluit ilio inoihil slmcUiiv is such tiuil the higher somul pressure 
levels ;ire itt modes close to eulotf. This result is not ullopether ime\|K'eted. It was shown. Cor 
example, by I'iekett (ref. lo) that for rotor-lurbulenee interaetion. the higher amplitude 
random modes are tliose close to eulolT. A similar model (used lor stator-wake-rott'f inlerae- 
lion) lor a steady distortion Ilial excites all the duct modes, sliows essentially the same result. 

.\ iH'ssible explanation for the 5. ' dlt ditTcrence between the luvdicled and measured smind 
Held at tnicroplione location N(V S is that the measured resultant amplitude at that location 
is tire smallest of the measured amplitudes at the three cheek tnicroplione locations. The 
reliability in inedicling low resultant sound Helds must itselt' be low as can be seen Horn the 
I'ollowing argument: If the low resultant sound Held amplitude is considered to be a small 
complex vector that is the re.sultant of a series of large vectors tiiat represent each mode, 
then if the large complex vec'.'>is are slightly perturbated in amplitude and phase, it is not 
difficult to see tlial large percentage changes can be experienced by the small result;. nt vector. 

In an efforl to obtain a better estimate of the sound mode structure in this case, it was de- 
cided to apply tl'^‘ matrix inversion method to determine the ;implilude and phase of the re- 
Hci led inodes that may liave a significant portion of the acoustic energy. The renecied (-.T 
(D and reflected (.T ID moiles were ineluiled .ii the selection due to an observation that t-.T 
0) and (.TlD modes had a signific;inl portion of tlie acoust c energy and with a cutoff ratio 
near one were expected to be tellected. .Again seven nuules were selected so ;is to group 
;iboul the dominiml nuules near cutoff. I he modes selected include; ( .T 01 rellecteil. (.T 0). 
(-.MD, 1 1. 0). (T 0). i.T m. and t.T ID tellected. 

.Acoustic pressure me;isuremenls ;il seven microphone loc.ilions ;ire iviiuired for the deter- 
mination of the sound field l>;iscd on the alunc mode g.roiip. I'or this study, five sets <0' seven 
microphone locations were chosen Hxnn the ten av;iil;ible loc.itions so that tlie matrix invet 
sion method could be applied five times. I his was done In running the Ml I' loi v.n ions com 
lunations vif se\cn out ol ten until ;i low eonditionmg numiu'i w.is obliiined. It was Idund 
that the live sets ol seven mictoplume loc.itions coiulilioned the ciiuation s\slem so ,is to 
\ ield .1 coiulil ioniuj’. numl'er I'cl ween I .> to I . '. 


oo 


Tho caU iihili’d ;inipliliuK's t>r llic soviMi moiii's I'ov oadi of IIk* Hvo cast's arc listed in I aide 
M- ” In addition, tlie mean and statulard deviation of the amplitude Joi eacli mode were 
caletilaled and presented in this table. It can be seen that the mean amplitude of the mcidont 
modes with a cutoff ratio near unity have the highest sound pressure level. ( oinparison with 
the mode structure obtainetl without includini* rellected modes shows tlrat the 0) mode 
compares well to within 0.1 dH. Alsti. the amplitude of the (-.V 0) mode is I0.S.4 dH lor the 
case that asMimed rellected waves were insipiiincanl aiul is within the standard tieviation ol 
the mean amplitude for the (-.1. 0) mode when two rellected modes were includeil m the 
matrix inversion melhoil. A final point is that the mean amplitude ol the rellected modes is 
an order of 1.1 db below the sanre incident modes that contiolled a sittnilicant portion ol 
the duct soumi field. I he inclusion of reflected modes in the calculation procedure did not 
exidain the observed discrepancy between the predicted and measured sound field at the 
check microphones. However, in each of tire live cases the accuracy at a check microphone 
deteriorated when the measured resultant pressure was small, and thus, realiirmed the above 
explanation. 


TABU .1.4-::' 


CALCHLATI D MODE STRUCTURli 
(SEVEN MODE INPUT TO MCP) 

REi Li:n ED modi: study 

('ONI K'.U RATION NO. .S 


Mode 


asc 

(-1.0) R 

(-.1.0) 

(•:.o) 

(1.0) 

C.O) 


'■“(OTir 

1 

X(> s 

lO.S.X 

101..1 

X3.6 

xn.i 

10(v4 

‘>1 .1 

T 

xs a 

10.S.(> 

101. t 

x.i.n 

8‘).0 

106.4 

‘>0.2 

1 

i)s 7 

io:.o 

lOl.X 

X().7 

‘>4.0 

104.‘> 


A 

U4 4 

107.(> 

101. .'1 

XX. 1 

‘>.1.7 

104.2 

‘>.1.0 

5 

ni,: 


101 0 

X4.5 

‘>4.5 

104.‘> 

‘>4.6 

Mean 

‘HU 

lO.l.x 

lOl.s 

X5..S 

‘>2.4 

105.4 

‘>2.6 

Staiutard 

IX'viation 

.1.1 

l.X 

0..1 

l.X 

1 ■> 

0.‘> 

1.6 


Mode Structure Conriguration No. (i 

The final modal structure considered was that generated at HIM- (l‘>'>5 II/) by the wake ol 
the single rod interacting with tiie rotor operating at a speed of 3740 rpm. 1 his configuration 
was chosen to determine the effect of a ten percent speed change on the method to deter- 
mine mode Structures with microplione locations chosen for the original speed and Ircpuency. 

T he same seven modes as in configuration No. 1 and No. 5 were generated in the duct, and 
at tlie frciiuencv of II/. tlie seven microphone locations as used in the previous conligura- 
lion were found bv the Ml I’ to yield a conditioning number of 3.1 . This value is comparable 
to those ..btained where the Ml P is used to select the microphone locations; so no 
in method accuraev were anticipated Data were taken at micnmhot.es No. I through Nt>^ 10 
•md are presented in T able 3.4-33. T he cc.lciilated iiunlal amplitudes, phases, aiu standaid 
deviations arc presented in fable 3.4-34 .md the amplitudes and ampblude standard 
devi.ilioiis arc illusiratcd in f igure 3.4-1 (v f inally, the comparison of predicted and measured 
sound field at the three check microphones is shown in fable .14 :5. 
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TAIU.K 3.4-23 OP POOR QUALlTYi 

Mi-Asuki<n niKT acou.stic rRiissiiui- 

(K1‘M- 3740 UPF- 1')‘)5IIz) 


CONFKJURATION NO. b 


No. 

DIB) 

Colierent Resultant Amplitude ^ 

(dyiics/cin") (lhf/iii.“) 

I’liiisc 

(dcK.) 

1 

91 8 

7.769 

1.127 X 10’'* 

49.3 

2 

95.7 

12.24 

1.776 X lO''^ 

2.3 

3 

85.9 

3.940 

0.5715 X 10''* 

357.7 

4 

104.2 

32.48 

4.71 1 X lO''^ 

42.4 

5 

100.4 

20.92 

3.034 X lO^" 

167.4 

6 

100.4 

20.89 

3.030 X 10"^ 

23.0 

7 

103.3 

29.32 

4.253 X 10-4 

172.5 

8 

102.7 

27.34 

3.965 X 10-4 

27.0 

9 

91.7 

7.707 

1,118 X 10-4 

358.3 

10 

94.4 

10.53 

1.527 X 10-4 

107.0 


TABLE 3.4-24 

CALCULATED f lODK .STRUCTURE AND ERROR 


CONFIGURATION NO. 6 



Amplitude 

Error 

Phase 

Error 

Mode 

(dB) 

(dB) 

(deg) 

(deg) 

(-3,0) 

103.2 

0.8 

36.8 

4.0 

(-2.0) 

101.6 

0.8 

209.2 

3.0 

(-1.0) 

97.1 

1.0 

357.0 

7.3 

(0,0) 

91.1 

1.5 

177,8 

15.6 

( 1 .0) 

94.7 

1.5 

175.9 

15.7 

(2,0) 

99.7 

1.3 

180.3 

8.0 

(3.0) 

97.3 

1.6 

49.5 

7,0 
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Figure 3.4^16 Sound Modal Amplitudes and Standard Deviations - Mode Strueture Configuration No. 6 - 
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Ursiilis loi lliis i;isc :iiv similiir (o Ihoso olisi'iwtl in n^nriiuiration N<». .S. altlioii(-li il can be 
M'cn tioni Tables .T4-l‘) and (bat llu’ ten iiereent rotor speed dilference between lesls 

residted in sipnirieant i haii|>es in measiiu'il amplitudes and phases at the 10 inicroi>hone 
lociitnms. Apain it is dirilciilt to assess (he accuracy of the method lor (his case, but it can 
be seen to be no worse than (or conllpuration No. .S. despite (he speetl clianpe. Apair., a 
luissibte eNplanation tor the ,S.4 dll dilTerence between the predicteil and measmed sound 
held at microphone location No. 0 is that (he measured resultant amplitude at (lf.it location 
is (he smallest of (he measiireil amplitudes at (lie three check microphone locations. 

.T4.7 Siiniinary of Test Results and As.sessmcn( of Impact of I-xtranemis Modes on Method 
Accuracy 


In the previous section, it was shown how the modal calculation metliod was applied to six 
didcrenf Ian duct modal strucUircs. In the first case, all incident propagatinp modes that 
could be supported by the duct were included in the modal calculation. In subsequent cases 
a subset ol tlie total number of inodes that could be supported by the duct was selected as 
predicted by the Tyler-SolVin theory and was assumed to contain tlie bulk of tb.e propapating 
acoustic energy. In the second and third cases considereil, this subset was enlarged to include 
nine propagating duct modes in an effort to improve the accuracy of the results. The ampli- 
tudes and phases of the selected number of modes were ciilculated using as input the meas- 
ured sound field at an equal number of microphones in the rluct as modes to be calculateil. 
t he calculated modal .structures were then used to predict the amplitude and phase of the 
sound Held at other duct locations where measurements were also made. 


|■■or the first test case where the structure of all seven incident modes was calculated and for 
the additional test of the second case where the structure of nine modes was calculated, it 
was .seen tb.at the re.sultant sound field was predicted within -H .7 dB and -0.8 dB for case 
one and within -0.5 dB and -l.‘) dB for case two. An alternate separate check could be pro- 
vided in each of these two cases, since the amplitude and phase of the single dominant mode 
could be calculated accurately based on an average of all the in-duct microphone ineasua'- 
ments. A comparison of the calculated amplitude and phase of the dominant mode using 
the modal calculation method with that obtained from an average of all the microphones, 
showed much better agreement. For this method of comparison, agreement was within 0.2 
dB in amplitude and 1.2" in phase for case one and within O.b dB in amplitude and 10.5° 
in phase for case two. Although the phase is somewhat off in case two. it is apparent that 
for the first two cases, the amplitude and phase of the dominant mode was calculated (]uite 
accurately using the modal calculation method. Since the prediction of the resultant sound 
field ba.sed on the calculated modal structure, however, is seen to differ from the measured 
sound field by almost two dB. it can be seen that using the calculated mode structure to 
predict the total sound field at other *luct mea.sim'inent locations is a stringent test of the 
modal calculation method. 


In case two, tlie amplitude and phase of the dominant mode was accurately determined des- 
pite, in the first attempt, the pre.sence of F' extraneous modes not included in (he modal 
calculation and in (he second attempt, 1 1 extraneous movies. In the third test ease eoiisiil- 
ered, the amplituile aiul phase of two movies. (-2.0) atul (-2.1 ). were calculated ami fourteen 
extraneous modes were present I vir this case the (vital soutul llelvl at other vluct locations 
was predicted within three vIH. with a stamlarv! deviativm of |..s dB. When tiie number of 


lO.f 


moili's iiK'Imk'ii in llu* liikuliiliun w:is ini'ivnsc’cl to niiu'. so that only si'vi'n uxtrani’ons modfs 
woiv pivsnit, IIk’ aivnracy in pivilictinH tlu* sound (u'ld at otluT duct locations was essen- 
tially unchanf.eil. Hie calculaletl amplitude ol the (-2.0) mode was 1.^2. 5 ilH lor the two 
mode calculation and |.M,I dli lor the nine mode calculation, and the amplitude of the 
(-2.1 ) mode was 1 2K.X dli for (he two mode calculation and 1 20.4 dl) for the nine mode cal- 
culation. file differences lielwen tlie (wo and nine mode calculations are fairly small (witli- 
in a standard deviation as olilained from (he nine mode calculation), indicatinj; that either 
calculation could he used to determine the structure of the (-2.0) and (-2,1 ) modes with 
reasonahle accuracy. 

l-or the fourth test case, where the .strucUire of three modes was determined, the ami>litude 
ot the total sound lieki at otlier duct locations hased on the calculated mode structure, was 
preilicted to be as mucli as 12.2 dB different from that measured, with a standard deviation 
of 3.5 dB. The significant difference between test case four and the other three cases was 
that in test case four, the number of extraneous modes supportable by the duct, but not in- 
cluded in the modal calculation was much greater than in the other three cases. In test case 
four, fifty incident extraneous modes could propagate in the duet in addition to the three 
used in the modal calculation method based on an application of the Tyler-Sofrin analysis 
to the fan design. Clearly the impact of these extraneous modes greatly affected not only 
the calculation ol the structure of the three modes, but also the prediction accuracy of the 
total sound field at other duct locations. 

The accuracy of the method when applied to test cases five and six is more difficult to deter- 
mine. In these test cases, a single rod was placed upstream of the rotor with the purpose of 
generating seven comparable incident modes at two slightly different frequencies. Modes 
witli a cutoff ratio near unity were expected to contain higher sound pressure levels than the 
other incident modes. It was shown by Pickett (ref. 1 6) that for rotor-turbulence interaction, 
(he higher amplitude modes are those close to cutoff. A similar model (used for stator-wake 
rotor interaction) for a steady distortion that excites all the duct modes shows essentially the 
same result. I'or the fre(|uencies considered, the duct was able to support 7 modes, so there 
were no extraneous moiics for these cases. The calculated mode structures were used to pre- 
dict the sound fields at three other measurement locations. It was found that in both cases, 
two out of the three predictions were reasonably accurate, but the third prediction was on 
the oriler of 5 dB different from the measured level. Because of the widely differing results 
from only three check locations, it is not possible to make a definitive statement on the 
method accuracy for these two cases. As explained in the last section, however, a possible 
reason for the large difference at a single location for each case, is that at those locations, 

(he amplitude of tiie total measured sound field was smaller than for the other two check 
microphone locations. When the resultant amplitude at any duct location is .small but cau.sed 
by the confiuence of seven modes, as in this case, and when the amplitudes of lhe.se modes 
are cottiparable, inaccuracies can be expected in predicting the total sound field at those 
duct locations based on the calculated mode structure. This is suggested as a possible expla- 
nation ofthe poor agreement between measured and calculated levels at the third micro- 
phone location. 
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In nitli.’r In oxpiori' I’urtlicr the cITcct on the method accuracy of extraneous modes, an 
order-ol-mannitude analysis will now he considered. The contribution of extraneous modes 
to the amplitude and phase of the resultant sound field at any duct location can he consi- 
dered e<iuivalent to a measurement system error. The magnitude of this error at any location 
is the resultant of all the extraneous modes at that location. The order-of-magnitude e.stimatc 
of the variance of the restiltant of the extraneous modes is obtained by as.suniing that the 
number of extraneous modes is large and invoking a statistical argument. If the total number 
of modes propagating in the duct is N-p and the number of modes used in the modal calcula- 
tion is N, then the number of extraneous modes is (Nj N). Since the mean amplitude of 
each extraneous mode as measured at the various duct measurement locations is 
and if (Nj N) > > 1, it is shown in Appendix H that so long as the phases of the extrane- 
ous modes are considered to be uniforml;^ distributed in [0, 27 t 1, the variance of the resultant 
ol all the extraneous modes is (A|.]^H|^)“ (Nq- - N), where the average is taken over all the 
extraneous modes. Thus, it can be seen from this analysis that the expected resultant ampli- 
tude of the extraneous modes at any duct location increases significantly as the amplitudes 
and number of extraneous modes increases. 

In summary then, the modal calculation method presented in this report can be expected to 
be accurate only for test cases where either the amplitudes of the extraneous modes are much 
less than those being calculated, or the number of extraneous modes is small, wliich is equi- 
valent to selecting the number of microphones used to obtain data during the test to be com- 
parable to the total number of modes propagating in the duct. 

In the next section the experimental requirements are determined for the application of the 
modal calculation method tospedlic applications. An approximate procedure for defining 
these experimental requirements is described and is applied to an example of a fan design 

that has been run on the NASA 1/2-meter (20-in.) Fan Rig. The general utility of the method 
also is discussed. 


3.5 EXPERIMENTAL REQUIREMENTS FOR APPLICATION OF METHOD 
3.5.1 Overview 

In order to apply the methods developed it is desirable to have, a priori, a procedure which 
for a given problem and desired calculation accuracy, would provide a guide to the number 
of microphones and measurement accuracy required. Such a procedure should be simple to 
apply and should not require input from either the MCP or MLP. I'o develop such a proced- 
ure. approximate expressions based on an order-of-magnitude analysis for tlie MCP inlluence 
coellicients were derived. The.se approximate expressions were used to generate a set of para- 
metric curves from wliich an estimate can be obtained of the retiuircd measurement accur- 
acies that will yieUI a itesired calculated modal amplitude and phase accuracy. Additionally, 
an apftroximate expression was tierived tlial relates moilal amplitude aiul phase errors to (he 
number of extraneous modes li.e., the difference between the total number of modes that 
can propagate in a dud at a given fre<|uency minus the number of modes (i.e.. microphone 
inputs) used in the modal calculation j . 
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i'lu' iiliji'iMivi’ ol (his pail of Ilu* piot’.iam is siiinmaivoil as lollmvs; 

In ili’vi'liip a piiH’i'ilmv I hat I’an lu‘ iisi'il lo piiiviiliv 

1 ) an l.•stilnaU‘ ol' iho nnmiHM oT iniiToplioiu's ivt|iiiioil to caUnlalo cohcivnl Can 
inoilal ^tl■nl•ln 1 v lo wilhin a spivifiod accuracy lor a specific applicalion. 

an csliinalc of the experimental accuracy reiiuiivil lo ensure tluil calculated modal 
amtilitudes and phases are wilhin a specified accuracy for a specific application. 


.1.5.2 Analysis for Inllucnce (.'oefneionts 
Objective 

To derive a set of on.ler-of-ma!'iiitude expressions for the intluence coefficients, equation 
(5. 2- IS) and to check them ayainsl tlio.se values obtained from the analytical and experimen- 
tal lest cases ilescribed in secdons .f 2 and .V4. respectively. 


.Analysis 


In section .1.2.4 and Appendix 1). closed form expressions were derived fora set of inllucnce 
coefficienls that are used lo deicnnine Ihe first -order errors in calculated mode amplitude 
and phase due to errors in Ihe knowledge of microphone location and errors in acouslic 
pre.ssure ineasnremenls. It is shown in Appendix n dial by assuming: I) the numbei of 
modes. N, being determined is large ti e.. N I and 2) the large set of deterministic phase 
angles associated with the complex matrix elements of the modal calculation algorithm can 
be formallv replaced bv a set of statislicallv indepeiuleiit and random phase angles uniformlv 
distribuled in |0. 2ir|, then oidei-of-magnitude expressions for Ihe variances 

of the calculated modal ampliludes t.A[^|) and phases (i,''m) can be deiaved in lerms of the 
variances in measured acoustic pressure ann>liludes and phases tujj.u^). anil Ihe variances in 
localing microphones in the duel tu“. ujr. 

I he assumptions u.sed in the analysis would seem inluitively reasonable so long as the num- 
ber of modes being calculated is large, but die jiislification required a comparison with ex- 
perience, which will be considered presently. 


Using. -Appendix l>. we ean write for (he expected values of die variances of Aj^j and 
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tii'M' .<nnpanso„N wn.- ,„;uk' i,Mni. ilu- wn analylioal losl oasos nsal in Soolion .k.' lo oluvk 
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IS a ,al,o ol tho amplittnlo at iho Onol wall dno to tho ,„|i; „u,do a,ul tho'.noall amllliuuio of 
lu usnlianl somul l,okl al iho ntoa.sniomonl looalions, Snpoiposoii on ihis plot is Iho oqna- 

‘’l-map,„|,ulo, M,uo mosl oh tho Oala falls within an o,Jo,-ol-ma!miliulo of Iho lino uivon hv 
uination ( ( oi„pa,i.so„.s lor phaso nioasnronionl onois anti mioiophono looalion onois 

S o\\ S„n,la, ,osi,lts. I i;’.n,o shows tho oo„ipa,iso„ of app,o\i,nalo a,ul ovaol \aluos 

mTm"? s"iV,"' ‘’I^Mln.lalotl nuulal phaso and n,oasn,otl pha.so. hor Ihis oaso, oqua- 

I o, . s- 1 ) im- tho phaso va„a„oo hasotl on Iho tutloi-olAnapnilndo anah sis ,otlnoos lo: 
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ligure 3.S-1 Comparison Between Calculated and Order of Magnitude Analysis of Influence Coefficient 



l igure .i.3-3 Comparison Between Calculated and Order oj Magnitude Analysis oflnjluenee ( oefficieni 
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Api>i'iii!i\ I'. I lk' ;m.il\si\ in Ai’i'i-iiiliN I' shows lhal (lii’ vaiiiiiu'o ol tiu* tvsiihaiu ainpIiliuU' 
lit Ilk' I’xiraik’oiis iinuii's is piopm lional (o llio mraii stptari' ami'litiiih' ol Ilk' oMianooiis 
iiukli's link's Ilk' iniinlvi- ol' i'\liaik'i>iis mi>tii'S, I'llal is: 
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(N,.-N) 


13.5-4) 


wIkTO N |- is tik’ lolal mmilk'i ol' inoilos snppoi laNo In lltc N ^ s IIk' miinhor ol moik's 
as.snnk'ii (o roiilain Ilk' hulk of Ilk' aionslir ciu'iny, aikljl'm" Ap^" { ‘ is the rms vahio ol 
llu* oxtiaik'oiis nukU’s al llio nu’asuiviik'iU locations. 


Now. the rcsiiltanl ainpliUklc of the cxliancoiis inoilcs can be cmisiiteivii to be equivalent to 
an eiror in acoiislic ineasinvinenl. Thus, as.suniintt perfect measuring system, (he variance 
of the inessure ineasiiremcnl. U|^ in the sensilivity analysis (equation 3 . 5 - 2 ) can be consideivil 
to be a pressure nieasureincnt error causeil by ihe number rrf extraneous imnles. Substituting 
into equation (3 5 -.?) yicKIs eipialion ( 3 . 5 - 5 ) 





( 3 .. 5 - 5 ) 


In etiualion ( 3 . 5 - 5 ). can be considcrcil to be the square r)l' the siiUial to tu>ise 

* til '^l in 

ratit>. since I'l^A ^is (he .icoiislic .uupliluiie al the nicasurcinenl localion atiribulable to Ihe 
imules lhal arc due lo be ealculaleil (ami winch arc .issumcil to contain the bulk of the 
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acoustic energy) and y' I • "^M 


is (he overall conli ii'ulir>n due lo the exlrancous modes. 


I xperictue gained fioin ihc cxpei imeiilal portion of Ihe lU'ogr.iin demonstrates Ih.it Ihe 
mellukl aecuraev detenorales as the numlk'i of exlrancous inodes increases, providing quali- 
tative substantiation of eqii.ilion (.'.5 5 ). I he ext'eeted error ol the ealculaleil mode ampli- 
lude in decibels e.in be obl.iitied fiom 
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whereupon, using equ.ilion (3 s st. 
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.<•5.4 .\|)|m)\iin:ilo (oi Di'liiiilioii ol Kt'(|iiiivil Mr;iMiivnuMi1 AiviiiiU'y 

Tlu' niili ‘1 ol in.n'iiiliiilr I'Nlinuih ^ ol ilu- iiiHiu'iKi' lOi lTh u iils i;!!! also 1 h“ iisi\l lt» (k'U'rmiiU' 
llii' oxiH'iiiiU'nlal ;u\Mii;u \ u-(]mionu'iiis lo I'lisiiro n.’sulls uiihiii ,i spooirii-ii I'lior ranj’.i' lor a 
siH\irii' ai'i'liialiim. I lu- loiimii'ii ai mislii' i>ivsmiu' nu'aMiivmi'iil aivurav v will lu' I'on.siiloivd 
I'iisl. .\sMimiiiti oiil\ aoouslii' |m\'smiii‘ ami'liliuli' nuMMivi-moitl I'liois, i’i|iialioii (.kS-l) bc- 
itinu’s; 

£'.15=!] L 

il^ ' 2. !\;^i 

I'hi' o\|K'i‘toi] onoi in nunial aiiipliliuU’, I'xpivssi'ii iu ilocibols. can be slu>\vn in a similar man- 
ner to that pivscnlcd in ciinalion (.Vx-o) to be: 

Ad8[^- Za L»3„ [ 1 

In aniicii'alion of itcvelopinu a proccilinv I'oi esiimalin.e llie instrnmcntalion rcqiiiremenls. a 
desiivil accinacy tor masiimim ctioil I'oi ilic calciilalcd iilojal amplilmie taiiil possilde 
phase) will be assumed iiixeii. I bis same error will Ih- used lo esliiiiate Ibe number of allow- 
able exlraneous nuuies (or. b\ appliealion ol l ii'me .v.s-.'. llie luimber ol' mieropbones required), 
the allowable error m pressme ampliHub- and phase measiiremeni. and the allowable error in 
inieroplume location Hie inleipieialion ol ihese allow able errors is diseus.sed in section ,£5.5, 


Proeeedinir we can combine eqiulions i.£,s-(i) and (.£5- to pive; 
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t.).5-S) 


.Since N | and N have alreaiiv I'cen .issimied iaii’c, it will be tnitber assumcil that the mean- 
square ampliliide ol ibe measmed piessiiie. is ol the same order-or-maimitude as Ibe 
sum ol' llie expeetv'd me.in square resuiianl due onl\ lo (lie exii.meons modes and Ihe e.x- 
peeled mean-sqiiaie lesnll.ml due oiiK lo Ihose modes Ivmi' ealenl.iled. I lial is: 
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A') /^, 



(.£5 ‘M 


t learlv. ibis appioxim.ilion w ill be leaMiiiable tl Ihe eoninbulion liom (he exiianeous mode 
is mueli smalh‘1 lb, in the i onitibiitioii liom the modes I'eme used m the modal ealeiilalion. 


I I I 


This will lU-li-riornli- lunvt'vcr, if tho two conliilnilimis luv coini>:ir;ibk’. Suh- 

sliuitinn oiiualion (.V5-‘)) in filiation (3.5-S) yields an equation I'rom whieii a set ol paia- 
metric curves can be obtained demonstratinj?. lor a range of signaMo-noise ratios (i.e., ratio ol 
the amplitude of dominant modes to the rms acoustic amplitudes ol extraneous mode . when 
extraneous mode contributions and pressim' measurement errors cause an equal error 
calculated mode amplitude. This set of parametric curves is shown in I igure .T5-4. Also, it 
will bo seen that, following a similar analysis, the same curves can be used to relate extraneous 
modes and measurcil phase errors in tlie same manner. Thus, by way of example, il the con- 
tributions from the average of the extraneous modes are about 30 dB down from the domin- 
ant modes, and if the ratio (N^ ■ N)/N takes the value 50, the error introduced m ihe calcula- 
tion of the mode amplitudes will be equivalent either to a 1 .7 dB error m amplitude measure- 
ment or a 13° error in phase measurement by the microphone system. 



o: 

o 


o 


a: 

o 

a: 


{/) 

< 

I 

& 

n 

u 

CL 

O 

v> 

< 

UJ 

S 


rif;urf T Itvquin’itu'iih < >fl’n-s«irc Mcasiin'inritts 


ii: 


V 


I'Mii!’ .1 siinil,ii .iii.ilwis, .1 ol I'.iunu'lnr i i.m Ih' )'viu'i . ili'il iti.il Ir, loi .1 

I'.ivi'ii Mt'ii.il I.' lU'iM' i.iliiv wIk'M I'MuiK'oiis moiL' vimU 1 iImiiu>ii> iiiui asi.il aiui i iikiiimi’ii'iitial 
mik'ii'i'lu'iu' lo»’,ilti>n oih'iN iMiiv .m i’i|u.il I’lioi m ilu- i ,iK iil.ilk‘d uiotlk' .imi'liliuli- .iiul j'Ii.imv 
I'Iu'.huKmn Ioi U'l.itni!'. .im.iI UuMlioii I'liots s\ iili tiu’ iuiuiIhm I'l oxh.iiu'oiis nuul-.'s i> a> 
fiillows: I'k’ilVil puMsiiivnu'iii aiiJ oiu>is oiilv in avial liH-alion. oitnaliim 1 ) 

Iw’oiiU's I'oi llu’ Ptioi in k.ili. iilaloil inoilal aini'liUnU' 
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It' I \ l*’ ' i> I'lmnnalod t'loin I’n nations t..'..'' aiul t.v > 10) a ivlation is oblaiik\l Ih‘ 

tn 

IWkVn Ilk' pmn m axial inivioplioik- Kvalkm ainl itk' nninlvt o!\'xttaik'ons nunk’s ilial 
i*aNp l1k‘ sank" onot in iMkailaloil iikklal am)'lilikio 1 liiis. 


' " t ( V ) 


In oixkn to I'nisiip tlk- oialoi of mat'iiiliklv' analysis, it will l\’ assuiik-il tlial aikl A^l ^ aiv 
stalistii'ally iikli'|k'tkk’nl. I tu'ii '' 
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w lu'ii'iipon. 




I 01 a i.iink' ot signal to noisp i.iiios / 1 -I, ). i'>)ualion ( v> 1 ’) is iisi',| |o ok 

lam I I inn p ' > wlmh is .1 plot »'t llio nonnali iki axial kkation piioi \i-tsns ilk' lalio ol 
I’xliaik'oiis iik'ik's lo I'ak'iilalk'J moik's Siippiposikt v>n tin' xaiik- tii’im- aii- tik' I'liviinPvi 
I'lilial loi'alion piiois ximmis tIk' talto ot pxtiaik'ous inoik-s to l ak ulalpil inoik-N tii.it was 
ol'laitk'ii I'x a ^lmllal .iiuIxm' 


It!'. 1 1 IV i;m tiuis bv iisvil Id sIidw. lor ;i niiiv t>rsi!!ii;il lo iioiM- nilios, whvii i'\ii;mvDiis 

Muuli’ i'Diilrilnilions :iiul a\i;)l aiul i-iivmuri'iviiliiil mii-ropliDiu' loialiDii i-i iors j-iw lisi- lo an 
viliial viTor in llu- valviilalod imuK' ainidilinlo. A similar ivlaliDiisliip iDr mrors in radial 
iniunpliDiU' lovalion is iu>i ivadily oblaiiii'd lu’vansv ilu* approprialv inlinviu’v viH'lTii’iniis 
VDiitaiii radial tlmivaliws of Hvssvl rimvlions. Based on e\pviiviue iroin ihis pioj-.iani, how 
I’vor. Hie errors inlrodiieeil inio liie nuulal ealeidalion by radial mieroplKnie hualion errors 
are always imieli sinalier lhan tlu)se due itr ilie other measurement aiul loealion erians eon 
sidereil here. 

It is now possible iisinjj l••i^;ures .<,.S-4, and to deim>nslrate a proeednre for esti- 
mating the e.\periniental accuracy that wouhl be ier|uired in applying the method lor de- 
termining cohercnl Ian sound nuulc siriiclure. 


ACCURACY REQUIREMENTS OF MICROPHONE LOCATIONS 
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IIATIO / 


NUMRFni>F 1X1 HANLOUS Moms to I’llINCII’Ai MOOFS, 
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tiiinniinrs .y,/;, i; 
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vs fi III I'liMi'tlmo Im Di'niiinj- I UiMiiiiii’iiu’iils 

Siiuv I I ;iiul > "i-ii' I'iiv’il on onlx'i ol min'iiilmK' .ni.ilysis, ihoy ^ ;il Ik'M 

hr iivii .IN .m .iiM'inMin.ilr iiiiinhris ol mu ioi'hoiu's nml nuMMiiviiiriil :iml 

posilion .u-air.K ii'N )n lion »'t ;iiu |'noi lAinTiriu r. llu- I'lOiYiliiii' lo lu* ilrsi rihnl riiii lu- 
UM'il .IN .1 I'litiK- in xlrlinm;' ihr insli iiiiu-nliilion iviiiiil'.'nu‘nh.. Oiuv ;i li,’sl 1 i;in Ivrn roini'lrU'il 
aiul Ihr NU’P has lu-rn iim'iI lo raKiilalr imuial ainiililiulrs aiul I'liasrs. ilir inllnriuY roi'llirirnls 
(llial aiv also ralrulalrd In ihr NU IM, ran hr nsnl lo inianlilV nioiv arniialrU Ilk' oriois tn 
rairnlalril inoilal .uni’liliulrs. 


I hr I'lorrilmr lo hr ilrsnil'nl slioi lK a>nsisls ol' sirps, l-irsl a valiu' of Ilk' ilrsirnl ao- 
rurary. or ma\iinuin mot. in llu' raK'ulalal inoilal stiiK-lmv is srlrrlni. I hr srroiul slop 
tuouiirs an rstiinah' of ilir lolal miinhrr o!‘ niirrophonrs ivipiiivil to proviilr llu- ilrsiivil 
rairnlalril inoilal ainplitmlr ami pliasr arrmary. assiiininii thriv air no inrasuivnirnt or posi- 
lion rrrors, I hr Ihinl slrp inovidrs an rsliniair of ihr allowal>lr rrror in mrasurnl anipli- 
liulr aiul phasr. ronsisirnl n ilh ihr iviinirrinrnts of slrp onr. assnmin.n all niirrophonrs arr 
prilrrlly loratnl ami Ihrrr arr no rxtranroiis moilrs. Stop 4 providrs an rstimair ol tlir al- 
lowahlr rrror in ihr mirroidionr posilions. ronsistrnl with thr rrqnirnnrnls ol slrp 1. assiitn- 
inp Hull ihrrr arr no rMranrous modrs or instrumentation nirasurrmri't rrrors. In step the 
allowahir rrrors in mrasnrrd prrs.surr ainplitmlr and phase are roinpared with thr estimated 
instrnmentalion sysirm inrasuremrnt rrrors that would he rxprrlrd in eondurtinji thr test. 
Similarly, thr allowahir position rrrors ralrnlaled in step 4 arr roinpared with the aetual 
position rrrors that would hr rxpretrd in loratinn thr niirrophonrs lor tlir lest. IT the aetnal 
rrrors arr an order or inaiiniU.dr (or morr) less than thr allowahir rrrors. then thr arrmary 
srleelrd in step 1 ran inohahly hr arhirvrd by use of the luinihrr of mirrophones drlinrd in 
strii If on thr oihri hand, (hr arliial rrrors arr romparahlr to thr allowahir rrrors that 
yield ralridalion rrrors of ihr modal rorlTieirnIs on thr order of (hat speeilird in step I. then 
these ralriilalion rrrors will arrmnulair from rarh sourer, and thr total error in the ralrnlaled 
modal aniolitiidr and phase will hr larger than that sprrifir.i in step I. In this rvrntuahly. an 
assessment of Ihr modal ralriilalion arruraey ran hr ohlainrd only after Ihe test has hern 
rondurtrd. 

In Ihr followini’,. Ihr above steps arr outlined'. 

Step 1 : SeIrrI drsirrd arrmary for modal ainplitmlr and phase ralrulalions. in dH and degrees. 

Slrp I or ,i given diirl grotm iry and sound Imiurnry. the total miinhri ol modrs suppoit- 
ahlr by the duel ran hr ralrulatrd’. l or the parlieulai arrmary in ralrulating the dominant 
mode sprrifird in Slrp 1. the nmnhrr of allowable extraneous modes ran hr drlrrmined 
from I'igmr asa funrlion of sig.nal-to-noisr ralio^ idrfmrd as the ratio ol (hr ainiMilmlr 
of (hr prinriplr modrs as mrasurrd a( (hr dur( wall (o (hr mean aiu|ditmlr of thr rx(ranrous 
modrs as mrasurrd at (hr diirl w.ilh. I hus. (hr numhrr of allowahir rxlraiiroiis modrs lo 
arhirvr (hr drsirrd arriirar\ in Ihr modal ralriilalion r,m hr drlrrimiird. Iroin whiih thr 
mnnhrr of mirrophoiir mrasiiiriiimis that nrrd lo hr made also ran hr ilrlrrmiiird by sub 
liarling (hr nuiulH-r of .illowahlr rxdanrous modrs (roin thr (o(al mmihrrs of modes sup- 
por(ahK‘ h\ thr diirl. 
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.V5.(> \p|>lir;iiinn nC Moili- (':i!i'iilaliii}% iVoiTihiri’ in a N ASA I an Ih'sinii 

I n i'i\n nil' an I'vampU- of llio api'lit'ali»'ii of llio pioiViluiv Inf iliMi'iiiiinin)’, llu' rvpi’rinunilal 
ai’i’uraiA loiiiiin’iiii'iils. M'l’lioii .'.S.S. usini*. an fan licsipn dial has Ih’i'm nin on llio 

NASA sOi iii ( ’() iiu’li) ilianuMor fan lii’. 

Applioatinn 

l o ilUisiraU' tlu- applioalion of Iho proa'diin' for diMonninin}’ llu' rcquiivd csporimonlal 
aioiiraoy for a spooifio api’licalion. a partiailar fan design will he used that has been run im 
Ihe NASA SOem l ’0 in.) diameler fan rig at Ihe Lewis-Research ('enter. I'he particular fan 
stage cho.sen, designated rotor .s.^, has a 1 bladed rotor. 25 fan exit guide vanes, and has a 
design speed of S020 rpin. At a typical approach speed, the Tyler-.Sofrin theory would pre- 
ilict that at twice blade passing frciiuency the rotor-stator interaction will generate five modes: 
(-^.0). t.s.l), (.S.2). (5..)). t.5.4). 1 \ pically. however, the duct could support at that freipiency 
on the order of one hundred incident aeoustie duct modes. 


If the desired accuracy for t'lvdicting the dominant modes in the iluct is set at approximately 
±2 dU. (Step 1 ). and if it is assumed that the mean amplitude of the modes omitted from 
the calculation procedure ti.e . Ihe extraneous modes) are about 20 dB itown from the 
dominant inoiles. then it can be seen from l-’igure that the allowable number of ex- 
traneous modes has to be about 1 2. I'his ipiick analysis illustrates that, providing the assump- 
tions made arc reasonable, a ' 2 dB accuracy in Ihe calculation of Ihe dominant inodes can 
be achieved only if 88 li.e.. 100 12) microphones are used to obtain data iSte)' 2). Note 

that with the same assumptions, a (d ilB accuracy in the modal calculation could be obtained 
with alHHil hO microphone measurements. 

It the rolor-sialor spacing 1t>r this fan design w;is reduced so that a strong interaction tone 
was geiiertled. and if care was taken lo control inllow distortions, it is conceivable that Ihe 
mean amplitude v>f Ihe extraneous inoiles will be 2.5 dB lo .)0 ilB ilown from Ihe five domin- 
ant modes ilue lo the interaction. In this case, it can be seen fiom I'igure ,V.5-d that an accur- 
acy with ' 2 dB coulii he expected with as few as five microphones. 


I’oi purposes vd proceeding with the illustration, it will be assumed for Step 1 that a ' 2 dB 
accuracy in the moilal calculation procedure is desired and that Ihe mean amplitude of ilu’ 
inodes is alnnit .U) dB down fu>m Ihe tiominant nuules, I'rom Step 2. using. I'igure .V.s-.t. it 
follows that five microphones should be used. I he allowable error of Ihe acoustic measure 
menis (Step .)) can be iletermincd trom l igure .V.^-L Lri'in the assumption slateil above. 


N| - N 
,_L — ) 


I'L Lrom figure .t..s-t it can I'e seen that Ihe allowable error in pressuie ampli- 


I I 


tihle measurement should be 1.1 ilB. and tlie alK^wable error in pressure phase measure- 
ment is al'oul eight degrees. 


IJr'nTlItTr." A?"'”!'' "'r'"'"','"’™ "'■ >"'■ «" i" tM* is „l,„.im.,l ,S,.,, 

k iK n I i' ui*'!'*’ “•‘‘'“"options s(:i(i‘il jihovc, jiiul siiu’c m is fivo iml 

i SnrV II '™ II" l"ini> consklcv.l, ll cm he seen Iroi,, 

!, M- 1 1 , '^i' »“<'»» 'Is- em.r in i.laeinit niieroplinnes nt speeirie eiienmrer- 

«Umr!s?,‘r',h ®l‘n’ '"““"""'"I i"'-l l*'“■•"■nl lire e,.,„„im.,l will, 

•o M .1 . / ‘"’ll pliiccmont errors (hiit wmilil be expected when 

f I n ^ location errors met In practice would be much 

Lss than the 0.1 cm in axial direction and the 1.6 degrees in circumferential direction. Also, 
tl.,' measured acoustic pressure amplitudes and phase errors sliould be able to be held within 
values wel. below 1. 1 dB and 8 degrees, respectively. Thus in the situation described in this 
al^d le’ d ■»‘^"-n'-ncnt and location accuracy requirements are «blc 

a d the desired modal calculation accuracy of +2 dB (Step I ) should be achieved If more 

Sum method, the cal- 

e nation aeeiiraey will improve. Ihis improvement would continue with the addition of 

^ the 

Of^it number P-oex-dum that can be used to provide an approximate guide 
t le number ot microphones and the measurement accuracy experimentally reouia-d for 
a speciUc appLcation. it is recommended that a test of the modal calculat o meSL^ 
hm fh t i“" \ NASA 50em (20 in.) diameter rig. The nuniber of mie.ophoncs use^I dur- 
1 I fe the test and the number ol modes u.sed in tlic modal calculation, could be increased fur- 
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himn.i' (liiN I'loi’i.im. twn sOuiu's wi'ic lot ihi' t'.ilriil.Kiiui ol llu' iiuiil.il slnu luu' 

III llu- iiiK'l 111 .1 Iiiilu’l.in I'MiiiiU' iiMM!'. Ill iliu'l iiuMMiK'iiu'iil ;is lIu' iMlrtil.iliiui m]Mil. I'm llu- 
liiNl Nliiih . i'si'i-nMU-iil,il .iiiil 11111111111.1(1011,11 I'uui'iliiii's .iiiil two lOini’iiUM ('io)'i.mi-. wou' 
ili-\ i‘lo|i|-il loi llu- I'.ili iil.ilioii 111 (111- .iini’liliiili-s .nul I'li.iM-s ol' llu' volu'ii'til i’iO)*;i! .iliiii' 
soiiiiil moili-s In .iililiiioii. .Ill i'\)U‘imu'iil.il |iioi'.Kmi w.i> iiiti lo pioMilo u'lilu-.ilu'ii .>1 llu- 
I'loii-iliiu- I Ol llu- otiu-r sliiiK . I'loi-i'ihiu’s woiv uli'iilirii'i.1 loi llu- ili‘U‘iiinii;ilioii ol' ilu- i-\ 
lu-i li'il ,inn'liiiiili-> ol llu- iliu’l moiU's u'siiliiiij'. I'loin .i niiulom iioi.st- lU-iuTiilinj', iiu'i-ii.misiii. 

In t'.u l. liu'M' I'loiViliiii’s ,nv surtu-ii'iilly )u-iu‘i.il .nul I'oiiM In- .ipiilii-i* (o .my iliu't smiiui 
tu-lil I lu’so iiu-i .ill pioj'.um ivsiills will In- ilisi ussoil si'p;ii;iti-ly I'm i-iu li ol' llu’ Iwo viisos 

•j.i rom-ui Ni MODi- si iuu ruRi-s 

l lu- v-iinipiil.ilion.il piov-i’iiiiri‘ lo i-.iU til;t(i' tlu- ;impliliuli-s .iiul pliasos ol (lu* i-ohoivnl iiuulo 
.slnu lmv in llu- iiiK-l iliu l ol a spooilii- l.m ilosiim is h.isal on iwo l omptiu-r pioi’,iam.s. Ilu' 
Minoplimu- I oi.ilion I’loiuam iMl aiul llu- Modal (.’.iK iilalion I'lop.iam (MCTV I'iisl. 
in lliosr i.isfs w hoiv (lu- niimhov ol' piopaii.ilini’. modos is luvaliT than tlu' mnnlH'r orniiiio 
pluuu's lo 1 h- iisi-d. .1 siil'si ‘1 ol' inodi's is si-li-i-ti-d liiim .ill tlu- ptopap.aliii!', iiiodos siiiipmlaMo 
In llu- inli-l diu l lhal .no oonsidou-d lo oonlain llu- bulk of llu- aoonslu- oiu-n'.y in llu- diu l 
sound Ik-ld. riiis si-Uu lu'ii l.m lu- b.isod on an appliialion o! tlu- l ylor Sol'iiii analysis. 

Willi llu- iuiiiilii-i ol modi-s lo bo o.ilonlalod ulonliliod. llu- Ml I’ is usod lo dolormino looa 
lions Im .111 oipi.il iiumboi of mioiophonos willim ibo duol llial oiisnios nnmotioal slabiliU 
ol tlio Mibsoiiiionl mod.il o.ilotil.ilioii .So loiii’. .is spooit'io oiiloria aio s.ilislioil llial oiismo 
iimiioiio.ll sl.ibililv m (ho mod.il oaloulalion. Ilio mioioplunio loo.ilions aio losliiolod only to 
a\i.il .nul I noumtoioiitial w.ill nu-.isniomonls I'or many sitnadons. lunvowi . lo piovulo snl 
lii 10111 sl.ibilil\ . mpiit mo.isnioinonls aro loiiniiod I'lom all llnoo dimonsions wiiliin (lio dnol. 
and 111 llu'so oasos Ilio Ml P will solool siioli po.silions. I lio MC'P, nsinp. d.il.i liom tlu-so nn 
oioi'lionos .is mpnl oaloiil.ilos (lio amplilndo and pliaso ol oaoli ol'lho nuulos oompiismi’. tlio 
olioson sol ol modos. 

In .uldilion. Ilio Ml P oalonlalos a sol ol inlliionoo ooortioionls dial oail bo iisoil lo o.iloiil.ilo 
llu- liisl null I 01IO1 III niod.il .nnpliUido and pliaso lli.il would losnil lioiii onois in llio 
knowlodi'o ol ilio nioasnimp. looalion and in Ilio moasinod pli.iso .mj’.tos and ptossino ampli- 
Indos I'liialb . il Ilio sl.itulaul dov ialions duo lo inoasmomonl aiul mioiophono loo.ilion 
oiiois .110 mpnl lo llu- Ml'P, tlu- siaiulaid dovialions ol llu- losullni!’. onois in oaoli ol llu- 
o.iloiil.ilod modos is oompulod. 

1 ho inolhod loi doloiiiimni!'. i ohoioiil niodo stiuolmos was ohookod out iismp Ion avbili.nily 
ohoson aiuK lio.il losi o.isos. .nul siibsoiiiiontU . ni a sol ol ovpoiimonls lo dolonnmo si\ dit 
loionl modal s| i lu linos I'onoi.ilod m llu- miol ol llu- Pi.itl X W Ininoi Anoi.iU .'>om ( Id 
nu hi loso.iu h Ian lii'. I \poiimoiil.il piooodiiios lhal woio osublishod on iho b.isis ol tho 
.iii.ib lu ,il inolhod woio also olu-okod init diniii)'. this sol ol ospoiimonis 

UoMills limn llu-so ovpotimonls domoiislialoil ih.il llu- inolhod loi doloiinnnii!' oohoioiil 
niodo siiiu linos IS siiil.iblo loi lolalnoK simplo inodo sliiu linos w lion Iho inimboi ol mu io 
|•IU1l1os iisi-d lo inovido d.il.i mpnl to llu- Ml P is oomp.n.iblo lo llu- tnmiboi ol nu utoiil 


1 l‘> 


m.uks m llu- .Im I \n o^w . mm-, llu' smu Un.' .>1 >vvc-n m>..W ss.is ul.'ndiu. 
uMm- so.-n m.nophoiu-v .uui .( «..s .U-lnimiuM tlu.l tlu> .lonmu.,1 uuui. ^as . al. iil.iUal 1. 
vMlh.i^ .m .ua-iiu. N ol O. ' m ..mphlu.U- a.ui 0.‘> .lo'uv •« I’tuso. I.. ott.n . asos. «1 k-u- |Ik‘ 

siuu (uu‘oiupt.Mmu- .muk•^^^as>U'^•^m,Kal.an.i (lu- nm.itv, ol rnu uM^s "Mai m llu 
w I \xason 111.' .M.U'i ol iMU'liall llu> (olal miinlvi o\ iiu uk-nl mo.K s ilun rouU pioi'a.M . . 
tlK> suiunnr was.K-l.-imnual will, an a.vuuu x ilial alloxval llu- somunu-Kl lo K imalu lyl 
in xMllmi a slaiukmi oiiot ot l o all al ollin .hu t kvalions Kisn.l on llu' .akiilaual 
smu-luK’ In \x-l anotlun naso xxlioiv tho stnu tmv nt (Imv nuuk's was .kMniimua iiMii} I ux 
;,„,,oplu>.u‘s iM.l wlK-.n llu- a.uM .onkl support >0 piopapatinp. nu uk'iit 

M otlun aiK-1 looauons vxas poorly pivvliotoa xutlun a sla.ulata nro, ol o,> aU lus.a 

on till' oak nlati’il inoik' slnu tiiro 

n isi’ii on llu-si- ovporinu'utal aiul atiaK tual li'si vasos a.ul on a oonipk'nu'ntarv analvtual 
suuU .1 was louna that tho aoouraoN o( tho mothoa aopo.uls on tho nnmho. a..a axoiapo 
:„„p,.ioao ot tho oxtranoous moilos. (i.o.. tho ihrioronoo ho.woon "V';'"'';;;:"'' 

pot tahlo In tho iluol . aiul tho numivi ol nioilos usoa as input to tho Mr I I as w i 1 as h 
moasmomoiil aiul i>osition aoouraoios An apiMOM.nalo tnooialmo was aovolopoi to lu p 
ruulo tho oxponmontalisi in ulonlilx mo. tho niimlvt ot niiorophonos that woiikl lu uiiuiui 
to aoluoxo a aosivoil aoomaox . Smoo tho ,MOOoaino was hasi-a on an oiaor ol niap.mliulo 
iMmIn SIS» ll sliiMlIil I'l* use'll I'uK * 1 ^ i*,lluU 

4 K\NIH>MI Y i;l Nl K \n l> MODI Sl KUl Il'Ul S 

Hx usmo Iho I'lopoiU that am aooustio liolii m a aiul riaiulom or oohoiontt ai am mstanl 
oMuno must oxlnlul 'a po.ioau ilx in Iho anonlai oooiainalo. it was slum n tlial a '"‘"'•'I 
Powoi .Spoolial Ik-nstlx Uinolion tl*sm oonki ho ponoialoa that ilotviuloa on lioiuioiu x aiul 
, 0100,01 vahios ol Iho wmo roi inoilol niiinhoi in Iho oiiouniloiont.al atioolion 1 liisaolmi 
turn ot Iho Moaal I’SH .ilso pormiltoa inloiiuolation ol laiuloni aoouslu tiokls hn Iho anmilai 
auols m loims orsiunninp. nuulo ooiuvpls snmlai to thoso oxolxoi! lot oohoioni aoouslu 

I lokls 

l ivo M'Pioaohos loi iiuasuiinp Iho iiioa.il I’SP m llun annular auols woio aomoa I ho lust 
ovolvoil iron. Mio aorolopmonl of Iho nuulal PSP oouooiM ami lo.tuuos oioss spoolial anaU s.s 
of m.o.o, . 1.0110 s.onals fo. its uu|.loiiionlalio.i \ sooo.ul aptuoaoh was also aomoa. usinp. 
Ohasoa a.iax inoasuioiuoiil ami aal.i loauot.on oo.ioopis I li.s sooo.ul mollioa aoos not io 
iliiiio OIOSS sivotial piopoilios loi Its uuplomoiilaliou aiul shoiikl Ihoioloio ho oasioi Io um 

\u oullino ol tho a.iax mothoa .s as lolknvs 1 o. a spoo.f.o a.u.ulai a.iol poouioln a.ul 
f.oiinonox lam-.o o! inloiosi. a minuuuiii tuuiiho. ol nuoiophonos is solooloa iii oulo. to s.ii.s 
tv ilio Nv oiiisl oMloiio.i 111. fio.iiionov ai.a nioao iiu.uhoi sp.ioo V I hoso m.oiophoiios au 
pbooa oi|U. a.sta.it 1.1 ono axial I'lano on Iho iluol wall Pmmp. Iho losl. Iho oulpuls 
1,0.11 Ihoso mioiopl.onosa.o tooo.aoa on analoo lai'o a.ul o. suhso.|uontlv a..'.li,oa llu 
Munul lioki .slhoiia.so.i.uu.aloa into aooi.st.o w avos w.lh spoo.lu spin lalos I hi a.s, , iiu.n 
,,„on oan ho aoh.ovoil who., t-iooiss,.... iho a,. la hv .usl p..ssuip iho s.p.ial lion, oao h mu lo 
t.hono Ihioiioh a iiai.ow ha.ui Hilo, l.' solool s.puals ,.i a piosoi.hoa lio.|uono\ . ami ihon s.ih 
soiiiioi.llv aolavm.'oaoh s,p.n..l oulpiil In a .-non ino.o.uoul ... tuiio ..ml Iho.i .ivo...piiu- Iho 
Mjril.lK lo \ U-lil .1 OllliUlt 
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III I’lU'i t. Iliis |n«n\‘'cs I’lili, mi l's .MUisitr w.ivi-s ol ii ('.ivni I'u'iiiu’iu v ii*l.nin(’ ;il ;i s|Vi’ihi 
Nj'in Kill' .iKmi* llu' iliii'l w.ill .iiul suppivssos tlu* iii|nil lo llu' sijuiiil lioni w.ivi's loiiitiiij*. al 
otiu'i si'iii uli's. Shut llio spin utc is I'lvon by n', ni iwluMV u’ is Ibo uilian lu’ipii'iu y »M 
iutiMi'si aiiil 111 IS du’ i iivumiViviUul nunli’ niimiH’i). oiu'i* Iho spin lali’ or ili'I.n linu’ has 
tvon si'li'i'U'it. thi' oinpul sic.nal consists i>l I'onliibtKions only rroin aomisiii' wavos w iih a 
spooitu' I'lii'inmi'ii'nlial nuuli' ninnlvr. l lu* pioross oan tlu'ii bo lopoaloil lor a lan.no ol'ilolax 
lalos lo j’.onoralo iho oonliibulion lo llu* nuulal l’SI> (al Iho soloolotl Iroipionoy'l of \va\os 
willi olhor oiivinnloionlial nuulo ninnbors Kopolilion ol' llu* proooss al olhor lroi|uoiuios 
loads lo a ooniplolo ilol'inilion ol iho nunlal I’SH. 

I bis inolhoil lor Iho dolorniinalion ol‘ raiuloinl> )*,onoraU*il niodo sirnoluios is in laol ipnto 
nonoial. aiul (ho molhoil oan bo applioil oijualh woll to oohoroni soinul lioKIs l or (hm annu 
lar ilnols. Iho spin lalo oan bo rolaloil lo Iho propagation anj’lo ol Iho inodo. IT (ho molhoil 
IS ai'plioil (o a oohoroni soinul liohl uhoso froipionov spooinnn oonsisisol'a rurulamonlal 
limo aiul ils harmonios. applioation of tho proooiinri' pormils Iho nioilal onoi(',\ of all moitos 
spinniii!’. out a o.ivon ralo to bo ilolorminoil by simply piooossinc. iho imliltoioil sij'aials in 
tho niannor osoriboil proMonsly. 

Whon apphoil lo thoo ilinionstonal ilnols, iho molhoil will ilolorinino only Iho oiromnloronlial 
inoilo siiiionno at iho vliiol wall, so in ils prosoni form il oan I'ivo no info:*nalion about iho 
i.tibal inoilal struotuto. 
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S.o CONCLUSIONS 


1 I'hc method has been assessed experimentally for sample modal structures and is reas- 
onably accurate when the number of available microphone measurements is comparable 
to the total number of propagating modes in the duct, b'or the experiments in tins 
program, up to ten microphones were used for the modal calculation. I he accuracy 
deteriorates as tlic number and average amplitude of those modes not included in the 
modal calculation increases. 


2 A procedure has been established to estimate for a specific potential application, the 
number of microphone systems that would be needed and the acoustic pressure itteas* 
urement and microphone location accuracy that would be required to achieve a desired 
accuracy in the calculation of the amplitudes and phases of the dominant propagating 

modes. 


3. Two approaches have been defined for measuring the modal PSD of a random sound 
field in two dimensional annular ducts; one requires cross-spectral analysis of the micro 
phone signals, whereas the other depends on processing through the use of incremental 
delays, the outputs from a circumferential array of microphones whose signals have 
been narrow-band filtered. 


’’’inraimG r\,-,R not bilmb» 
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6.0 ri:(ommi:ni)ations 


Hiiscd on till' work comiuctcil umicr tliis coiilraft, the following recomniemlalions arc sug- 
gested (or contimnng (lie development of methods for determining fan sound mode slriie- 
tures: 


I he metliod lor determining eolierent fan sound mode struetiires should be evaluated 
on a high speed fan rig, such as the NASA 50 cm CO in.) fan rig. The fan design should 
be .such that a strong rolor-stator interaction tone noise is generated, and provision 
should be made tor minimi/ing the levels of intlow turbulence and steady distortion. 

In this way. the amplitudes of extraneous modes not included in the modal ealeiilation 
will be minimized. 


2. Develop and experimentally assess circumferential microphone array methods that can 
be used to determine the pressure amplitude at the wall to define discrete coherent 
and/or random fan duct sound fields in terms of circumferential mode orders. This 
would provide guidance in selecting the dominant modes for use in the modal calcula- 
tion program (MCP) developed under this contract. 

3. Apply the circumferential array method, together with the existing P&W A/NASA pro- 
cedure (MLP and MCP), to the problem of determining the amplitude and phase of the 
dominant modes in terms of both circumferential and radial order when a relatively 
large number (i.e., 60) of discrete coherent duct modes are present in the P&WA 10- 
inch circular duct. 


4. 


5. 


I he method for determining the amplitude frequency spectrum of coherent and/or 
randomly generated .sound modes in a tliin annular duct, using a circumferential array 
of microphones should be developed and assessed experimentally on a small-scale low- 
speed fan rig such as the Pratt & Whitney Aircraft 25 cm (10 in.) fan research rig.’ 

The microphone array method should be extended to include axial arrays so that the 

radial structure of coherent and randomly generated modes ean be determined in three- 
dimensional ducts. 
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APPENDIX A NOTATION 


Symbol 

“k 

bn^bn,(t) 

B 

Bn 

bm(t) 

Cm 

C 

^mn 

C 

'-nni 


c 


dp,dz,dx,ctc. 
dAj^^, dZ|j^,etc. 

Ex ( ) 


English Symbols 
Description 

amplititde of modal pa'ssiire 
distortion harmonics of downstream How 
random time function 

- cross correlation functions 

number of rotor blades 

- amplitude of measured pressure 

distortion harmonics of inflow 

random time function 

drag coefficient 

complex amplitude 

modal amplitude complex coefficient 

cofactor of an element in llnm 
stator vane cord 
speed of sound 

errors in indicated variables 

error in variable due to error in superscript variable 
characteristic radial function 

expected magnitude 


i:‘> 
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Symbol 


APPENDIX A NOTATION (Continued) 
Description 


l‘(r. 0 , X. t ; w) 
y) 


f„(t) 




M 

Mx 


m 

in* 

N 

N 

N 

N.,. 


1,)0 


E’ouricr tnmsforni openitor 

incoherent noise filtered at frecpieney u 

wake delleit from airfoil 

array signal average 

complex matrix giving p in terms of Z 

amplitude of elements of H 
loss faeior 

integer 

charaeteristie railial eigenvalue 

ehaiacferistic radial eigenvalue 

defined by eqn. IM5 
axial wave number 

number of calculated modes 

specific value of m 

axial Mach iuimlx*r 

circumferential mode number 

highest propagating circumferential mode number 

number ol detectors in array 

number of calculated modes 

number of rotor revolutions 

number of duet supported modes 


AP1M;NI)1X a notation (Comimietl) 


Symliol 
l’(f, 0 , X, t;cj) 
< P> 

AP 

Q 


^nin 
R (a. 

Ro 

Rg(0;or) 


S(m,u)) 

(Itol). S,;, llojl) 

S^Cto) 

S^(o.co) 

t 


li 

ti 


Description 

tomplex nicasurctl (local) pressure 

average pressure for N averages 

probability distribution ilcnsity about mean A 

pressure difference 

Inverse of H giving Z in terms of Q 

dynamic head 

amplitude of elements of Q 

cross correlation function 

real part 

cross correlation function 
partial cross correlation function 

radial coordinate 

radius of narrow annular duct 

modal power spectral density 

imtdal component psd functions 

cross power spectral density 

array signal psd 

space average crov. psd functit)n 
time coordinate 
mean inflow velocity 
stator wake deficit 


I.M 


I 


I 


Syinlv)! 


V 

V 

X 

y 

2 

2 

b-l 


APPENDIX A - NOTATION (Continual) 

Doscription 

conibiiiiitions of spectra! functions defined l)v 
cqs. .r.Tl6 

flow velocity downstream of stator 
number of stator vanes 
coordinate in direction of axis 
transverse coordinate ot stator wake 

complex modal pressure 
array factor 

magnitude of array of factor; array Response or 
array response function 


Symbol 


0 

^mn 


'j 


A 


im 


Greek Symbols 
Description 

phase of elements of H 

array factor argument. /J = (m - tor?) 

phase of elements of Q 

sum of phase angles 

determinant of a N x N matrix 

delta tunction defined by properties; 

^ Cjij for m - 0 . |, 

«A 




5 ( ) 


(Jrci'k Symbols ((’oujiimctl) 


Symbol 


n 

0 

o 

a 

‘'A- ‘’x- Op 

T 

r 


<l> 


m 


(to) 


0 

0/,\lCO). 0,,*’ (cO) 

'J',,, (UJ) = i uo> 

0;','; (CO). 0l|;‘ Ico) 
i2 


ni’si'ription 

Slim of pliaso atiplcs 
ilcriiu'il by oqii. B-1‘) 
array ilolay rate, At;A(J 
atipular coordinate, radians 
radial mode index 
luib-tivtip ratio 

aiiiudar separation of detector pair 
standard deviation of errors of subscript variables 
period ol one revolution of the rotor 
delay time 

combination of psds 

phase ol inoilal pressure 

power spectral ilensities of a,,^(0. bj^|tt) 

combination of (^o). (cj) cross psds 

phase of measured pressure 

cross psd functions of a,„(t). b,„(t) 

mode spin velocity 

radian fret|uency 


(iii’i’k S\ml)<i|s i( onitmu'd I 


liulii'ii's 
S\ mliols 
I. k. I. n 
\| 

It. m. N 
\ 
tt 

\, f. 0 

No tilt ion 
S\ llll ^lol^ 

{T 

fj 


I )i‘Miipi ion 

inii'foi 

iniiKMiin;' N|'i\‘ilk’ imuli s 
luimhi-r of 

tliinil>i.'r t»f iniiTophoiK' loc.ilioii> 

spi\-iik- miaopliotu- Uvaiion 

.ivi.il. r.iiiial. .Hill ciramifi'ivnlial I'ooiilin.n 

nov Tiplion 

optional iisapi' for (.•ohimii matrix 


optional -i^ask' tor Mjiiaro matrix 
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(HM 


0 l;K-ih alo llH. .Ici.rminalion ,.r llu- s,vdnc ollKIs o„ mo.k- a,ul plinso 

I, N iik'»mv,l |.r«M,ros. ami ,H,o om.r. 

cq , nkmTlt ") wlll^''^ "“i ""'■■''T'"”''' Ilk- l«o k-nmo,, llu- RllSof 

B.l I-M-rcT Ol' TRLSSIIRH AMPLlTUDi; AN!) PHASl- hRRORS 


t'onsiiior 


c/if '■ 


(U:) 


.Ik 1 pik’.NMiik' P, ;it (lu‘ n mtk‘iophom> lokcilion iuis ampliimic B , ;nui i>liw tmil,. 

„ k-q. „kk.„o,„k.,„ o,rq. an., an., non,: / 

Icli iKiiul sikik- Ol k'kiiKition (B:) can bo oxpi-ossod in tonns ol'llioso k'nois iluoimli tho' 
appropriato parl.al dkMivalivos to give llio tollowino. ^ 


Th 




(B ,0 


|.m.,l and j .l/j, I am a.knm. „,alric.-s imlkaiinp tlio in.lq...„,l,,„ pr.-,,.„rq ..innliin.l.. 

"7 -ii'iiik u-nq-K-o a, ...v,TK.i..,„s. II, o "q, 

— I has eloiiK-nts — which ivpivsciil ti.o ct'lVcl of ihc t.“' pressure amplitude 


error on the M“‘ modal eoertu iem. I Ins matrix, and the Ckimi>anion | are to be 

e’n' Ml‘u' m"’' r' IB:>. Considerinp. without loss in 

t.kiK lalilN the chaiijie m one ol the components (the M^i) of the column m mix d/ 
equations (H:t and (lUtpive: ' 


-V 

Cf ' , y ^ / _ 

*■' . kj \fn c-f n 

nu 


lIU) 


n ' C_ , ^ ■ ) f>; / 


iB.m 
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in i-’qualioii (H4) and elimiiiatiiit? (l/.[^ belweeu ciiuations, 


(H4)aiul (B.S){?ivo: 




h 


"■■/ J'®" mi ” ” %! cf'I'n 


Since Ihc ilB,j anil the ili//j^ arc all indepcmicnt variables, the only way this equation can 
be satisfied for arbitrary values of dB,^ and is if corresponding coefficients of each 
dBj^ and di/^j^ on each side of equation (B6) arc equal. 

Thus , 

}fn J</r? 

Now expressing the complex quantities and p in polar form and taking the 
derivatives required in eipiation (B7): 

From - A^,,. 


(B7) 


^ " J Ay^ id*^ \A A 

7Bn Tsn^ 




f i Am J 

ff'*' ‘ u Am hn V ' 


and (B8) 


and from p„ - B„ e*'^'*' 


13 () 


J 

l^n 
^Pn 

J ' 


^ e. -'/>? 

/ /I • 'f'n 


(IFI) 


P (Bb) 


SubstiUilinnoilu:»>‘‘’'is(H«) ;nui (H'J) in i-iiualum (H7) aiul writing 




(imi) 


i'Ainating real ami imaginary parts gives tlie following rosnlls for the four inlluencc 
ooc flic ie n I mat ri ccs : 


[ f f/i jy 

( 7 ^' / -M] 

t'fjn] I A-' 
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lik'si' ;||V loin of Ihf inHiii-iui- riiiulioiis invsi'ntikl in i't|ii;ilioii (.V.’-IS). Tlu- :vm;iiiuU‘i 
iiii* inlliK'iuv luiu lions ivLiU'il (o oirois in min'0|'lu«tu' ItuMlioii jiiut (his siihji'oi is 
Oi»iisiili'ivil lU'M. 
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I lu' I'llivi ol I’lrors niiil in ilu' loonlinnli's of llu' mii'ioplioiu' locniion 

aiv ri’pivsi'iiu-il hv lIu' sivoiul (orm on llu- KIIS iX i-tinalion iHl). 



Il slioiilil hi- noli'ii thal llu- loliiinn matrix proiliia-il hy llu- pro-JiicI ilil*/'. can he 
consiilcicii lo he similar lo a sel o! pressure errors, since ilH-/ prenuiltiplied hy 0 is 
similar to llie sel of relations lor d/ Qdp) ilial has just heen analy/ed. TiiereVore. 

(he analysis {liven in the preceediii{i section apply ei|iially well here. The only com- 
phcalioiis result from tiavin{i lo formulaic explicit expres.sions for dll*/ in.slead ofsiarliiiii 
with {lixcii values ol dp. Insteail id seekinii ilirectly expre.ssions lor Ihe partial 
derivatives ol / with rcsneci lo r„. Xj,. and du. il is more si raiiihl forward lo find d/ in 
terms ol dr„. dx,^. and dd^^ ami lo IlnalK eipiale coefllcients of Ihe differcnlials in 
Older lo obtain expre.ssions for the derivalives. 

In ohlaminti an expression of dll. two features exisi thal simplily Ihe aljiehra. l irsl. 

Ihe eiiualion s\slem (Hi:) is linear in Ihe elements of dll. so that Ihe conirihniions 
ol errors dij^, dx,,. and dd|| may lu- ohiained separately and then superpo.sed. Secondly. 

Ihe df|,, dX||, and ddn errors of positions of the n**’ microphone location affecl onl\ llie 
elements dll in Ihe n'li row, Ihe olhei elements ol the dll matrix heinj', zero. I hese features 
allow explicit expressions for the partial derivalives to lu- ohiained. 

Theconlrihulions lodH due lo dr„. dx„, and dd„ will he denoled hy 'dll,,;.. \ill ^'dlt 
where ihe .stihseripl n denotes ihal onlv Ihe u"' row of |dll„,^n is affeeled, I he suhscripl's 
repl.ices Ihe M suhseripl lo denole successive eletncnls in lliis i?l' row. This leplacemeiil 
is necessary to avoid siih.seiiuenl confusion when it heconies necessary |o ivrform summatious 
with lespecl lo an index denoline. these suecessive n*’’ row elemenis the index M is 
reserved lo specify a particular mode coefficienl. and slu>uld not aKo he employed as a 
ihimmy or summation imlex. 

Ihe column malrix dll./ due lo errors in localiou of Ihe n'** microphone only can he 
Wiitlei' out in full .is: 
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If tho compk'X ipiantitics on tin.' UHS of oqiiation-U^lk) nro then cxiivcsseil in the polar 
forms: 
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On siihstilulinn (H17) ;nul (HIS) into (Hl(>) oqnatinit n-al :nul iinaitinaiy parts, aiul I'S- 
trai'liii!’ Ilio I'oi'ftu'ii'nts of ilin. tix„ aiiii Ilto followiiip oxpivssions aiv ol>laiiu-il 

for till' intliii'ikv I'oi'ffii'ionts or I'artial ilorivalivi's: 
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APPUNDIX C 

THF LI FUCT OF NONUNIKORM INFLOW ON ROTOR-IGV 
INTERACTION TONE NOISE 

The duct sound field gcnenited by rotor-KiV interaction consists of a set of duet modes 
that can be predicted by an application of the theory presented by Tyler & Sofrin 
(ref. I). Sofrin & McCann (ref. 2) extended this theory to include the effects of uniform 
axial flow. One of the main results frotn this work is that for a U bladed rotor and V 
inlet guide vanes, a set of circumferential acoustic duct modes at blade passing fre- 
quency harmonies are generated given by the expression m = nB + kV. In this 
expression, n denotes the particular harmonic of blade passing frequency and kV 
denotes the harmonics of inflow distortion due to the wakes of the IGV’s interacting 
witli the rotor. The range of k and the set of radial orders for each circumferential 
duct mode are bounded because at a given frequency, the duct can support only a finite 
set of propagating acoustic duct modes. 

The purpose of this appendix is to demonstrate how the kinematics of fan noise 
generation due to inlet-guide-vane wake-rotor interaction can change when steady 
non-uniform fan inflow modulates the vane wakes, resulting in the generation of 
additional acoustic duct modes. It is sufficient to consider the effect of steady 
azimuthally non-uniform inflow to a fan witii high hub-tij) ratio where no high order 
radial modes arc generated. 

It the inllow upstream of the IGV’s is given by U (0), the flow distortion downstream 
of tfie rotor will be of the same general shape but with vi.seous wake deficits super- 
imposed as shown schematically in I-igure ('I. Since tiie viscous wake ilefieit from 
each vane dejicnds on the free stream velocity in the vicinity of that vane, the wake 
deficit from each vane will be different and can be denoted by Uj where J 

denotes thej^’' KJV, and is the free stream velocity in the vicinity of the 

vane. 

If is a slowly varying function of 0 (relative to a vane gap), t/j can be taken 

to be a constant over the wake wiilth from each vane anil the non-uniform flow ilown- 
stream of the ICJV’s can be written as: 

VM= U ) 

j‘l 


Again, since H (0) is a.ssumed slowly varying, the wake deficit Uj(«i«(Jcan be reprc.sonled 
by empircally derived formulae for wake deficits behind isolated airfoils. An example of such 
formula is that given by Silvcrstein et al. (ref. 14). I lie form of these formulae are 
similar and the wake deficitnC^') can be written ('.enerally as 
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wlk'li' \ IS till' li.iiisv I’isi- kooiiliii.ik' i\'iU»'iikt .i( dll' miil|>i'iiU nl dii' w.iki', \ is ihr 
ilisUiU'i' ilnwnslio.iin t'l Mu’ UJ\ s, ^ is ii^- i iipiii oi I'.u h \ .iiu', Cp is du' ili.ii’, 

I'l'i'l lu'ii’iu loi I'iU'li \ .nu’. .itui t IS du' iiu'.m tlou lUi't tlii' .iiilnil. H.isi'd lUi Mu* .ilvui’ 
.issiiiui'doiis ilk' \i'lik'il\ ili'lii'il liik' In llk'j*^ \,iik‘ m dk’ loniiliikili’ s\s|i'in 
nl dk' dik'd I'.m I'l' w tilli’ii as 
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Now if tlic inflow distortion li (0) is cxpimdcd in tcnns of a I'cniricr series with 
period 2ir, U (d) can be written as; 

t JU> 

U - U S' 

iz -eP 


The inflow distortion interactini; with the rotor is thus 


V ) 


^ Z £ (CO) 

4,k » -*p 


When tlie distortion given by ^'(j) interacts witii the rotor, the cireuinferential duct modes 
generated are now given by 


= nI3 T KV t Jl 


where t! represents the distortion harmonies in tlic How upstream of the iClV’s. 

I'or low inflow harmonics of distortion (the ease most likely to be met in practice 
due to ground effects, stand effects, etc.) the as.sumptions u.sed to derive this result are 
good. The assumptions become weak for large values oft because this wonUl imiicale 
rapid changes of U (d) with d. If the absolute magnitude of the inflow distortion 
harmonics. Ib(i I. are monotonic ilecreasing. then the mode structure generated by rotor- 
I(1V interaction with non-tmiform inflow will consist of the primary circumferential 
mode given by (i.e.. with f=o). and successively decreasing side modes 

about the primary . given by . etc. This case is 

illustrated scitematically in I'igiire C.l. 

I bus, nonuniform inflow imposed on a slator-ndor interaetion. gives rise to side imules 
that would not exist if the inllow were unil\>rm. It can be seen that this process can give 
rise to significant pressure levels in extraneous modes, not included in the modal 
calculali<m. 
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CIRCUMFERENTIAL MODES GENERATED 
(WHERE m nB + kV + j) 
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APPl'NDIX I) 

()KDPR-Or-MA(iNITlJl)li ANALYSIS OF INFLDICNCF. COFFFICIFNTS 

Iti seel ion 3.2.4 :md Appendix B^dosed form expressions were derived for the first 
order errors in ciilcuhited mode amplitude aiul phase due to errors in the knowled^^e ol 
microphone loention and errors in acoustic pressure measurements. Tliese expressions, 
wliich are a set of partial derivatives, are termed influence eoefficients. Since these 
coefficients are calculated alonp with every mode structure calculation, it is possible to 
use the values of the eoefficients to determine the sensitivity of the method for fan 
sound mode structure determination for every modal calculation performed. 


A need exists, liowevcr, to have some estimate of the sensitivity of the method for any 
given situation before conducting a test and before running the MCP. This estimate can 
be obtained from both a consideration of an order-of-magnitude analysis of the influence 
coefficients and a review of the experience gained (and reported on in section 3.5) from 
conducting this program. The objective of this Appendix is to present this order-of- 
magnitude analysis. The analysis will be presented separately for i) pressure amplitude 
and phase measurement errors, and ii) microphone location errors. 


i) ANALYSIS OF PRESSURE AMPLITUDE AND PHASE ERRORS 

Lour influence coefficients ( ^ ) express the error in the 

calculated modal amplitudes and phases due to acoustic pressure measurement errors. 
The partial derivative, 

u ^ ^ (Di) 


obtained from e(|uations (3.2-18), wilt be used to illustrate the order-of-magnitude 
analysis. The analysis for the other three partial derivatives is similar and not shown 
here. Tile combined standard deviation for a given pressure measurement standard 
deviation, taken to be the same at all locations, while assuming all other measurements 
are perfect, can be rewritten from (3. 2-1')) as: 


_ 

^ CTf \ ^ 


where ^ deviations of the amplitude error in the calculated 

mode and acoustic pressure measurement respectively. 


It is convenient to initiate the order-of-magnitude analysis by examining the most 
probable magnitude of an elemeiit^>*^h iu the inverse matrix Q l irst, consider 
the original matrix H which can be written in full as: 
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^t), e. 


( 1 ) 3 ) 




wliere /^^^ -‘ Kxa^ //? '' Ofj and = 4?^ ( < ) arc 

consistent with the notation in Appendix A. 


Since in most cases considemd in this program, the duct microphones were located on 
the duet wall, and because of the simplification achieved in the ensuing analysis, in 
what follows, the microphones will be considered to be restricted to the duct wall. 
With this assumption, -- f L) . That is, the H-functions arc 

independent of microphone location. The determinant of the matrix 11 can then be 
expressetl as; 


c/e.-^ // 



(D4) 


where Ann is the NxN determinant of jj with the li-functions factored out. 
The elements of the detenninant are unit vectors. Similarly, the determinant of the 
co-factor transposed corresponding to the ///,// element, can be expressed as: 


Ljutn ' / - ^'7 


■M 
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(^.n) 

r\ -!.n t 


( 1 ) 5 ) 


where /\ n-/. m is the (N-1 ) x (N-1 ) determinant of the transposed co-factor of 
the nth row and mth column of the matrix II, The elements of this determinant are 
also unit vectors. 
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An dcmvnt fHn in the invo.se n.atrix jg is ilor.ncd hy the ratio ol ihe detor.nmant 
of the oo-laotor transposoil atui the determinant. Thus, usinn the iireeedm}t K-lalions 
6)4)aiul(n5)the modulus oC the element becomes, after some smiplifymg: 
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bmhining equatioiisQ)l)and(ll^ and using equationQ)b)for^/y^ ; 
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In order to obtain an order-of-magnitude analysis of equation 

<if the unit complex vectors in the two determinants AH*and a //-/.> / an i the phase 
aimles '6 ^ • denoted by 7^v • will bo assumed statistically independ- 
ent and 'uniformly distributed in £o,7.'lO ■ Also, the total number ot modes 

given by N will be assumed sufficiently large so that a reliable statistical analysis can be 
conducted. Although in fact the phase angle:, of the determinant and^^„are deter- 
niineil once a set of microphone locations has been fixed, it is conjectured that il N is 
large these angles will probably cover tlie range&#f in a reasonably uniform manner for 
most’ cases considered, nxperience trom the program conducted under this contract 
indicates that tlie above assumptions aa- good. A comparison of the order-ol-magmtude 
results from this Appendix with the exact values for the infiiiencc coellicients shown m 
I-igurcs 3.5-1 and 3.5-2, demonstrates that the order-of-magnitude results ol this section 

are representative of the exact values. 

Now. the determinant A A/, V . consisting entirely of unit complex vectors can 
be expanded formally and written in the form; 
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where J, j is assumed to be independent of n so tliat the expressions for the 

determinant can be taken out of the summation over n in equationO)"^. Substitutintt 
0)S)and(l)d)into equation (D7) and s(p)arin(> eaeli side, the expected value of the 
variance of can be written as: 
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The three expected value expressions in the RHS of equation (DIO) will be considered 
separately a'calling that all the ^ , f, and are statistically independent 

and uni formally distributed in ^ jfj . ' 
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Now, il w:is sluiwn in ApiHMulix I' Ihiit l-x ^ ^ ' {} N 

;iiul witli<7*.y , imitomily ilisirihiiiod in C^/'^TO I'luis, ilii' iti’mimimlor 

i';in Iv v’X|i;n 111 i'll, yk’lilint: lor larpo N: 






'■iV 

M- 



Vtj <1^2 




• 2 2 e. ^ 

J/X 


^ ^ }2 


J, J A/ ^ 


Jxa: 

J -^<2 


C'oiuinuini! this piocL'ss, for all N! inti'grations yields: 
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Substituting (D! H. (HI 2) & (D13) in (DIO): 
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I quatiitn (|)i n ivlali's the eNiH-eli-il variaiue ol'llu' iiuuial laleulalinn with the vaiiaitee 
assoeialeJ with llie eirtns iit aemistie pressuie ampliUule meastiivmeiil. In a similar 
manner, similar expressions ean be irblainerl for the other three inihienee eoellieients 
assoeiateri with measurement errors. I hey are: 
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where 


It' probe mieroplnmes haJ been nseil. the expressions (HI 4) aiul (1)1.^) wonhl be nuieh 
more eonipliealeil. lUrwever. lor eases where ‘ •. ..••■) :ire the same oriler-ol- 

mapnilnile. the above expressions ean be nseil to repieseni the oriler-or-mapnitmle althoiiph 
the api'it)ximatii»n is not as paioil as the ease lor just wall inonnleil mierophones. 
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Six inlliieiiee eoertieienis .34*' . ^ i- •' ) 
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express tlte error in the ealeiilateil moital amplitmles .mil pluses ihie to ertois in the 
knowleih'e i>t inieioplunie loeations. The partial Jerivalive Irom eipialion ( .V 1 S). 
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will lu’ UM’il lo illusli.iio ilu- uiili'i (<| nui'iiiitiiU' .iiKilysis. I'lu' ;iii;ilysis tor ilu- hiIut 
livi‘ mlliu'iuv is siiniLii. I iu. iMinhiiu'tl sl.iiul.iul iIi-vi.itMii lor 

Kv.ilioii r-iTois ill .ill miaiM'hotu' Uv.in.iiis whiK- .isMiiiiinj., .ill otiu'i errors ;iiv ato is 
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.SiiKslitmiiii! (Dio) in i|)|71 .nul iisiiii; (lu- ivsiill 
i\|ii.Uion ' 1.^): 
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.y.iiiii. the iiiiaoi'liotu-s will Iv iissmiivd ivslrivlal In lliv w;ill so tliai 87,5 - 
SiiiMiiiij' .nul liikiiii’. iliv i'\|vvU'd values ol' oavli side (assuiiiiii!'. as I'vtoiv dial lliv 
uiv slalislivally iiuK'ivtuU'iil and imifomih ilisliilnKcd in [;o, 
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VNwM.VSIS lOU nil- SI ANDAKI) DI VIA I ION Ol- I Mi: I KKOK IN rUI SSUKI- 
Ml ASI)KI MI:NT DUI- to I XTRANI ODS MODJ S 

ir IIU' to(;il miinivr of iiuuk's piop;i!>-atinn in :i l-m inlol Jik I is N |- and the ntimk'r of mode's 
usoii in tlu‘ modal lalailalioti is N. tium Iho mmilvi ol' cMiam-ons modes is (N p N). If llu- 
eonliihnlion to llu' tolnl measnivtl prossuiv amplilnde aiul phase at any piven loealion of the 
nth imule are A„ and ivspee(i\el\ . lor n 1. . N p. the resnilani of the total aeoustie 

Held at that point ean he written as; 
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I he ivsidtant ihie to the modes inelnded in the ealenlalion is: 
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and the resultant <lne l<i the evtr.ineons modes is: 
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riins. at the piven loealion 
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The sipiaie of the amplitude of the total sound field. It-', ean he iihtained In seiiaialiti!! the 
eoniple\ fnnelions into real and imai'iiian parts tml In Irieonomelrieal manipulation, thus. 
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I iiti.ition (I 4) thiMi ln‘oonios: j . 
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riu' .iNSumpiion is iun\ nudo ilul i\ , - \ ) is i.„C,o .,„d il.al ;,[ nnv diiot Icn-aiioi, llu- phases 
ol Hu- exli.UKVus modes e.m In- eonsideivd lo Iv umrormh distnluiled in |0. :^ | and e in 
Iv Heated :.s Malistieal xanabies. If the mnpliH.des of the extraneous modes are also eonsj- 
d.ud as st.mstie.)! vari.ibles suH, me.in A, and a iMobabiliix distribution funetion \\ 

i'll- "l t '!'■ "'.K |Vn.lom. ,m UM I,,, 

nuni ,m., .In.- u. lljn ..n.triln,,,,,,, „„„„l ni.i.k-v 1 h.„ „ 

iiiK l.ik.n "''■"iw i'i" I s i i .l|■m'l.■-Ml.lIllK■cMV.u■,l^.,ln,. 

I lio J~snii,|Mi,.n- nvd in llii. in.iK ms san Iv nisiiln-.i ,|,,n |,„ i, |n,|,„i 

o I 1 ,„l l„,i„ ,n,i,l,il i-alailali,.,, ,N - M nls„ mil Iv larnn. I Iv anH.|.ln,l,-s an, I |,|,as.-s 

.1 n-s,' moil,-, alll„.n,,h, s„Klh ,,vakin,i ,l,-vin„n,sln .an Iv .■„nsi.l..,.-,l Iv slalis- 

li.alls iml.'|vn,l.-nl laiial'l.'s ilv„. |, ivivralls in, siiiipl.' l.•a^,m lin -ivoilv .-sii i- 
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SiiKV tlu' Iasi iwo ti-mis nUr;aiatc to zero (lH\aiiso ol tlio miifonnly ilistnbiik'd phaso). 
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wIk'iv a is 11k- iiu'ail siptaiv anipliliuli- of tlk’si- oMiaiu-oiis iiunlos. 

rims. tlK- slaiulaul ili-vialioii in tlu' avamslic ampliUuk' nuMsiiiviiK-nl i-nor diu- to tlu- omiln 
btilion nl\'\traiu'ous moilos is: 
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This o\pivssion loi <i„ ilm' to cMiaiu-ous imulos is nsod m llu- mam l'od\ ot Itm u-poil ti> 
obtain oidvi-or imnmitmK' .'slimalos ol llu' ollV. t on llu- modal oaKnlalion nu'lhod ol a lai.v 
imiubcr ol i sliaiu-ous modv's. 


